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Annotation 


This  monograph  examines  frequency  methods  of 
analysis  of  modulators  used  in  optical  systems  for 
determining  the  angular  coordinates  of  moving  ob¬ 
jects  and  for  automatic  tracking.  The  feasibility 
of  modeling  optical  systems  for  determining  the 
angular  coordinates  of  moving  objects  by  means 
of  analog  computers  is  demonstrated.  A  method  is 
given  for  experimental  determination  of  the  spatial 
frequency  responses  of  modulators.  The  method  is 
based  on  the  use  of  a  laser.  The  results  of  exper¬ 
imental  and  theoretical  studies  are  covered. 

This  book  is  aimed  at  specialists  engaged  in 
the  development  of  optical  equipment.  It  may  prove 
useful  to  students  at  higher  educational  institu¬ 
tions  oriented  toward  instrument-making  specialties 

There  are  42  illustrations,  5  tables,  and  31 
bibliographic  entries. 


FOREWORD 


In  recent  years  several  important  works  have  been  published  on 
optical  information-processing  systems  which  make  extensive  use  of 
operational  problem-solving  methods  based  on  Fourier  transforms.  These 
methods  have  opened  up  the  potential  for  far-reaching  analogies  between 
optics  and  communication  theory.  Communication  systems,  like  systems 
which  produce  optical  images,  are  designed  for  the  processing  and 
transmission  of  information.  In  the  former  case  the  information  is 
temporal,  and  in  the  latter  spatial.  However,  this  difference  is 
insignificant  from  the  standpoint  of  using  common  mathematical  tech¬ 
niques.  This  is  especially  importemt  to  both  the  analysis  and  syn¬ 
thesis  of  optical  systems. 

One  of  the  most  important  elements  of  optoelectronic  designed 
to  measure  angular  coordinates  and  automatically  track  moving  objects 
is  the  unit  that  modulates  the  light  beam.  Combined  with  other  com¬ 
ponents  of  the  device,  it  converts  the  values  of  luminous  fluxes  from 
different  regions  of  the  thermal  field  being  scanned  into  correspond¬ 
ing  signals.  These  signals  are  used  to  determine  the  angular  coor¬ 
dinates  of  an  object  that  differs  from  surrounding  objects  (the  back¬ 
ground)  in  its  radiation  characteristics.  Moreover,  the  modulator 
performs  spatial  filtering  which  enhances  the  contrast  between  the 
object  and  the  background  and  which  is  a  powerful  tool  for  improving 
the  noise  immunity  and  operating  range  of  an  optoelectronic  device. 

In  this  book  we  present  the  mathematical  foundations  of  the 
theories  of  Fourier  transforms  and  spatial  filtering  and  the  basic 
theorems  and  Fourier  transforms  that  are  used  later  to  solve  various 
problems.  We  introduce  the  concept  of  the  spatial  frequency  transfer 
function  and  define  it  for  the  most  widely  used  optical  filters. 

We  also  excunine  analytical  methods  of  determining  the  frequency- 
versus-time  characteristics  of  optical  modulators  and  demonstrate  the 
feasibility  of  computer  simulation  of  optoelectronic  devices  that 
contain  modulators. 
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Theoretical  and  experimental  methods  are  presented  for  deter¬ 
mining  the  spatial  frequency  transfer  functions  of  modulators  used 
in  optical  (infrared)  systems  for  angular-coordinate  determination 
and  automatic  tracking  of  moving  objects. 

Chapter  1  was  written  by  Prof.  L.  Z.  Kriksunov,  honored  scientist 
and  engineer  of  the  Ukrainian  SSR  and  a  doctor  of  technical  sciences. 
Chapter  2  was  written  by  O.  N.  Naigovzin,  candidate  of  technical  sci¬ 
ences.  Chapter  3  was  written  by  Assoc.  Prof.  V.  I.  Mekhryakov,  can¬ 
didate  of  technical  sciences.  The  overall  arrangement  of  the  book 
and  editing  were  handled  by  L.  Z.  Kriksunov. 

The  authors  thank  G.  M.  Mosyagin,  candidate  of  technical  sciences , 
for  reviewing  the  book  and  for  his  constructive  advice,  which  was 
taken  into  consideration  in  preparing  the  book  for  press. 

The  authors  ask  that  all  comments  and  requests  be  addressed  to 
Mashinostroeniye  Publishing  House,  Bldg.  3,  Pervyy  Basmannyy  Pereulok, 
B-66,  Moscow,  USSR. 


CHAPTER  1 


THE  FUNCTION  OF  MODULATORS  AND  FREQUENCY  METHODS  OF  ANALYSIS 
1.1.  The  Function  and  Types  of  Modulators 

Optical  (infrared)  systems  for  angular-coordinate  measurement 
and  automatic  tracking  of  moving  objects  make  extensive  use  of  modu¬ 
lation  of  the  light  flux  .that  passes  through  an  optical  system.  The 
modulating  disks  or  modulators  Used  for  this  purpose  (which  also  are 
known  as  analyzing  diaphragms  in  the  foreign  literature)  are  thin 
sheets  having  some  specific  law  of  transparency  distribution  which 
are  located  in  the  focal  plane  of  an  objective. 

When  the  modulating  disk  moves  relative  to  the  optical  image  of 
an  object,  current  pulses  that  bear  information  on  the  object's  loca¬ 
tion  in  the  system's  field  of  view  periodically  occur  in  the  circuit 
of  the  radiant-energy  detector  located  behind  the  modulating  disk. 
Electronic  processing  of  these  pulses  is  used  to  identify  signals 
that  depend  on  the  angular  coordinates  of  the  moving  object.  These 
signals  are  used  for  automatic  tracking  based  on  the  angular 
coordinates. 

An  optical  system  accepts  radiation  not  only  from  an  object, 
but  also  from  the  surrounding  background.  The  radiant  characteris¬ 
tics  and  the  linear  dimensions  of  the  object  and  the  background  gen¬ 
erally  differ.  Modulation  of  the  light  flux  that  passes  through  the 
optical  system  makes  it  possible  to  discriminate  between  signals  due 
to  a  small  (point)  object  and  the  extended  background. 

Thus,  the  following  three  tasks  are  handled  by  a  modulator: 

(1)  conversion  of  light  fluxes  from  different  regions  of  the 
scanned  thermal  field  into  corresponding  signals  that  are  functions 
of  time; 
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(2)  determination  of  the  angular  coordinates  of  an  object  that 
differs  in  radiant  characteristics  from  the  surrounding  background 
by  electronic  processing  of  signals  resulting  from  modulation  of  the 
light  flux; 

(3)  spatial  filtering,  which  helps  discriminate  signals  from 
a  small  (point)  object  and  an  extended  background. 

The  modulating  disk  .performs  these  functions  together  with  other 
components  of  the  optical  system  and  of  the  entire  optoelectronic 
device. 

A  modulating  unit  usually  is  implemented  as  a  modulating  disk 
(a  modulator)  that  is  placed  in  rotation  at  a  constant  angular  ve¬ 
locity  by  an  electric  or  pneumatic  motor.  Schemes  using  a  fixed  mod¬ 
ulating  disk  relative  to  which  the  image  of  the  object  moves  are  en¬ 
countered  less  often. 

The  configuration  and  combination  of  transparent  and  opaque  ele¬ 
ments  of  the  modulating  disk  determine  the  type  of  modulation  of  the 
light  becun.  Shown  in  Fig.  1.1  are  scandard  schemes  of  the  modulating 
disks  of  infrared  systems  for  determining  the  angular  coordinates 
of  moving  objects  [14]. 

A  disk  with  two  series  of  opaque  bands  applied  with  different 
frequencies  (a)  provides  frequency  modulation  of  the  light  beaun  and 
makes  it  possible  to  determine  the  location  of  a  radiating  object 
in  some  half-plane  from  the  change  in  signal  frequency  at  the  output 
of  the  radiant-energy  detector  of  photocurrent  amplifier.  By  chang¬ 
ing  the  character  of  the  placement  of  the  modulating  strips  on  the 
disk,  it  is  possible  to  obtain  frequency  modulation  of  a  light  beam 
such  that  the  modulation  frequency  is  linearly  dependent  on  the  co¬ 
ordinate  of  the  object  in  the  plane  being  analyzed  (b) ,  or  such  that 
the  frequency  varies  according  to  a  meander  law  with  arm  ratios  de¬ 
termined  by  the  present  position  of  the  radiating  object  (c) . 
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Fig.  1.1.  Standard  schemes  of  modulating  disks:  a,  b,  c  - 
disks  that  provide  frequency  modulation  of  a  light  beam; 
d,  e,  f  -  disks  designed  for  phase  or  amplitude  modulation; 
g,  h,  i  -  disks  that  provide  pulse  position  modulation; 
j,  k,  1  -  disks  for  pulse  eimplitude  modulation. 

A  modulating  disk  with  a  spiral  cutout  (d)  successively  scans 
the  image  of  a  given  field.  If,  when  the  disk  is  rotated  by  some 
angle,  the  cutout  intersects  the  image  of  an  object,  then  a  current 
pulse  appears  at  the  output  of  the  radiant-energy  detector  located 
behind  the  modulating  disk.  The  position  (phase)  of  this  pulse  rela¬ 
tive  to  the  selected  origin  of  coordinates  characterizes  the  coordi¬ 
nates  of  the  image  of  the  radiating  object  in  the  image  plane  of  the 
field  being  scanned  (phase  modulation) . 

A  radial  sector  disk  (e)  produces  amplitude  modulation  of  a  light 
beam.  The  percentage  modulation  is  the  parameter  that  determines  an 


object's  position  in  the  plane  being  analyzed.  A  disk  with  spiral- 
type  modulating  strips  (f)  performs  an  analogous  function. 

During  one  half-rotation  a  modulating  disk  with  a  sector  cutout 
(g)  transmits  radiation  from  the  object  to  the  radiant-energy  detec¬ 
tor,  and  during  the  other  half-rotation  it  obscures  it.  The  phase 
of  the  photocurrent  pulse  depends  on  the  object's  location  in  some 
quadrant  of  the  field  of  view.  If  the  cutout  in  the  disk  is  made 
shaped  (h) ,  then  the  pulse  length  moreover  turns  out  to  be  propor¬ 
tional  to  the  displacement  of  the  object  (and  its  image  in  the  focal 
plane)  from  the  position  adopted  as  the  zero  position.  A  variable 
pulse  length  also  can  be  obtained  by  using  a  disk  (i)  that  has  a  cer¬ 
tain  number  of  windows  with  parallel  transparent  and  opaque  modulat¬ 
ing  strips. 

The  most  widely  used  disks  in  modern  optical  systems  for  deter¬ 
mining  the  angular  coordinates  of  moving  objects  are  modulating  disks 
that  produce  light-beam  modulation  in  which  the  signal  at  the  output 
of  an  electronic  circuit  is  sinusoidal.  The  signal  amplitude  is  pro¬ 
portional  to  the  movement  of  the  object's  image  from  the  zero  point, 
and  the  phase  determines  the  direction  of  motion  (the  plane  of  dis¬ 
placement)  .  Such  disks  are  shown  in  Figs,  l.lj,  1.1k,  and  1.11.  In 
particular,  the  transparent  elements  of  disk  j  are  arranged  in  a  chec¬ 
kerboard  pattern,  and  their  dimensions  are  equal  to  the  image  of  a 
point  emitter  in  the  plane  of  analysis  of  the  radiating  field.  As 
a  result,  the  light  flux  from  the  object  is  modulated  and  a  signal 
is  transmitted  from  the  output  of  the  radiant-energy  detector  by  an 
electronic  circuit,  while  the  light  flux  from  the  spatially  extended 
background  is  not  modulated  and  does  not  cause  a  signal  to  appear 
at  the  circuit's  output. 

These  examples  show  how  important  a  modulator  is  in  forming  the 
working  signal  of  an  optical  system  designed  to  determine  the  angular 
coordinates  of  moving  objects.  Hence  particular  attention  is  paid  to 
the  modulator  in  the  development  of  new  systems,  as  this  selection 


largely  predetermines  the  block  diagram  and  specifications  of  the 
system. 

In  recent  years  the  frequency-versus-time  and  spatial  frequency 
responses  of  modulators,  based  on  Fourier  transforms,  have  been  used 
for  analysis  and  synthesis.  The  use  of  multidimensional  Fourier  trans¬ 
forms  makes  possible  the  solution  of  a  number  of  important  problems 
that  arise  in  the  design  of  optoelectronic  systems.  These  include; 

— the  determination  of  the  characteristics  of  signals  when  they 
are  converted  by  different  optical  elements; 

— the  synthesis  of  optical  filters  that  filter  legitimate  signals 
in  the  presence  of  noise  signals; 

— the  determination  of  the  characteristics  of  electrical  signals 
at  the  output  of  the  optoelectronic  system  when  it  converts  optical 
signals  into  electrical  signals,  and  so  forth. 

Before  moving  on  to  the  solution  of  some  of  the  problems  listed, 
let  us  examine  the  block  diagram  of  a  standard  optoelectronic  device 
designed  to  determine  the  angular  coordinates  of  radiating  objects 
(Fig.  1.2).  The  input  element  of  the  device  is  a  lens-  or  mirror- 
type  objective  in  whose  focal  plane  there  is  a  rotating  modulating 
disk  6.  If  we  assume  that  the  modulating  disk  is  implemented  after 
the  diagram  in  Fig.  1.1k,  then  a  sinusoidal  signal  with  a  period 
T=(27r)/n,  where  n  is  the  angular  rotational  velocity  of  the  disk, 
periodically  appears  in  the  circuit  of  the  radiant-energy  detector 
located  behind  the  modulating  disk. 

This  signal  carries  a  "high-frequency  filler"  that  is  a  sequence 
of  square  current  pulses  t^=ir/u)  long  with  a  period  t=  {2t\)  / uj ,  where  co 
is  the  modulation  frequency  of  the  light  beam,  with  uj=zQ,  where  z  is 
the  number  of  transparent  (or  opaque)  sectors  of  the  modulating 
disk. 
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Fig.  1.2.  Block  diagram  of  an  optical 
system  for  determining  the  angular 
coordinates  of  radiating  objects: 

(1)  optical  axis  of  objective;  (2)  di¬ 
rection  to  object;  <p  -  error  angle; 

(3)  objective;  (4)  image  of  object; 

(5)  focal  plane  of  objective;  (6)  mod¬ 
ulating  disk;  (7)  radiant-energy  de¬ 
tector;  (8)  preamplifier;  (9)  tuned 
amplifier;  (10)  detector;  (11)  error 
signal  amplifier;  (12)  phase  detector 
and  filter;  Up  -  error  voltage;  U^jj  - 
reference  voltages;  (13)  reference 
generator  (FOH) ;  (14)  drive  motor  of 
modulating  disk  and  reference 
generator . 


The  signal  from  the  output  of  radiant-energy  detector  7  is  ampli¬ 
fied  by  electronic  precunplif ier  8  and  tuned  amplifier  9,  which  is 
tuned  to  a  frequency  u).  Then  the  signal  is  detected  and  the  frequency 
envelope  Q  is  isolated.  This  is  a  sinusoidal  voltage  with  an  amplitude 
proportional  to  the  error  angle  (the  angle  between  the  optical  axis 
of  the  objective  and  the  direction  to  the  object)  and  the  phase,  which 
is  determined  by  the  plane  in  which  the  object  is  located. 


Reference  generator  (FOB)  13  and  phase  detector  (<t>fl)  12  are  used 
to  determine  signals  proportional  to  the  components  of  the  error  angle 
in  two  mutually  perpendicular  planes.  The  rotor  of  the  reference 
generator  rotates  synchronously  and  in  phase  with  the  modulating  disk 
(one  motor) ;  two  sinusoidally  varying  voltages  (of  frequency  fi)  that 
are  phase-shifted  by  an  angle  of  2/tt  are  picked  up  from  the  terminals 


of  the  reference  generator.  These  voltages  are  supplied  to  the  phase 
detector  and  are  the  reference  voltages,  since  the  phase  of  the  error 
voltage  from  the  output  of  the  photcurrent  amplifier,  which  also  is 
fed  to  the  phase  detector,  is  determined  with  respect  to  them.  DC 
voltages  proportional  to  the  components  of  the  error  angle  in  two 
mutually  perpendicular  planes  are  picked  up  from  the  filter  of  the 
phase  detector. 

Basic  data  for  9ome -foreign  optoelectronic  devices  are  presented 
in  Table  1.1  [31] . 

1.2.  The  Fourier  Series  and  Integral 

Any  periodic  function  f(x)  of  a  variable  x  of  the  period  T  which 
satisfies  the  Dirichlet  conditions  (i.e.,  the  function  is  bounded, 
piecewise  continuous,  and  has  a  finite  number  of  extrema  over  the  per¬ 
iod)  is  known  to  be  representable  as  a  sum  of  trigonometric  functions: 

/W=<**+ V(a»cos2n*y-f-^*sin2«*-|-J  ,  (1.1) 


where 


1 

J  /(x)dx; 
~T 


(*=  I.  2,  3, 


2nk—dx 

T 

...)r 


(1.2) 


(1.3) 


This  series  is  called  the  Fourier  series  of  the  function  f(x); 
the  quantity  expresses  the  average  value  of  the  function  over  the 
period  and  is  called  the  constant  component;  the  individual  terms 
are  called  harmonics  and  have  frequencies  that  are  multiples  of  the 
fundamental  frequency  1/T. 
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Table  1.1.  Basic  Data  on  Foreign  Optoelectronic  Devices  with 

Modulators 


„  (a) 


F.M-2 
Bujrd  — 
Aiomc 

Reistone 
.\rsen.  1 


Aeroj.'i 

Qoneri  1 

Bendix — 
Scyscrnpcr 


Omb  H. 
Heidelberg 

T-oIA 

Aerojet 

General 

;mio 

Barnes 


rr  (b)  , 

no  xnMUCimii 


Moxv^m* 

PJiouinA 

AMCK 

(n) 

r.lOCKO« 

aepKA^io 

r 

repnaiHe* 
aaa  ojit* 
CTMlia 

UOAVJIMo 

pvICtUMA 

"ahck 

pytotunA 

.IKCK 

(k) 

Moty.iM- 

pytoutHA 

*  .tKCK 

Moiy.iH« 

pyiouiHA 

INCK 


13^ 

Merp 

aaap* 

TypM 

MM 


305 

127 

S9 


254 

254 


166 


153 


SI 

SI 


(d) 

Hoa* 

spasHfl 


(e) 

Ttm 

n;ia 


30* 


1*,3 

S*,7S 


30» 

6' 


30'— 1* 


S  Mf«« 

(bb) 


PbS 

[PbTc  (ou.) 
(r) 

PbS 


PbS  (ou.) 
InSb 

PbS 
a«£u 
(ou 


PbS 
PbSc  (o: 

PbS 


)(I) 

JV 


PbS 

PbS 


um  aoia 


O&iy^M. 

g>«*Tb. 

ToonoT- 

CTByioau 

lavMan 

Bi.'cm* 


0,S-2,7 

2.3-i.5 


0,6— 2,7 

2- 3 

3- S 

1.2— 2.7 
t.S— 14 

1,7— 2.7 

3-S 

0.7— 2.7 

l.»-2,7 

l.S>2.7 


10-U 

S-Ifl-ll 


1.7.10~*> 

7- 10-12 
7-10— » 

1.6.10-® 

7,4-10-*0 


i.s-iori». 

$.10— w 

4-10-12 

1,5-10-1® 

2,0-10“l> 


TIT 

nojoca 

•taCTOTl 

ru 


1.2 

10 


r^l^accarpeaOKicaa«  0« 
naacTHRA  m  ooaeMi 
aMRsa 


2  (y) 


1  ( 

500 

500 


(i) 

Oirrnceui  chctchi 


CoocoS  <buoTpaiuiii 
a  BpocTpaacTaaaaaa 

Hacma 
I'paa 


(j) 


n  (1) 

napamRoa  napafioaR- 
aecjroe  aapiraao.  aro* 
pRMioa  naocKoe  aepKa- 
ao  M  noaaaaa  annsa 


Ka> 


tonTpanacKan 


Kaccerpeaoaacaa 

noaeaafl  aaaaa 


KaeeerpeaoacKaa 
ROaeaaa  .iMnaa 


y )  Kaecerpeaoaeicaff 
noaeaaa  anasa 


Zlae  50  mm  .iNMaw  r 
noaeaaa  aN>aa 


To  art 
(CC) 


(m) 

npocTpaMCTBaraai.  20 


n  (PJ 
npocTpaRcracraaa. 
520 

CiMKTpaaaaan  ■  npo« 

eTpascTMRRaa 

(w) 

Cii«npa;ikBta  a  npv. 
crpaacTaeaaaa,  6003 


r 

CaeKTpaaaMa  n  npo* 
CTpaBcmaMaa.  lOOO 


aaaTpa.ta«aa  m  npo* 
cTpaaemaMaa.  I3fi0 

pajanaa  m  npo- 
CTpaaeTaeaMBb  12S 

npeerpaacTRauna 


Berpaat 

(dd) 


Key;  (a)  model,  company;  (b)  movable  element;  (c)  aperture 
diameter,  mm;  (d)  field  of  view;  (e)  type  of  PLE  [expansion 
unknown];  (f)  wavelength  interval,  microns;  (g)  irradiance; 
corresponding  to  noise,  W/cm^;  (h)  frequency  band,  Hz; 

(i)  optical  system;  (j)  filtering  method  and  spatial  fre¬ 
quency,  1/rad;  (k)  modulating  disk;  (1)  primary  parabolic 
mirror,  secondary  plane  mirror,  and  field  lens;  (m)  spatial, 
20;  (n)  plane  mirror;  (o)  catodioptric;  (p)  spatial,  250; 

(q)  germanium  sheet;  (r)  cooled;  (s)  Cassegrain,  Ge  sheet, 
and  field  lens;  (t)  spectral  and  spatial;  (u)  mirror,  mod¬ 
ulating  disk;  (v)  Cassegrain  and  field  lens;  (w)  spectral  and 
spatial,  6000;  (x)  spectral  and  spatial,  1000;  (y)  spectral 
and  spatial,  1360;  (z)  two  50-mm  lenses  and  a  field  lens; 

(aa)  spectral  and  spatial,  125;  (bb)  mrad;  (cc)  the  same; 

(dd)  spatial.  [Table  continued  on  p.  9.] 


If  the  function  f(x)  is  even,  i.e.,  if  f(-x)=-f(x),  then 


r 

T 

a,=— p  ^  f{x)co5  2nkydx\ 
« 


For  an  odd  function  [f (x) =-f (-x) ] ,  we  have 


Table  1,  continued 


(a) 

(b) 

noiBNJKmA 

E&ri 

(d) 

nojie 

(e) 

Tka 

JIUBatOB 

^OOiiyata* 
i)  mocih* 
coener* 

(i) 

■i'll}.  Ma 

i.icHen 

wep- 

9p«m 

n.n3 

AIMR  90»n 

RtCXOT 

OimMccKai  cHCTeHa 

rypu 

MKM 

ttayMBM 

Pa 

Bt/cm* 

(k) 

(ee) 

.'Mil 

Motyaa- 

25 

I® 

FbS 

1.8— 2.7 

500 

HM  25  HM  .1MRSU  n 

Bornes 

pyiotuHA 

o,*s-o.e 

noJieaaii  Aunaa 

^MCK 

:i-iic 

Mo  CV.TM- 

30S 

PbS 

1,8— 2,7 

1.10-13 

20 

KXecerpejioB- 

RoriH'S 

pviomaft 

CKAI  H  BOXBMX  .IHRSI 

AHCK 

(r) 

3.10—13 

ATR-l 

140 

12  Hpt« 

PbS  (ou.) 

2-3 

0.6  ( 

(7)KaccerpeiipiicRaB  h 

E'f'  if'inics 

(k) 

(bb) 

InSb  (oxa.) 
■  (r) 

3-5% 

2.10—13 

DOieMR  .iMBSa 

TGA-1 

M0.1VJH- 

76 

70  lipax 

PbS 

2.1-2,7 

1.10—10 

(j 

iTT 

pyfomifft 

(bb) 

ONl.IO 

180 

10*X2(P 

PbS 

1. 8-2.6 

9-10—10 

5 

(kk) 

PeibpaKHNonao 

Ji  •riii'S 

pymuiHA 

84 

2* 

0,4S->0,6 

910— 11 

MaKcyTOBCxaH 

AMCK 

(r) 

PbS  (oxji.) 

204 

2* 

0,7— 2.7 

1.5.10— 13 

?5 

Kaccerpaiioa‘ 

y<  ‘fncs 

pyiOiUMA 

nn'»3y 

0.35-0.f5 

2-10— IS 

cxaa  N  no.itaaa  .‘wasa 

30* 

PSSe  (^x^) 

2, 7-6.3 

3- 10— 11 

CllQe06  ^NATptnilH 
a  apocTpMncTMmnaM 
(A\  «facTOTa 


(ff) 

npocTpancTMiniaR. 
aajrxaaeTMa.  <>2 


Cmmrf^a.ibRiff  m  np^ 
CTpaiicTacfiHaa.  70  (140) 

(ii) 

Cneirrpa.ibHaa.  Mrno* 
BtKNoa  aoae  aptBRa 

6  MptA 


CTpaacTMHHaa.  357 
(11) 

npocrpaNCTNimaa. 
TayxuwTNaa.  210 

TpaHcraemiaa, 
rpaxuMrtcaa.  270  w  1090 


Key;  (a)  model,  company;  (b)  movable  element;  (c)  aperture 
diameter,  mm;  (d)  field  of  view;  (e)  type  of  PLE  [expansion 
unknown];  (f)  wavelength  interval,  microns;  (g)  irradiance 
corresponding  to  noise,  W/cm*;  (h)  frequency  band,  fiz; 

(i)  optical  system;  (j)  filtering  method  and  spatial  fre¬ 
quency,  1/rad;  (k)  modulating  disk;  (ee)  two  25-mm  lenses 
and  a  field  lens;  (ff)  spatial,  two-color,  62;  (gg)  300-mm 
Cassegrain  and  a  field  lens;  (hh)  spectral  and  spatial,  70 
(140) ;  (r)  cooled;  (v)  Cassegrain  and  a  field  lens; 

(bb)  mrad;  (ii)  spectral,  instantaneous  field  of  view  6  mrad; 
(jj)  spectral  and  spatial,  357;  (kk)  Maksutov  refraction 
system;  (11)  spatial,  two-color,  210;  (mm)  or;  (nn)  photo¬ 
multiplier  [FEU];  (oo)  200-mm  Cassegrain  and  a  field  lens; 
(PP)  spatial,  three-color,  270  and  1080. 


«»=0; 

r 

T 

j*/  (jc)  sin  2n*  -y  dx. 

Let  us  consider  as  an  example  the  expansion  into  a  Fourier  series 
of  a  rectangular  sawtooth  function  f(x)  that  assumes  values  of  "+1" 
or  "-1"  with  a  period  x  (Fig.  1.3a). 


We  find 
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Fig.  1-3.  Expansion' of  a  rec¬ 
tangular  sawtooth  function  f{x) 
into  a  Fourier  series:  (a)  f(x); 
(b)  first  harmonic;  (c) ,  (d) , 
and  (e)  -  curves  calculated  from 
the  first  two,  first  three,  and 
first  four  terms  of  the  expansion. 


* 

2 

j*  /{x)cos2nk-^dx=0. 


^  /  {x)sm2nk-^dx= 


10,  if  *=2rt  (even) 


=  4 


,  if  *=2rt-|-l  (odd) 


Thus,  the  function  f (x)  can  be  represented  by  the  following 


/  W- J_fsin  -5«.+-L  sin  3 

»  V  T  3  t 

+4-  sin  5 4.  sin  7  ^4-. .  .'i . 


series: 


As  the  number  of  terms  in  the  series  increases,  the  representation 
of  the  function  becomes  increasingly  accurate,  as  we  can  see  from 
Fig.  1.3,  which  shows  curves  corresponding  to  the  first  term  of  the 
series  (b)  ,  to  the  sum  of  the  first  two  terms  (c)  ,  to  the  stun  of  the 
first  three  terms  (d) ,  and  to  the  sum  of  the  first  four  terms  (e) . 

A  Fourier  series  also  may  be  written  in  a  different  form; 


m 

•/(•*)= flo + ^  cos 


where 


(1.5) 

?»=arc1g-^.  (1.6) 

The  aggregate  of  the  amplitudes  of  the  frequency  components 
is  called  the  amplitude  spectrum,  and  the  aggregate  of  the  phases 
is  called  the  phase  spectrum.  In  most  cases  the  values  of  Cj^  are 
only  real  or  only  imaginary.  This  makes  it  possible  to  describe  the 
function  as  just  a  single  spectrum. 


The  spectrum  of  a  periodic  function  can  be  represented  graph- 

2  ttIc 

ically  in  a  coordinate  system  (Cj^;  w=±  (— ip") )- .  The  amplitudes  of  indi¬ 
vidual  frequency  components  commonly  are  marked  by  vertical  lines 
of  the  appropriate  length.  Thus,  the  spectrum  of  a  periodic  function 
is  discrete  (or,  as  we  say,  it  is  a  line  spectrum,  by  analogy  with  an 

optical  spectrum) .  The  individual  components  of  a  discrete  spectrum 

2it 

are  separated  from  each  other  by  — The  spectrum  is  symmetric  rel¬ 
ative  to  the  vertical  axis,  which  passes  through  the  origin. 

Let  us  examine  a  frequently  encountered  complex  formula  for  writ¬ 
ing  a  Fourier  series.  Using  Euler's  formulas,  we  can  represent  a 
Fourier  series  as 
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where 


(1.8j 


(k=±l.  ±2,  ±3,...).  (1.9) 


The  quantity  Cj^  is  called  the  complex  amplitude  of  the  k-th  harmonic, 
and  the  aggregate  of  the  quantities  Cj^  is  called  the  complex  spectrum 
of  the  function  f  (x)  .  The  quantity  Cj^,  which  is  the  modulus  of  Cj^, 
is  called  the  amplitude  of  the  k-th  harmonic;  as  noted  above,  the 
aggregates  of  the  quantities  Cj^  amd  are  called  the  amplitude  spec¬ 
trum  and  the  phase  spectrum. 


/(*) 

ilrfm  rri  r. 


z 


— 

— 

3 

44 

j  - 

0  n  .Til 


Fig.  1.4.  A  rectangular 
periodic  function. 


Fig.  1.5.  Graph  of  the 
function  (sin  z)/z. 


As  an  example,  let  us  expand  the  periodic  function  shown  in  Fig. 
1.4,  which  has  a  width  6  amd  a  repetition  period  T,  into  an  exponen¬ 
tial  Fourier  series.  On  a  one-period  interval  this  function  may  be 
written  as 


/W= 


[a  when  — Y  • 

0  when 

2  2 


6  5 

Selecting  limits  of  integration  from  to  T“2'  find 
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A  square-bracketed  function  of  the  type  (sin  z) /z  is  shown  in 
Fig.  1.5.  It  oscillates  with  a  period  2tt  ,  decreasing  in  amplitude 
with  increasing  z  and  passing  through  zero  at  the  points  z=±Tr,  ±2n, 
±3tt,  etc. 


Thus, 


(l.lOa) 


OQ 


(1.10b) 


It  follows  from  formula  (1.10a)  that  c^^  is  real;  therefore,  a  single 
spectrum  (the  amplitude  spectrum)  suffices  for  a  frequency  represen¬ 
tation  of  the  function  f(x).  Since  k  can  assume  values  ±1,  ±2,  ±3, 
etc.,  the  spectrum  found  is  a  discrete  function  that  only  exists  at 
frequencies  w-O,  etc.  with  the  corresponding  amplitudes: 


j.  ..  »m—  sin  — — 
,41  r  r 


Shown  in  Fig.  1.6  are  the  spectra  for  the  following  values  of 


T  and  6: 


(a)  f  ^~20' 


(b)  T=|,  (c)  T=l,  S—. 


It  follows  from  an  inspection  of  Fig.  1.6  that  the  frequency  —  de¬ 
creases  with  increasing  period  T.  As  a  result,  the  number  of  frequency 
components  per  unit  frequency  increases.  In  other  words,  as  the  period 
T  increases,  the  spectral  lines  converge  but  the  amplitudes  of  the 
frequency  components  decrease.  The  shape  of  the  frequency  spectrum 
remains  unchanged,  i.e.,  the  envelope  of  the  spectrum  depends  solely 
on  the  pulse  shape,  but  not  on  the  repetition  period  T. 


■5 


111 


Fig.  1.6.  Spectra  of  a  rec¬ 
tangular  periodic  function  for 
various  repetition  periods. 


Expressions  for  the  Fourier  series  may  be  written  analogously 
for  a  function  of  several  variables.  For  example,  if  the  function 
f(x,  y,  z)  satisfies  the  periodicity  criterion,  i.e., 

f{x,y,z)=f{x-\-kTjp  y+mr„  z-\-nT,), 

where  T^,  T^,  and  are  the  signal  repetition  periods,  then  the  Fou 
rier  series  assumes  the  form 

fix,  if, 

The  coefficients  of  the  series  are  determined  by  the  formulas 


-1 


f  ffi  =  1.  2,  3,  ,  .  .  \ 
\«  / 


The  Fourier  series  is  written  in  complex  form  as  follows: 


where 


o  ,  ^r.m  ,  \ 


(1. 13) 


T  j  J 


“T  T  “T 


X«  ^  *  "  *  dxdydz. 


(1.14) 


The  real  coefficients  and  of  series  (1.11)  are  related 

to  the  complex  coefficients  c,  of  series  (1.13)  by  the  equation 


(1.15) 


A  Fourier  expansion  also  can  be  applied  to  nonperiodic  functions, 
Any  one-dimensional  nonperiodic  function  f (x)  that  satisfies  the 
Dirichlet  conditions  and  that  is  integrable  over  an  infinite  range 
can  be  represented  as  the  sum  of  an  infinite  set  of  harmonic 
components; 

V- 

/{^)=^c,e  =  , 


-iSictv  >3**4- 


dx  e  (j  jg) 
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Replacing  the  repetition  period  of  the  function  by  the  first 
harmonic  frequency  (u^=l/T) ,  we  have 


/W=2  (1.17) 


-«  r 

L-t 


Upon  passage  to  the  limit,  when  T-+<»,  the  discrete  values  of  the 
frequencies  ku  are  converted  to  the  present  frequencies  u,  and  the 
first  harmonic  frequency  is  converted  to  a  differential;  hence. 


)=  J/(Jc)*(  — 2ni«)dxj»(2aayc)d«,  (1.18) 


where  (Jf). 


The  inner  integral  of  expression  (1.18) 

?(/«)=  J/(x)*(-2nox)dx 


(1. 19) 


is  the  direct  Fourier  transform  of  the  function  f(x),  designated 
<I>[f(x)],  and  the  outer  integral 


/(•«)=  ^  g(lu)%{2nax)du 


(1.20) 


is  the  inverse  Fourier  transform,  symbolized  $  "^[fCx)].  Let  us  note 
that  the  mathematical  operations  of  the  direct  and  inverse  transfor¬ 
mations  differ  only  in  the  sign  of  the  exponent  in  the  integrand. 

The  difference  between  a  Fourier  series  and  a  Fourier  integral 
is  that  a  Fourier  series  is  a  periodic  function  in  the  form  of  a  sum 
of  an  infinite  number  of  sinusoids  with  frequencies  that  have  defined 
discrete  values,  whereas  a  Fourier  integral  is  an  aperiodic  function 
of  a  sum  of  periodic  components  with  a  continuous  sequence  of 
frequencies. 
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The  fiinction  g(iu)  is  called  the  complex  spectrum  of  an  aperiodic 
function;  it  may  be  written  as  the  sum  of  two  functions; 

g (iu) =A(u) +iB (u)  (1.21) 

or  as  the  product 

g(iu)=g{u)e^’^‘'^^  (1.22) 

where  A(u),  B(u),  g(u),  and  il;(u)  are  the  real,  imaginary,  amplitude, 
and  phase  characteristics  of  the  spectrum,  respectively. 

It  follows  from  the  formula 


f{x)^ 

that  the  aperiodic  function  f (x)  may  be  considered  as  the  continuous 
sum  of  exponential  functions  having  frequencies  from  -<»  to  +<».  The 
aimplitude  of  the  component  at  any  frequency  is  proportional  to  g(iu); 
hence  g(iu)  is  the  spectrum  of  the  function  f (x)  and  is  called  the 
spectral  density  function.  Let  us  note  that  the  spectrum  of  an  aper¬ 
iodic  function  is  depicted  as  a  continuous  curve,  and  not  as  discrete 
points  or  lines,  and  is  said  to  be  a  continuous  spectrum. 

The  Fourier  transform  enables  us  to  represent  a  given  signal 
in  terms  of  exponential  components.  The  function  g(iu),  the  direct 
Fourier  transform  of  a  signal  f (x) ,  characterizes  the  relative  ampli¬ 
tudes  of  different  frequency  components,  i.e.,  it  represents  the  sig¬ 
nal  f (x)  in  the  frequency  domain.  In  the  general  case  the  function 
g(iu)  is  complex,  and  graphs  of  the  amplitude  spectrum  g(u)  and  phase 
spectrum  ip(u)  are  needed  to  represent  it.  In  many  cases,  however, 
g(iu)  is  either  a  real  or  an  imaginary  function,  so  its  representation 
is  limited  to  a  single  graph  of  the  spectrum. 

It  follows  from  the  formula  for  a  Fourier  transform  that  it 

00 

exists  if  the  integral  /  f  (x)  <  (-2ttux)  dx  has  a  finite  value.  Since 
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the  modulus  of  k(-2ttux)  is  equal  to  unity,  the  condition  for  the  ex¬ 
istence  of  a  Fourier  transform  of  the  function  f(x)  is  that  the  inte¬ 


gral  /  |f(x)|dx  have  a  finite  value.  This  condition  is  sufficient 

—CO 

but  not  necessary.  Such  functions  as  the  sine,  cosine,  and  Heavi¬ 
side's  function  do  not  meet  this  criterion  and  strictly  speaking  do 
not  have  a  Fourier  transform.  However,  in  the  limit  there  also  is 
a  Fourier  transform  for  these  functions  [16,  20]. 

For  a  two-dimensional  aperiodic  function  f(x,  y)  the  Fourier 
transforms  are  expressed  by  the  following  relations: 

g  (to,  tv)  =  fj  f  (jc.  (1. 23) 

— a» 

/(X.  y)=JJ  g{iu,  (1.24) 

A  number  of  theorems  that  are  used  extensively  in  solving  various 
practical  problems  stem  from  basic  definition  (1.19)  of  a  Fourier 
transform  [3,  16,  20]. 

Theorem  of  Linearity: 

®  l«/i  (Jf. ») + */t  (Jc.  y)l = l/i  (Jf .  y)l + 

■  •  (1.25) 

i.e.,  the  Fourier  transform  of  the  sum  of  two  functions  is  equal  to 
the  sum  of  their  transforms. 

Similarity  Principle 

If  [f  (x,  y)  ]  =g(iu,  iv)  , 

*  When  only  one  limit  of  integration  is  indicated  above  or  below  the 
double  integral  symbol,  this  limit  applies  to  integration  over  both 
variables. 
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iLlJ 


then 

('-26) 

i.e.,  stretching  of  coordinates  in  the  space  domain  (x,  y)  leads  to 
contraction  of  coordinates  in  the  frequency  domain  (u,  v)  and  to  a 
change  in  the  overall  amplitude  of  the  spectrum. 

Translation  Property 

If  $[f(x,  y)]=g(iu,  iv) ,  then 

y—b)\= g{la,  iv)^{  —  2:iau  —  2:ibv),  (1.27) 

i.e.,  a  shift  of  a  function  in  the  space  domain  causes  a  linear  fre¬ 
quency  shift  in  the  frequency  domain. 

Parseval *  s  Theorem 

If  -^tflx,  y)]=g(iu,  iv)  ,  then 

OD 

JJ !/(•*•  =  lv)\*dudv.  (1.28) 

Convolution  Theorem 


If  'D[f^(x,  y)]=gj^(iu,  iv)  and  y)]=g2(iu,  iv)  ,  then 

^  y)Mjc,  y)]= 

m 

(«-«'.  v-v’)g^{a',‘tf)da'dv',  (1.29) 

i.e.,  the  Fourier  transform  of  the  product  of  two  functions 
fj^(x,  y)  f2(x,  y)  is  expressed  in  terms  of  the  transforms  g^(iu,  iv) 
and  g2(iu,  iv)  by  means  of  an  integration  operation  called  the  con¬ 
volution  of  the  functions  g^^  and  g2- 
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In  frequency  analysis  the  convolution  theorem  is  an  important 
theorem.  Let  us  recall  that  the  convolution  of  two  given  functions 
f^(t)  and  f2(t)  is  defined  by  an  integral 

/(0=  jA(i:)/.(<-T)dT, 

which  is  symbolically  written  as 

By  analogy  with  multiplication  algebra,  we  can  derive  the  following 
laws  of  convolution  algebra: 

commutative  law: 

distributive  law: 

flit)  ^  [/i(04-/,Wl=/iW  ^  Ut)-hA{t)  ^  /,(0; 

associative  law: 

A  (^)  ¥:  [/,  it)  ^  A  it)]  =  [/x  it)  ^  /,  (01  ^  A  (0. 

To  provide  a  graphic  illustration  of  the  convolution  operation, 
we  will  assvune  that  the  function  fj^(t)  has  a  rectangular  form,  while 
f2(t)  has  a  triangular  form  (Fig.  1.7a).  Let  us  find  the  convolution 
of  these  functions  graphically.  By  definition, 

A(0^A(O=  jA(t)A(^-t)rft. 

Here  t  is  an  independent  variable. 

The  functions  fj^(T)  and  f2(-T)  are  shown  in  Fig.  1.7b.  The  graph 
of  the  function  f2(t-T)  for  t=tj^  is  the  function  f2(-T)  shifted  by  an 
amount  t=tj^  in  the  positive  direction  (Fig.  1.7c).  The  convolution  of 
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Fig.  1.7.  Graphic  representation  of  the 
convolution  operation. 

the  functions  fj^(t)  and  £2^^^  ^”^1  determined  by  the  area  under 

the  curve  of  f ^ (t) f 2 (t^-r) .  This  area  is  hatched  in  Fig.  1.7d  and  is 
plotted  as  a  dashed  line  in  Fig.  1.7e.  The  values  of  the  convolution 
can  be  obtained  in  an  analogous  fashion  for  different  t.  As  a  result 
we  obtain  the  graph  shown  in  Fig.  1.7e. 

Let  us  also  note  that  the  convolution  of  a  function  f (t)  with 
a  Heaviside  function  gives  the  function  f(t)  itself,  i.e., 

/(/)^8(0=  J /(/)  5  (/-T)rfT  =  /(/). 

The  Fourier  transforms  of  some  functions  are  presented  in  Table 

1,2,  and  the  values  of  the  functions  sin  c (x) and  that 

X  z 

are  needed  to  calculate  the  Fourier  transforms  are  given  in  Tables 
1.3  and  1.4  [16,  20,  27] . 


Let  us  recall  that  the  first-order  Bessel  function  is  de¬ 

fined  by  the  relation 


m 

j'cos(^s^^  ?  — sltfw, 


and  can  be  calculated  by  using  the  series 

/,  (zl= - - -  .  .  .  , 

.  '  ^  2  18  ^  384  18432 
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Table  1.2.  Fourier  Transforms  of  Frequently 
Encountered  Functions 


Key:  (a)  function;  (b)  Fourier  transform; 
(c)  rectangular  (step)  function;  (d)  if; 

(e)  two-dimensional  rectangular  function; 

(f)  elsewhere;  (g)  circular  function. 

while  for  large  values  of  z. 
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Table  1.2,  continued 


Key;  (a)  function;  (b)  Fourier  transform; 
(c)  triangular  function;  (d)  elsewhere; 
(e)  Gaussian  function;  (f)  exponential 
function;  (g)  sign  function;  (h)  if. 


The  zero-order  Bessel  function  is  expressed  by  the  formula 


m 

/,(z)=— ^  cos(crcos  f)d<o. 
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Table  1.2,  continued 


(a) 

(b)  npeo^ajoSoMut  ‘Pypte 

(c)  /leflbma.-(fiyHMnuii  Mupaxa 
/(*)=5(x)-0  nput*0 

J  S(x)ix-1 

-oo 

/(*)♦ 

^(11) 

-/ 

-3T-2T-T  >0  T  ZT  3T  X 


XQcuHijcoudajiuiaii  qnjMHuun 
f(x)~casnix 


^  ^  ^  l!L  “ 

’t't't  t  T  t 


Key;  (a)  function;  (b)  Fourier  transform; 
(c)  Dirac  delta  function;  (d)  comb  func¬ 
tion;  (e)  cosine  function;  (f)  sine  func¬ 
tion. 

The  function  I«(2)  can  be  calculated  by  using  the  series 


/,(*)  =  .  .. 

^  4  ~  64  2304~ 


Table  1.3.  Values  of  the  Functions  sine  x 
and  sine*  x 


X 

sincjc 

sinc^jr  1 

1 

sinc(jf+l 
+  n)  ’ 

slnc-fjfT 
+  •■1) 

sine  (jr+l 
+•2.) ! 

smc-(jc  + 

1  +  2n) 

0 

1.0 

I.o 

0.0 

0.0 

0,0 

0,0 

n 

4 

0.9003 

0.810S 

—0.1801 

0.0324 

0,1001 

O.OlOO' 

n 

2 

0.6366 

0.4053 

-0.2122 

0,0450 

0,1274 

0,0162 

3 

—  .t 

4 

0.3001  ' 

0.0901 

1 

-0.1286 

0,0165 

0,0819 

0,0067 

21. (z) 

Table  1.4.  Values  of  the  Functions  — = - 


I 

2/t  (*) 

z 

Z 

2M£) 

z 

j^2/,(z)y 

0 

1.0 

1.0 

5.0 

-0.1310 

0,0172 

1.0 

0,08801 

0.7746 

6.0 

—0,0922 

0,008.5 

2.0 

0.5767 

0.3326 

7.0 

—0,0013 

.  0,0000 

3.0 

0,2260 

0.0511 

8.0 

0,0587 

0,0034 

3.832 

0.0 

0.0 

9,0 

0.0545 

0,0030 

4.0 

-0.0330 

0,0011 

10.0 

0,0087 

0,0001 

while  for  large  values  of  z 


The  functions  l^iz)  and  Iq(2)  are  interrelated  by  the  equation 

1.3.  Spatial  Frequency  Spectra  of  Radiating  Objects  and  Backgrounds 

Let  there  be  a  flat  emitter  that  has  a  radiance  distribution 
B(y,  z) .  The  variables  y  and  z  are  the  coordinates  of  some  point  of 
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of  the  object  in  a  Cartesian  coordinate  system  yOz  selected  so  that 
the  plane  yOz  is  perpendicular  to  a  specified  direction. 

For  objects  with  noncoherent  self-radiation,  by  radiance  we 
mean  the  power  per  unit  area  of  emitting  surface  in  a  solid  angle 
of  1  steradian  that  corresponds  to  a  unit  wavelength  interval  near 
the  wavelength  X . 

The  Fourier  transform  of  the  function  B(y,  z)  will  be  written  as 
b(la,  /«)= 2)it[— 2a(By-f«z)lrfyrfz;  (1.30) 

5  (y,  *)=  JJ  *(<«.  t‘o)*l2n(uy-\-vz)]dttdv.  (1.31) 

The  quantities  u  and  v,  the  inverse  periods  of  the  harmonics 
that  describe  the  distribution  of  an  object's  radiance  along  the 
directions  y  and  z,  are  called  the  spatial  frequencies  by  analogy 
with  frequency  in  the  generally  accepted  sense  of  the  word  which  we 
obtained  [sic]  if  the  function  depended  on  just  one  variable — time. 

If  the  dimension  of  the  temporal  frequency  is  1/s,  then  the  dimension 
of  the  spatial  frequency  is  a  unit  of  length  to  the  power  -1. 

Thus,  the  distribution  of  an  object's  radiance  may  be  repre¬ 
sented  as  the  superposition  of  eui  infinite  set  of  various  harmonic 
components  each  of  which  can  be  characterized  in  terms  of  amplitude, 
phase,  and  spatial  frequency.  Such  a  representation  of  a  two-dimen¬ 
sional  function  of  radiance  makes  it  possible  to  bring  into  consid¬ 
eration  the  notion  of  the  spatial  frequency  spectrum,  which  plays 
an  important  part  in  the  analysis  of  modulators. 

Let  us  examine  the  main  differences  between  spatial  and  temporal 
harmonic  oscillations.  The  frequencies  of  temporal  oscillations  may 
range  from  0  to  +“>,  since  time  cannot  be  negative;  spatial  frequencies 
can  assxime  values  from  -<*>  to  +<». 
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Fig.  1.8.  Change  in  the  frequency 
of  a  spatial  harmonic  upon  passage 
through  an  optical  element. 


Fig.  1.9.  Diagram  illustrating 
the  conversion  of  a  spatial  fre¬ 
quency  into  a  temporal  frequency 
of  changing  electrical  current. 

(1)  radiant-energy  detector; 

(2)  objective;  (3)  plane  of 
object . 


Upon  passing  through  linear  dynamic  sections,  temporal  (elec¬ 
trical)  harmonics  undergo  a  change  in  amplitude  and  phase  while  main¬ 
taining  unchanged  frequencies.  Upon  passing  through  optical  elements 
spatial  harmonics  may  undergo  a  change  in  both  amplitude  and  phase 
on  the  one  hand,  and  spatial  frequency  on  the  other.  For  example, 
if  the  magnification  of  an  optical  element  is  3,  and  the  spatial  fre¬ 
quency  at  its  input  is  u  ,  then  the  value  of  the  spatial  frequency 
at  the  output  of  the  element  will  be 


u  -u  I 

BbOC  bx8 


(1.32) 


where 


according  to  the  symbols  used  in  Fig.  1.8. 


The  spatial  frequency  u  can  be  converted  into  the  temporal  fre¬ 
quency  f*  of  an  electrical  current,  as  shovm  in  Fig.  1.9.  At  every 
instant  a  section  of  the  image  having  an  area  Aa  is  projected  by  the 
optical  system  onto  a  radiant-energy  detector  [PLE] .  The  load  re¬ 
sistance  of  the  radiant-energy  detector  will  cause  an  electrical  cur¬ 
rent  to  flow  with  a  frequency 

f*=Vu  (1.33) 

where  V  is  the  scan  rate. 

The  spatial  frequencies  of  different  regions  of  the  spatial  fre¬ 
quency  spectrum  play  different  parts  in  the  construction  of  images 
of  radiating  objects.  The  low-frequency  components  carry  primary  in¬ 
formation  on  extended  distributions  of  radiance,  on  large  details, 
and  on  uniformly  radiating  regions  of  the  object.  Medium-frequency 
components  of  the  spatial  frequency  spectrum  of  signals  are  important 
for  correct  reproduction  of  tones,  and  high-frequency  components  are 
important  in  the  construction  of  small  details,  sharp  transitions, 
and  image  contours.  If  low-  and  medium-frequency  components  are  re¬ 
produced  without  distortion  while  the  high-frequency  components  are 
attenuated,  small  details  and  sharp  borders  of  the  object  will  be 
conveyed  poorly  in  the  resultant  images. 

Returning  to  formulas  (1.30)  and  (1.31),  let  us  note  that  the 
function  b(iu,  iv)  is  a  complex  function  of  the  spatial  frequencies. 
The  aggregate  of  the  moduli  of  b(iu,  iv)  is  called  the  spatial  fre¬ 
quency  spectrum  of  the  function  B(y,  z) . 

The  function  b(iu,  iv)  is  the  spectral  density  of  the  function 
B(y,  z)  and  is  expressed  by  a  direct  Fourier  transform.  The  expres¬ 
sion  of  the  function  B(y,  z)  in  terms  of  its  spectrum  is  the  inverse 
Fourier  transform. 

The  function  B(y,  z)  often  may  be  assumed  symmetric  relative 
to  the  origin;  in  this  case,  assuming 


2=QC0S^,  y=esin(!>, 
2nB=wcos»,  2nx)=TOsina 


■5^ 


and  inserting  the  new  variables  into  Fourier  integral  (1.30),  ,  we 
obtain 


m 

b{lw,  <p)=  f  f  fl(0.  t)*l  — e®{cos? 
6  6 


X  QdQd'i. 


cos + sin  ?sin  6)]  X 

(1.34) 


We  write 


4— ?=6,  biiw,  o)=b(w)e‘'^, 

B  (c> 


then 


all 

b(w)=t-‘'  j  j  5(q)*(  — CTOCos6  4-rt6jQrfQ</9,  (1.35) 


or 


where 


ft(w)=2n  j'e^((?®)s(Q)rfo. 
6 


11.36) 


— 

/,  (qw)=-^^  g-iewco.4  +  .>i»^9 


0 

is  the  integral  representation  of  a  Bessel  function. 

Expression  (1.36)  is  called  the  Hankel  transform  of  the  function 
s(p).  The  inverse  Hankel  transform  is  defined  by  the  integral 

«• 

s(«)=  f  w/,(Qw)d(w)rfw,  (1.37) 


or 


^  (c. 


•• 

J 


wf^{Qw)b{w)dw. 


(1.38) 
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If  the  function  B(y,  z)  is  symmetric  with  respect  to  the  origin, 
then  n=0,  s(p)=B(p),  and  a  zero-order  Bessel  function  must  be  inserted 
into  the  Hankel  integral.  We  obtain 


6lTO)=2.n  J  Q/o(QV)B(Q)dQ, 


B(Q)=2n  (  c/o(cw)^  («)</». 


(1.39) 

(1.40) 


The  Hankel  transform  is  distinguished  by  the  following  property: 

*  ("t)  c/o((lw)S(aQ)dc.  (1.41) 

0 

whence  it  follows  that  when  the  function  B  is  expanded  in  space,  its 
spectrum  narrows. 

If  the  object's  radiance  varies  not  only  in  space,  but  also 
in  time  (as  is  characteristic  of  "flashing"  or  fast-moving  objects) , 
then  it  is  described  as  a  function  of  three  variables,  B(y,  z,  t) , 
for  which  a  Fourier  transform  can  be  found.  In  the  particular  case 
where  the  function  B  depends  only  on  time,  we  have 

/,(/*)=  \  (1.42) 

The  function  b(f*)  is  the  spectral  relation  that  describes  a  temporal 
spectrum. 

Now  let  us  consider  the  spatial  frequency  spectra  of  backgrounds. 
The  radiation  of  clouds  and  of  topographic  details  of  the  relief  of 
the  earth's  surface,  inhomogeneous  atmospheric  radiation,  and  reflected 
or  scattered  radiation  from  various  sources  have  the  feature  that  the 
radiation  field  they  produce  can  only  be  represented  as  a  random  func¬ 
tion  of  radiance  that  obeys  some  statistical  laws.  Therefore,  iso¬ 
lating  useful  information  on  an  object  that  is  against  a  radiating 
background  entails  major  difficulties. 
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The  full  diversity  of  backgrounds  usually  is  divided  into  "en¬ 
sembles"  (e.g.,  forest  or  steppe  landscapes,  large  urban  areas,  con¬ 
tinuous  cloud  cover,  etc.)  which  have  certain  statistical  character¬ 
istics  of  a  random  two-dimensional  radiance  function.  Moreover,  a 
random  radiance  function  is  assumed  stationary  in  many  cases  in  order 
greatly  to  simplify  the  mathematical  conversions.  For  example,  it 
has  been  established  that  the  random  function  that  describes  a  bright 
background  is  nonstationary  in  the  general  case.  However,  the  back¬ 
ground  of  a  cloudy  sky  nonetheless  may  be  assumed  steady  at  small 
angles  from  the  directions  to  the  sun  and  horizon  [18,  25]. 

A  peculiarity  of  random  functions  that  must  be  dealt  with  in 
the  design  of  optoelectronic  devices  is  the  fact  that  they  usually 
are  not  sign-alternating,  and  their  values  always  are  positive.  The 
corresponding  signals  at  the  output  of  the  radiant-energy  detector 
also  are  unipolar. 

The  set  of  statistical  characteristics  of  a  random  function  that 
is  simplest  for  mathematical  processing  includes  the  mathematical 
expectation,  the  dispersion,  and  the  correlation  function. 

The  mathematical  expectation  of  a  random  two-dimensional  radiance 
function  B(y,  z)  when  there  are  several  realizations  of  the  function 
is  the  "average"  function  about  which  specific  realizations  of  the 
random  radiance  function  vary  in  different  ways: 

m(y,  z)=M[Bj^(y,  z)  ]  (1.43) 

where  n  is  the  number  of  the  realization. 

The  dispersion  of  a  random  two-dimensional  function  characterizes 
the  spread  of  possible  realizations  of  the  random  function  with  re¬ 
spect  to  the  mathematical  expectation: 

D(y,  z)=M[B^(y,  z)-m(y,  z)  ]  ^  (1.44) 

The  correlation  function  shows  the  degree  of  interrelation  of 
the  quantities  of  a  two-dimensional  radiance  function  at  two  adjacent 
points  having  coordinates  (y,  z)  and  (y+Ay,  z+Az) : 
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^\^y<  {y-\-^y,  z-\-^z)]  —  M  {[B{y,  z)  —  m(y,  2)Jx 

X  + A//.  2-rA2) — 'n(y-{-^.y,  z-j-A^'']'.  (1.45) 

When  Ay=0  and  Az=Q ,  the  correlation  function  becomes  the  disper¬ 
sion  of  a  random  function.  Thus,  the  need  vanishes  for  the  dispersion 
as  a  separate  characteristic  of  a  random  function;  the  mathematical 
expectation  and  correlation  function  may  be  adopted  as  its  basic  char¬ 
acteristics  . 

If  the  random  function  is  stationary,  then  its  correlation  func¬ 
tion  depends  solely  on  the  quantities  Ay  and  Az  and  is  independent 
of  y  and  z. 

The  spatial  frequency  properties  of  stationary  random  functions 
are  characterized  by  a  so-called  Khinchine-Wiener  spectrum  (X-B) , 
which  is  the  dependence  of  the  spectral  density  of  the  dispersion 
of  a  random  function  on  the  spatial  frequency. 

The  Khinchine-Wiener  spectrum  of  a  random  function  B(y,  z)  of 
the  spatial  distribution  of  radiance  is  defined  as  follows.  Let 
the  random  function  B(y,  z)  be  defined  in  the  domain 

IyUA;  |z|SC 

The  Fourier  transform  of  this  function 

A  C 

b{u,  “0)=  J  dy  f  B{y,  z)%\—2n[ay-\-vz)\dz  (1.46) 

—4  — c 

has  a  property  defined  by  the  equation: 

AC 

J  flS(y,  z)fdydz=>JJlb(a,  vfdudv.  (I.47) 

i.e.,  the  integral  of  the  square  of  the  function  in  question  over 
its  domain  is  equal  to  the  integral  of  the  square  of  the  modulus  of 
the  amplitude  spectrum  over  the  entire  spatial  frequency  domain. 
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The  correlation  function  Kg (Ay,  Az)  of  a  random  function  B(y,  z) 
is  defined  by  the  expression 

A  C 

A2)=Um^r 

J-oo  J 

c—  ~A  ~C 

X3(cf-^Ay,  z-{.\z)dydz,  (1.48) 

where  *^(y+Ay,  z+Az)  is  the  complex-conjugate  function  of  the  function 
B (y+Ay ,  z+Az) . 

The  inverse  Fourier  transform  of  the  function  b(u,  v)  has  the 

form: 

00 

B{^y,  Az)=|  b(tt,v)%[2:^(u^y-\-v^z]\dudv.  (1.49) 

Inserting  expression  (1.49)  into  (1.48),  we  obtain 

A  c 

/Cb(A«/,  Az)=Um  ^  ^B{y,z)dydzX 

C^m  —.4  — C 

«• 

X  Jj  S{u,  v)»  {  —  2:i[u{y—Sy)-\-v(z—:iZ]\}dadv. 

(1.50) 

Since  the  function  B(y,  z)  is  square- in tegrab le ,  the  integration 
order  in  equation  (1.50)  may  be  changed: 

«• 

c-«»  — «• 

A  C 

Xf  { B{tj,  z)%[—2n(uy-\-vz)\dyUz.  (1.51) 

ZaJc 

Using  equation  (1.46),  we  find 
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A'fllAy,  A2)=Iini  b{u,.v)b{u,  v)X 

WC  J  J 
C-^  .  — «o 

X  X  (2n  (ady  4-  v^z)\  dudv= 

m 

“**'"7^  \i*  !*(“•  '°)T'^{^^{’^b.y-\-v^z)\dudv.  (1.52) 

A-rnm  4i4C  JJ 


C-*ao 


The  limit  lim  — ^  | ft(a, 'DJl*=ftx-fl  (*•  ““)  the  Khinchine-Wiener  spectrum  or 

,4-00  iAC 

C-oo 

the  spectral  density  of  the  dispersion  of  the  random  spatial  distri¬ 
bution  of  radiance. 


Obviously, 


f(f,(^y,  i2)=  J  bx~B(a,  v}*{2n(uiiy-{-v^z)\dudv,  (1.53) 

bji^B  (a,  fA" B ( Ay,! A«)  X 

X  ■'[—2n  (aAy 4-«Az)l  d(^y\d  (A«),  ( 1. 54) 

i.e.,  there  exists  a  relation  describable  by  a  Fourier  transform  be¬ 
tween  the  Khinchine-Wiener  spectrum  of  a  random  function  B(y,  z)  and 
its  correlation  function;  the  correlation  function  of  a  random  func¬ 
tion  B(y,  z)  can  be  obtained  as  the  inverse  Fourier  transform,  of  the 
Khinchine-Wiener  spectrum. 

Formula  (1.54)  can  be  used  to  determine  the  Khinchine-Wiener 
spectrum  of  an  ergodic  background  process,  i.e.,  when  the  correlation 
function  can  be  obtained  from  one  realization  of  the  random  function 
B(y,  z)  for  a  sufficiently  large  area  in  which  it  is  defined. 

For  stationary  and  isotropic  random  radiance  functions,  the 
correlation  functions  and  Khinchine-Wiener  spectra  are  expressed  by 
functions  of  two  variables  that  have  circular  symmetry  in  the  coor¬ 
dinate  systems  (Ay, Az)  and  (u,  v) .  This  makes  it  possible  to  depict 
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these  functions  graphically  as  corresponding  one-dimensional  cuts. 

A  complete  two-dimensional  representation  of  the  functions  can  be 
obtained  by  rotating  the  one-dimensional  cuts  about  the  vertical  eixis. 

Examples  can  be  cited  of  functions  b„  n(u,  v)  that  are  encountered 

X“I3 

in  infrared  equipment  design  practice  [5,  25] .  A  spectrum  of  the 
following  type  is  used  to  represent  an  anisotropic  cloud  background: 

biu,  +  (I.ooal 

where  Ay  and  Az  are  the  so-called  correlation  intervals  for  the  y 
and  z  axes,  i.e.,  the  distances  at  which  the  correlation  function 
amounts  to  37%  of  the  dispersion  of  the  background  luminance. 

In  a  polar  coordinate  system  we  have 


b{w,  y)  sin  29(1 +■»*/=•( (1.55b) 


where 


A4(=A8Cos9;  Az=A(}sin9: 

/■  (  Y)= 2  (sin*  Y  sin*  5 -|- cos*  Y  cos*  9); 

AQ= V' Ajc*-i- A«/*  is  the  correlation  radius  (a  quantity 

analogous  to  the  correlation  interval) . 

The  deviation  of  the  angle  9  from  11/4  characterizes  the  degree 
of  anisotropy  of  the  background  for  a  constant  passband  with  respect 
to  the  modulus  w  of  the  vector  of  spatial  frequencies. 

When  Ay=Az=^  (or  when  9=tt/4  in  a  polar  system),  equation  (1.55b) 
reduces  to  the  form 


*('ir)=2aAQ*Z3*(l-f T»*Ac*)  (1.55c) 


which  corresponds  to  an  isotropic  cloud  background, 
in  Ref.  30  has  an  analogous  form: 

A(«)=2.TAe*lJ*i  1 


The  formula  given 


:i .55d! 
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At  high  frequencies  function  (1.55d)  tends  asymptotically  toward  the 
quadratic  hyperbola 

b~2TTD2w  * 

whose  parameters  are  independent  of  the  correlation  radius. 

It  proves  convenient  to  use  the  Khinchine-Wiener  spectrum  of 
a  heuristic  model  background  to  estimate  the  noise  immunity  of  opto¬ 
electronic  devices  [4.18]: 

b(tt,  v)=BD^  I^y  (1.55e) 

When  Ay#Az,  this  approximately  gives  background  anisotropy,  while 

A  0 

when  Ay=Az=— 2  it  corresponds  to  a  quasi-isotropic  background. 

1.4.  Optical  Signal  Conversion  by  Multidimensional  Filters 

Signal  conversion  in  optoelectronic  systems  is  carried  out  by 
means  of  various  elements  called  filters.  A  distinction  is  made  be¬ 
tween  one-dimensional  and  multidimensional  filters.  One-dimensional 
filters  transform  the  spectral  composition  of  transmitted  radiation. 
Multidimensional  filters  are  used  to  change  the  spatiotemporal  form 
of  signals.  Below  we  will  consider  only  multidimensional  filters. 
These  include  lenses,  condensers,  diaphragms,  modulating  disks,  and 
other  components  used  to  focus,  change  the  direction  of,  cut  off, 
modulate,  separate,  and  add  light  becims. 

A  signal  that  carries  information  on  an  object  or  the  surrounding 
background  is  represented  by  a  two-dimensional  radiance  function; 
the  components  listed  above  therefore  are  called  spatial  filters  as 
well . 


Spatial  filters  differ  essentially  from  electrical  filters.  One 
of  the  differences  is  that  the  independent  variable  of  functions  that 
describe  processes  in  electrical  circuits  is  time,  which  always  is 
greater  than  zero  in  real  processes,  whereas  in  spatial  filters  the 
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independent  variables  of  the  functions  are  the  space  coordinates, 
which  are  unbounded  to  either  side  of  zero.  The  second  difference 
is  that  electrical  filters  are  characterized  by  both  amplitude  pa¬ 
rameters  and  energy  paraimeters  of  signal  transmission;  the  descrip¬ 
tion  of  processes  by  means  of  amplitude  parameters  is  used  more  of¬ 
ten.  Only  the  energy  parameters  of  a  signal  can  be  registered,  by 
contrast,  in  spatial  filters,  which  work  with  noncoherent  radiation. 
Finally,  the  signal  at  the  output  of  an  electrical  filter  is  time- 
delayed  with  respect  to  the  input  signal;  no  such  delay  (in  space) 
occurs  in  spatial  filters. 

Despite  these  differences,  the  methods  employed  in  the  study 
of  electrical  filters  can  be  used  in  the  theoretical  investigation 
of  spatial  filters. 


Fig.  1.10.  Toward  the  deter¬ 
mination  of  the  weighting  and 
spatial  frequency  transfer 
functions  of  a  very  simple 
optical  system:  (1)  object 
plane;  (2)  objective;  (3)  im¬ 
age  plane  (focal  plane) . 

The  most  important  characteristics  of  multidimensional  filters 
are  the  weight  function,  the  unit  step  function,  and  the  spatial  fre¬ 
quency  transfer  function.  We  will  give  the  definitions  of  these  func 
tions  by  using  as  an  example  a  very  simple  optical  system  that  con¬ 
sists  of  an  objective  which  converts  the  radiance  distribution 
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function  for  the  object  plane  nO^  into  the  irradiance  distribution 
function  for  the  image  plane  yOz  (Fig.  1.10). 

Let  no»  Co  and  y,  z  respectively  be  the  coordinates  of  some  point 
of  the  object  and  of  its  ideal  image.  For  numerous  reasons  (aberra¬ 
tion,  diffraction,  etc.),  the  image  of  the  point  extends  to  some  small 
region  whose  irradiance  distribution  can  be  represented  by  the  func¬ 
tion  Ej(y»  z)  if  the  object  point  in  question  is  located  at  the  ori¬ 
gin  0,  and  by  the  function  E^(y-gno»  z-BCo)  if  the  object  point  is 
not  located  at  the  origin;  3=f^g/L. 

The  function  E2(y,  z) ,  which  describes  the  irradiance  distribu¬ 
tion  in  the  image  of  an  isolated  point,  is  the  optical  system's  reac¬ 
tion  to  a  point  source;  it  is  called  the  weighting  function.  To  avoid 
the  difficulty  associated  with  the  presence  of  a  point  source  that 
has  a  finite  radiance,  we  use  the  concept  of  the  spatial  delta  func¬ 
tion  (also  called  the  Dirac  function) : 


^iy—yo^  2-z.)= 


09  if 

0  if 


y=y,and  z  = 

or  ^^0. 


The  spatial  delta  function  is  distinguished  by  the  following 
properties; 


^(y-y„,  z  —  z„)dydz=l;  (1.56) 

00 

jj  ■2-r,)/(y.  z)dydz=f{!/^,  zj,  (1.57) 

where  f(y,  z)  is  a  continuous,  bounded  function. 

Depending  on  the  type  of  weighting  function,  multidimensional 
filters  are  divided  into  stationary  and  nonstationary  types.  Whereas 
the  weighting  function  of  a  stationary  filter  is  independent  of  its 
movement  along  the  coordinate  axes,  a  nonstationary  filter  has  such 
a  dependence. 
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The  function  E^(y»  z)  is  an  important  characteristic  of  an 
optical  system.  If  we  know  how  an  optical  system  converts  the  radi¬ 
ance  of  a  point  source,  then  we  can  determine  how  it  alters  the  ob¬ 
ject's  radiance  distribution  function,  as  the  object  can  be  represented 
as  an  infinite  number  of  point  sources  that  have  a  radiance  equal  to 
the  radiance  of  the  corresponding  point  of  the  object. 

By  writing  Sno=y'  and  SCo=z'  and  representing  the  object's  ra¬ 
diance  distribution  function  as  B(rio,  Co)=B*(y'/  z'),  we  can  obtain 
the  irradiance  distribution  function  of  the  object  image  E{y,  z)  as 
the  sum  of  the  irradiances  of  the  images  obtained  from  different 
points  of  the  object: 


E(y,  z)  =  jj  z)Ei['i-y\  z-z')di/ dz'.  (1.,58) 

Relation  (1.58)  shows  that  E(y,  z)  is  the  convolution  of  the 
functions  B*  and  E^,  which  may  be  conditionally  written  as 

£(//.  z)=/9»  ^  0.59) 

The  limits  of  integration  of  (1.58)  are  determined  by  the  bor¬ 
ders  of  the  object,  i.e.,  by  the  range  of  real  values  of  B*(y,  z) 
or  by  the  system's  field  of  view. 

We  will  use  a  Fourier  transform  to  solve  equation  (1.59).  In¬ 
troducing  the  variables  u  and  v,  which  are  spatial  frequencies,  we 
may  write  the  Fourier  transform  %(\x,  v)  of  the  function  E(y,  z)  as 

^(l*.  '')=]*!  E(y,  z)*[—2n{v-y i-vz]\dy dz\  (1,60) 

=  v) X V.  (1.61) 

Replacing  E(y,  z)  by  its  value  according  to  (1.58),  we  find 

z-z')x 

X*[  —  2:i{^y-\-\z)\dy'dz'dydz,  \  1,62) 

or,  using  the  new  variables  Y*y-y'  and  Z=z-z', 
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X£.(K,  Z)*[-2:i(iiy i.\Z)]dy'dz'dy'dZ.  (1.63) 

This  expression  can  be  easily  reduced  to  the  form 

«((»,  v)=ft*(ii,  v)^fs(ti,  v).  (1,64) 

i.e.,  the  Fourier  transform  of  the  irradiance  distribution  function 
of  the  object  image  is  equal  to  the  product  of  the  Fourier  transform 
of  the  object's  radiance ' distribution  function  by  the  Fourier  trans¬ 
form  of  the  weighting  function. 

The  direct  Fourier  transform  of  the  irradiance  distribution  func¬ 
tion  of  the  image  of  an  isolated  point  (of  the  weighting  function) 
is  called  the  spatial  frequency  transfer  function  [n^n®]  of  the  op¬ 
tical  system.  If  the  weighting  function  E^(y,  z)  depends  on  two  co¬ 
ordinates,  then  the  spatial  frequency  transfer  function  is 

m 

Ei[y,  (1.65) 

The  spatial  frequency  transfer  function,  which  is  the  analog 
of  the  transfer  fiinction  used  in  automatic  control  theory,  can  be 
represented  as 

r(/(i,  /v)= />({»,  v)  +  /Q((i,  V),  (1.66) 

or 

«?'{/(»,  iv)=»'(l»,  (1.67) 

where  P(u,  v) ,  Q(u,  v) ,  W(u,  v) ,  and  flu,  v)  respectively  are  the 
real,  imaginary,  eunplitude,  and  phase  frequency  responses.  The  fol¬ 
lowing  relationships  exist  among  these  characteristics: 

(|»,  v)=V  (^{|*.  v)l»+l(>(|i,  v)l*;  (1.68) 

V  (1»,  v)=arctg  .  (1.69) 
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If  the  weighting  function  is  axisymmetric ,  then  the  amplitude 
frequency  response  can  be  found  by  using  a  Hankel  transform; 


where  Q=Yy*+z*-,  V  !»•*+>*• 

In  this  case  W(w)  has  rotational  symmetry. 


When  the  optical  signal  depends  on  three  variables  (the  space 
coordinates  y  and  2  and  the  time  t) ,  the  spatial  frequency  transfer 
function  of  a  multidimensional  filter  is 

«•  iv.  (y.  2.  (1. 70) 

In  order  to  give  a  physical  interpretation  of  the  spatial  fre¬ 
quency  transfer  function,  let  us  examine  the  expressions  for  the  in¬ 
put  and  output  signals  of  a  multidimensional  filter.  As  noted  above, 
the  light  signal  (y,  z)  at  the  filter's  input  can  be  represented 
as  an  infinite  sequence  of  spatial  delta  functions.  Based  on  property 
(1.57)  of  the  delta  function,  we  have 

•Pofy-  •2)=JJ  z  —  n)dmdn.  (1.7n 

To  find  the  signal  at  the  filter's  output,  let  us  recall  that 
the  weighting  fvuiction  is  its  reaction  to  a  signal  in  the  form 
of  the  spatial  delta  function.  Therefore,  each  of  the  infinitesimal 
light  signals 

V'o  (m,  n)  5  (y-m,  z-n)dm  dn 

causes  an  elementary  signal  at  the  filter's  output; 

(m,  n) Eg  (y-m,  2-n)dm  dn 

The  resultant  output  signal  is 


<Pi(y.  z)= JJ  — m,  z—n)<f„(m,  n)dmdn.  (1.7l.h 


The  integral  in  the  right  side  of  expression  (1.72)  is  a  two-dimen¬ 
sional  integral  of  convolution,  which  is  written  as  follows  in  vector 
form: 


•Pile)  I  Et(Q—7)<f^Cr)dr. 


(1.73) 


I  Taking  the  Fourier  transforms  from  both  sides  of  equation  (1.73), 

we  f ind 


f,  (0  — J<^0.  (1.74) 

Introducing  the  new  variable  A=p-r,  we  have 

(tm)  =  f  e-'"^  {r)dl  f  c—  (r)  dr^  (I.  75) 

Hence  it  follows  that 

F^ilZ)  =  W{tZ)F,{U), 

or 


i.e.,  the  spatial  frequency  transfer  function  of  a  multidimensional 
filter  is  the  ratio  of  the  Fourier  transform  of  the  output  signal 
to  the  Fourier  transform  of  the  input  signal . 


Before  moving  on  to  a  description  of  some  properties  of  spatial 
frequency  transfer  functions,  let  us  make  two  remarks.  The  first 
pertains  to  the  optical  system's  reaction  to  the  spatial  delta  func¬ 
tion.  Because  of  diffraction  effects,  even  an  ideal  optical  system 
in  which  there  are  no  aberrations  produces  the  image  of  an  isolated 
radiating  point  in  the  form  of  a  spot  of  finite  dimensions  whose  ir- 
radiance  decreases  rapidly  from  its  center  to  its  periphery. 


for  the  image  of  an 
isolated  point. 


The  distribution  of  the  irradiance  in  the  image  of  an  isolated 
point  that  is  produced  by  an  ideal  instrument  is  described  by  the 
function  [19,  22] : 


where 


£i{y.  z)= 


(1.77) 


a' — half  the  angle  at  the  vertex  of  the  cone  of  rays  that 
forms  the  image  (Fig.  1.10); 

\ — the  wavelength; 

— a  first-order  Bessel  function  of  the  first  kind; 

E  — the  maximum  value  of  the  function  E.(y,  z)  that  corresponds 

M  0 

to  the  values  of  the  coordinates  y=0  and  z=0. 

The  function  (y,  z)  can  be  represented  as  a  surface  (Fig.  1.11) 
which  bounds  a  volume  called  the  "diffraction  body."  If  the  fxinction 
E^(y,  z)  is  symmetric,  then  it  is  represented  as  a  one-dimensional  cut. 
Table  1.5  gives  the  values  of  E^/E^^  for  various  n.  The  function  E^ 
vanishes  for  n  equal  to  3.83,  7.02,  10.17,  13.32,  16.47,  etc.,  and 
has  a  maximum  at  the  following  values  of  n;  5.13,  8.42,  11.62,  14.80, 
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3 


1 


Table  1.5.  Values  of  Ej./E  for  Various  n 

sj'  M 


I 


n 

n 

Eu 

n 

E, 

0.0 

1.000 

2,8 

0,0856 

5,6 

0.0153 

0.4 

0.961 

3.2 

0,0267 

6.0 

0.0085 

0,8 

0.850 

3.6 

0,0028 

6,4 

0.0032 

1.2 

0.690 

4.0 

0,0011 

6.8 

0,0004 

1.6 

0,507 

4.4 

0.0085 

7.2 

0.0003 

2.0 

0.333 

4.8 

0.0155 

7,6 

0,0018 

2.4 

0.188 

5,2 

0.0174 

8.0 

0.0034 

I 


17.96,  etc.  Ratios  E^/E^  of  100%  are  respectively  equal  to  1.75, 
0.416,  0.160,  0.078,  and  0-044  [22]. 


I 


The  function  E^ (y,  z)  has  no  analytical  expression  for  real  op¬ 
tical  systems.  As  a  result,  experimental  methods  of  determining  this 
function  and  approximate  theoretical  methods  that  take  into  account 
the  effect  of  spherical  and  chromatic  aberrations  only  have  been  under 
development  recently. 


The  second  remark  relates  to  the  optical  system' s  response  to 
a  signal  in  the  form  of  a  unit  step  function.  If  the  object  plane 
contains  a  sharp  boundary  between  an  illuminated  quarter  of  the  plane 
and  the  other  quarters,  which  are  dark,  i.e.,  if  B*(y,  z)  is  equal  to 
unity  when  y>0  and  z>0  and  to  zero  otherwise,  then  the  irradiance 
distribution  in  the  image  of  this  picture  can  be  found  from 
formula  (1.58): 


J  J  Ei[y-y\  z-z')dy'dz‘.  (1.78) 


Integrating  over  y-y'=Y  and  noting  that  dy’=-dY,  we  obtain 

«• 

^*(y.  *)=“  [  o{z—z')dz\  (1. 79) 

0 
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The  function  a(z)  is  equal  to  the  cross-sectional  area  of  the  dif¬ 
fraction  body  where  it  is  intersected  by  a  plane  having  an  ordinate 
2  (Fig.  1.11) . 


The  function  ^hiy,  2).  which  represents  the  reaction  of  a  multi¬ 
dimensional  filter  to  a  unit  light  signal  propagated  in  one  quadrant, 
is  called  the  unit  step  function  (the  Heaviside  unit  function) .  It 
is  analogous  to  the  unit  step  function  of  the  linear  section  of  an 
automatic  control  system — the  section's  reaction  to  a  unit  step. 

Returning  to  our  examination  of  spatial  frequency  transfer  func¬ 
tions,  let  us  note  an  important  aspect  that  follows  from  its  I  sic] 
mathematical  definition.  If  we  Icnow  the  complex  spatial  frequency 
spectrum  b_„(iu»  iv)  of  a  signal  at  the  input  of  the  transfer  section 
and  its  spatial  frequency  transfer  function  W(iia,  iv)  ,  then  the  sig¬ 
nal  b„„„(iU/  iv)  at  the  output  of  the  section  is  defined  by  the 
relation 

ty).  (1,80) 

If  several  sections  are  connected  in  series,  the  signal  at  their 
output  will  be 


^v)  n (i.sn 

n^l 

where  m  is  the  number  of  series-connected  sections. 

Expression  (1.81)  follows  directly  from  expression  (1.80)  if 
the  output  signal  of  each  preceding  section  is  considered  the  input 
signal  of  the  next  section. 

In  most  cases  we  are  concerned  only  with  pealc  signals  that  can 
be  represented  by  their  amplitude  spatial  frequency  spectra.  In  these 
cases  expression  (1.81)  assumes  the  form 

I  *„„(/(»,  /v)|=16^(ll».  fv)|  fj  VI.  (1.82) 


Fig.  1.12.  Evaluation 
of  the  properties  of  a 
transfer  section  based 
on  its  spatial  fre¬ 
quency  transfer  function. 

The  spatial  frequency  transfer  function  completely  describes 
the  transfer  properties  of  individual  sections  and  of  an  entire  op¬ 
tical  information  transmission  system.  It  also  makes  it  possible 
to  synthesize  a  system.  Efforts  are  made  to  ensure  that  the  curve 
of  the  spectral  frequency  transfer  function  comes  as  close  to  unity 
as  possible  in  the  range  of  spatial  frequencies  used.  To  compare 
two  sections,  we  usually  select  characteristic  points  on  the  curve 
of  the  spatial  frequency  transfer  function  that  correspond  to  the 
levels  0.8,  0.5,  e  ^=0.37,  and  so  forth  (Fig.  1.12).  The  selection 
of  a  given  level  as  a  criterion  for  evaluating  a  section  depends  on 
its  function. 

In  synthesizing  two-dimensional  systems  from  the  spatial  fre¬ 
quency  transfer  functions  of  individual  sections,  one  should  keep 
in  mind  that  the  spatial  frequency  transfer  function  of  each  section 
is  defined  in  a  certain  plane;  in  recomputing  the  spatial  frequency 
transfer  function,  allowance  should  be  made  for  scale  transformations 
of  the  sections. 

The  mathematical  tools  considered  above  for  determining  the 
spatial  frequency  transfer  function  of  two-dimensional  systems  are 
valid  only  if  these  systems  are  linear  and  if  the  light  beam  that 
carries  information  regarding  an  object  is  produced  by  a  source  of 
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noncoherent  radiation.  Let  us  determine  the  spatial  frequency 
transfer  functions  for  some  elementary  filters. 

In  the  absence  of  diffraction  effects,  an  ideal  nonaberrational 
lens  is  an  amplifier  section  whose  weighting  function  is  expressed 
in  terms  of  the  spatial  delta  function  as  follows; 

Et.{y,  z)=k^{y,  z). 

Obviously, 


(/J*.  /v}=k  JJ* «  (y.  z)  'e-^‘  11^1'+  «)  dy  dz.  (1. 83) 

It  is  not  hard  to  see  that 


h)=k.  (1.84) 

An  optical  shifter  changes  the  direction  of  propagation  of  light 
rays  without  changing  the  waveform.  It  has  the  weighting  function 

2)=l{y—y^  Z  —  Zo). 

Based  on  property  (1.57)  of  the  spatial  delta  function,  we  have 

W  (/(»,  ( 1. 85) 

A  nonaberrational  lens  stretches  the  image  of  a  point  source  of 
optical  radiation  into  a  line.  Its  weighting  function  is  represented 
by  a  one-dimensional  delta  function; 


£‘i(y.  Z)=9(y-y„), 

where  yg  is  the  displacement  along  the  y  axis. 

The  Fourier  transform  of  the  weighting  function  is 

/v)=»»(v)e-**''‘*».  (1.86) 
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A  Diaphragm  with  a  Rectangular  Aperture.  The  weighting  func¬ 
tion  in  the  form  of  a  rectangular  parallelipiped  has  the  following 
analytical  expressions: 


E,.{y,  z]  = 


1  when  — -^and— 
CWhen</>^.  ^<—7  or  ^ 


£1 

•> 


In  this  case 

-iL  £l 

2  J 

(/(», /v)=  1'  j  1.871 

"yt  it 

~T-r 

Separating  the  integrals  and  integrating,  we  find  [7] : 

f/ji,  i\)=  —^  sin  n(ti/,sin  .nvcj,  (I.881 

A  Diaphragm  with  a  Round  Aperture.  Setting  E^(p)=l  and  using 
the  Hankel  transform,  we  find 


W  (ui)=2n  J  e/o(C“>)‘^C. 

where  R  is  the  radius  of  the  circle. 


,1.891 


Taking  into  consideration  that 


d  [(<?•)"/, to*"' 


we  obtain 


A  Lens  with  Aberrations.  The  light-scattering  function  for  the 
aberrational  point  of  a  lens  can  be  approximated  as  follows  [6]: 
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I.‘)l 


where  Pi=0.4769 — a  constant 

— the  radius  within  which  50%  of  the  energy  in  the 
light  signal  is  concentrated. 


Obviously, 


W  (/[i,  l\)= 


e? 

77? 


—Ol  — 


~0l  — X-  — 

e  'i 


dz.(1.92 


Without  giving  a  detailed  solution,  we  will  present  the  final 
result  [6 ,  7] ; 

'■'i 

\ru>,  .  (1.93) 

A  Filter  with  Complex  Weighting  Functions.  If  the  weighting 
function  of  a  filter  has  a  complex  expression  or  is  determined  exper¬ 
imentally,  then  it  may  be  divided  approximately  into  regions  whose 
frequency  responses  are  readily  determined  analytically.  For  example, 
when 


m  m 


Ei{y,  z)=»2  2 


we  have 


/v)-.  22  /V).  (1.94) 

kml 


Dividing  the  complex  weighting  function  into  elementary  rectangles 
(Fig.  1.13a),  we  find 


/v)=^S-sin  n|ib^u-sin  nvAz-e* 
nV* 


<(.**,+«») 
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Fig.  1.13,  Calculation  of  the  spatial 
frequency  transfer  function  of  a  fil¬ 
ter  with  a  complex  weighting  function 


where  — the  amplitudes  of  elementary  components  of  the 

weighting  function; 

Ay,  Az — the  dimensions  of  the  elementary  rectangles; 
y^,  — the  coordinates  of  the  centers  of  the  rectangles. 

Obviously, 


'^sin  n|iAy-sin  nv^zx 


(1.95) 


Since 


llm  — —  sin  npAtf-sin  nvA^  =  Atf•A^, 


we  have 


«  m 

W  (/,*.  iv)=,2  2  ( 1. 96) 

r-l  *•! 

If  the  weighting  function  has  rotational  symmetry,  then  it 
advisable  to  divide  the  function  into  cylindrical  rings  (Fig.  1 
for  each  of  which 

^W,  {»)-  >l,2n  [  RUt  _  R>  ii^l  ,  1 1 . 97( 


is 

13b)  , 
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where  and 


are  the  radii  of  the  ring. 


In  this  case 


(1.98) 


Passing  to  the  limit  for  small  thickness  of  the  cylindrical  ring 
we  obtain 


W  («)=: 2n  y  A,/^  (A*,-)  H .  (i-h 

( 

A  Filter  with  Time-Varying  Weighting  Functions.  First  let  us 
consider  the  case  where  the  weighting  function  undergoes  translation 
in  space  and  has  the  form 


Ic+^Wl,  (1. 100) 

where  r(t)  is  the  displacement  vector. 

Based  on  the  property  of  the  Fourier  transform  of  a  displaced  signal 
that  if 


/(q)=A(c^'). 


then 


we  find 


W'(»i./)=\r  (to) (1.101) 

where  W(iu))  is  the  spatial  frequency  transfer  function  of  a  filter 
whose  weighting  function  is  found  at  rest. 

If  the  weighting  function  is  rotated  by  an  angle  Y*Y(t),  the  spa 
tial  frequency  transfer  function  of  the  filter  has  the  form: 


t)=W  \iM, 


(1. 102) 


where 


Mud 


cosy(0  —  sinY(0 
sinY(/)  cosyW 


is  the  rotation  matrix. 

If  we  know  the  spectrum  of  the  input  signal  and  the  spatial  fre¬ 
quency  transfer  function  of  a  multidimensional  filter,  then  we  can 
calculate  the  spectriam  and  spatial  form  of  the  output  signal .  As 
an  example,  let  us  calculate  the  spectrum  and  shape  of  the  output 
signal  of  a  multidimensional  filter  with  a  bell-shaped  weighting  func¬ 
tion  to  whose  input  a  unit  signal  propagated  in  a  quadrant  of  the 
half-plane  is  supplied  (Fig.  1.14a). 


Fig.  1.14.  Calculation  of  the 
spectrum  and  shape  of  the  output 
signal  of  a  multidimensional 
filter. 


Analytically,  such  a  signal  is  expressed  in  the  following  form: 

f  (V  z)=l^ 

0  '  |0  for  other  values  of  y  and  z 

The  input-signal  spectrum  is 


^ o(*>.  <v)=  I  r  dydz- 


4ji!^v  * 


(1. 103) 
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The  amplitude  characteristic  of  the  spectrum  is  shown  in  Fig.  1.14b. 
By  stipulation,  the  spatial  frequency  transfer  function  of  a  filter 
is  defined  by  the  relation 


W  {fy,  iv)=e 


e? 


Therefore,  the  output-signal  spectrum  is 


or 


- 5- (!*•+.*) 


(1. 104) 


Fx (/|*.  (/>,  /v)  W  (/^  /v).  ( 1.  105) 


(1. 106) 


The  shape  of  the  output  signal  is  determined  by  the  inverse 
Fourier  transform: 


m 

Jjir 


(»»’+»•> 


Calculating  the  integral,  we -obtain 


(1. 107) 


where 


*‘W— 


e  rfjc  is  the  reduced  error  function. 


As  a  second  example,  let  us  consider  the  passage  of  a  rectangular 
signal  through  an  optical  shifter  whose  weighting  function  rotates 
at  some  constant  angular  speed  [6,  7], 


The  spatial  frequency  transfer  function  of  an  optical  shifter 
with  a  stationary  weighting  function  according  to  (1.85)  has  the 
form 
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k 


i 


Upon  rotation  of  the  weighting  function. 


WW{t),  VW] {1. 109) 

The  component  frequencies  are  determined  by  the  rotation  matrix  and 


j»'=jtcosY<— vsin 
v'=i»sin  ^-|-''Cos^, 


(1.110) 


so  that 


W  (/(»,  jv,  t)—er-'^*  (►».+«.)+  la j  j 


The  shape  of  the  output  signal  is  determined  by  the  inverse  Fou¬ 
rier  transform; 


«i 

X-;^  sin  nf^y^  sin  7ivz^d\^dv.  (i.  112) 


After  transformations  which  we  omit,  we  find  [7] ; 


/j(/)*^(y,cos  ^+z,sin  iJ<)t»(*,cos  ^-y.sin  yT), 

(l«  113) 


where 


0  if  y,cos-v<+z,sin  •y#<  — -fiand 
y,cos^+«,sin  Y/>-^; 

1  if --y<y*cos>i+*,sin 
0  if  z,  cos  ■{<  —  y,  sin  yi  <  —  and 
Zacosy  — y,  sin  ^  ; 

1  if  — ^<z,cosi^-y,sin^<-^. 
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At  the  filter's  output  the  signal  has  the  form  of  rectangular 
pulses. 


CHAPTER  2 


TIME-AND-FREQUENCY  CHARACTERISTICS  OF  MODULATORS 

2.1.  Analytical  Determination  of  the  Time-and-Frequency 

Characteristics  of  Modulators 

If  a  radiant  flux  FQ(t)  is  incident  on  a  modulating  disk,  then 
the  value  of  the  modulated  radiant  flux  is 

F(t)=FQ(t)T(t)  (2.1) 

where  T(t)  is  the  tranmission  function  of  the  modulating  disk. 

The  time-and- frequency  characteristics  of  signals  picked  up  from 
the  radiant-energy  detector  depend  on  the  configuration  of  the  raster 
of  the  modulating  disk.  By  expanding  the  transmission  function  of 
the  modulating  disk  into  a  Fourier  series,  we  can  estimate  the  char¬ 
acteristics  of  signals  going  from  the  output  of  a  radiant-energy  de¬ 
tector  to  the  selective  elements  of  an  optoelectronic  device. 


Fig.  2.1.  Modulating  disks: 

(a)  with  a  transparent  sector; 

(b)  with  alternating  transpar¬ 
ent  and  opaque  sectors. 


One  of  the  simplest  modulating  disks  is  a  slit  modulator,  which 
is  implemented  as  a  circle  with  a  sector  that  differs  in  its  trans¬ 
mission  coefficient  from  the  rest  of  the  surface  (Fig.  2.1a). 


The  transmission  function  of  a  slit  modulator  in  relation  to 
the  angle  of  rotation  can  be  represented  as 


t(a)  = 


0<a<o„  aj<a<2.n. 


t(0= 


If  ^  iL 

2  ^  T  IT  ’ 


t„.  if  0<^<^,  T<f<^> 


where  il  is  the  angular  rotational  velocity  of  the  modulating  disk. 

Using  the  relations  presented  in  Chapter  1,  we  can  represent 
the  coefficients  of  the  expansion  of  a  slit  modulator's  transmission 
function  into  a  Fourier  series  by  the  following  expressions: 

a>=  (sin  *a,-sin  *a,); 

**  =  ( cos  *a,  —  cos  ka.);  ( 2. 4) 

^02  —  01  , 

^  2  * 

?,=arctg  ^  °2  +  g|  ^ 

sin  tao  —  sin/hX|  2 

According  to  formulas  (2.4),  the  transmission  function  of  a  slit 
modulator  will  assume  the  form 


q? — Oi 


- - cos  *  (ut  -f  I . 


(!'.  ■■5) 


For  a  modulating  disk  patterned  after  Fig.  l.lf,  a2-ai=7r.  If 

we  assume  x  =>1,  t  .  =0,  and  ai=0,  then  the  transmission  function 
in^x  mxn 

of  the  disk  is 


2.6 


For  even  k,  sin  k|=0,  and  expression  (2.6)  simplifies,  assuming 
the  form 


9.0 


(2.7) 


A  "uniform  grating" -type  modulating  disk  is  shown  in  Fig.  l.le. 
Depending  on  the  angle  a  of  rotation  of  the  disk,  its  transmission 
function  is  determined  by  the  relation 


r(a)= 


2r 


I  if  ~  ('«  — 1)  <<*'< -j- (2/n— 1), 
0  if  —  (2m-l)<a<5Lff,, 


(2.8) 


where  z — the  number  of  transparent  (or  opaque)  sectors; 
m — the  sector  number  (m=l ,  2,  3,  ...,  z) . 


Since  the  modulating  disk  rotates  with  a  content  angular  velocity 
Q,  the  transmission  function  can  be  expressed  by  the  formula 


T(/)  = 


lif  ^(/n-l)<-J<  J-(2m-l). 
Oif  m. 


(2.9) 


za 


za 


2Tr 


The  period  of  the  function  (t)  is  the  quantity  and 

is  the  modulation  frequency  of  the  radiant  flux.  The  expression  for 
the  transmission  function  of  a  modulating  disk  of  the  "uniform  grat¬ 
ing"  type  is  analogous  to  expression  (2.7),  but  with  the  difference 
that  the  angular  frequency  to  of  modulation  of  the  radiant  flux  will 
figure  in  this  equation  in  place  of  the  angular  rotational  velocity 
of  the  disk. 


The  raster  of  the  modulating  disk  shown  in  Fig.  2.1b  has  a  more 
complex  pattern.  The  disk  consists  of  two  halves  each  of  which  trans¬ 
mits  50%  of  the  incident  radiation,  while  the  other  half  consists  of 
alternating  transparent  and  opaque  sectors.  The  transmission  function 
of  such  a  disk  is  represented  by  the  expression 


(2.  10) 


1  if  —  (m- 

T(a)=  0  if  l)<a<^  fft, 

if  n<a;<2n. 

where  z — twice  the  number  of  sector  pairs; 

m — the  number  of  a  sector  pair  (m-1,  2,  3,  z/2)  . 


(a) 

(b) 

(c) 

(d) 

Fig.  2.2.  Waveform  at  the  output  of 
a  radiant-energy  detector  (a) ,  and 
the  type  of  auxiliary  functions  that 
allow  representation  of  the  trans¬ 
mission  coefficient  of  a  modulating 
disk  as  a  linear  combination  (b,  c, 
d)  . 

The  waveform  at  the  output  of  a  radiant-energy  detector  in  front 
of  which  the  modulating  disk  in  question  is  positioned  is  shown  in 
Fig.  2.2a.  For  Fourier  expansion  of  the  transmission  function  T(t) 
of  a  modulating  disk  with  a  complex  raster  pattern,  it  is  convenient 
to  represent  T(t)  as  a  linear  combination  of  simple  functions: 

T(t)=Ti(t)T2(t)+T3(t)  (2.11) 

The  form  of  the  functions  Ti(t),  T2 (t) ,  and  Tslt)  is  shown  in 
Figs.  2.2b,  c,  and  d.  The  elementary  functions  in  expression  (2.11) 
have  the  form 

2^^ sin  (2^+  (2. 12) 
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m 

tiW=Y+-7^^^^sin(2?+l)2/;  (2.13) 

T,(0= Y  •  (2. 14) 

The  transmission  function  of  the  modulating  disk  is 

»(<)— 5-+-^  (2«+  'I  •<]  X 

P.I5) 

Some  optoelectronic  devices  use  a  "uniform  grating"  modulating  disk 

that  has  variable  optical  density  over  its  diameter  (Fig.  1.1k).  This 

is  accomplished  by  using  a  material  whose  transparency  obeys  a  wedge 

law.  The  transmission  function  of  such  a  modulating  disk  (provided 

that  and  t„.„=0)  is  expressed  by  the  series 

max  min 


"'>“[t+T  §  TJTT''"  <^+  ”  ■']  X 

As  we  can  see  from  formula  (2.16),  the  factor  in  brackets  is 
the  transmission  function  of  a  "uniform  grating"  modulating  disk, 
while  the  second  factor  determines  the  change  in  the  optical  density 
of  the  disk  material  according  to  a  wedge  law. 

The  problem  of  precision  of  modulator  fabrication  arises  in  solv¬ 
ing  practical  problems  of  the  design  of  optical  apparatus  for  measur¬ 
ing  angular  coordinates.  All  changes  in  the  configuration  of  the 
raster  pattern  can  be  reduced  to  a  change  in  the  dimensions  of  the 
transmitting  (or  nontransmitting)  sectors  or  to  a  change  in  the  trans¬ 
mission  function. 

The  transmission  function  of  a  modulating  disk  with  deviations 
from  the  specified  raster  pattern  can  be  represented  as  follows: 
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Fig.  2.3.  Waveform  at  the  output  of 
a  radiant- energy  detector  when  a 
modulating  disk  having  a  nonideal 
raster  pattern  is  used  (a) ,  and  the 
form  of  auxiliary  functions  that 
produce  AT(t). 

T*(0='fU)+At,/),  (2.17) 

where  T(t) — the  transmission  function  of  an  ideal  modulating  disk; 

At (t) — the  change  in  the  transmission  function  due  to  deviations 
of  the  actual  raster  pattern  from  the  ideal. 

Fig.  2.3a  shows  the  waveform  at  the  output  of  a  radiant-energy 
detector  for  a  modulating  disk  implemented  as  in  Fig.  2.1b  that  has 
an  additional  transmitting  sector.  Figs.  2.3b  and  2.3c  show  the  aux¬ 
iliary  functions  T*(t)  and  x*(t)  which  are  used  to  facilitate  the 
representation  of  At (t)  as  a  Fourier  expansion. 

The  width  of  the  additional  transmitting  sector  is  assumed  equal 
to  0 ,  the  distance  between  the  primary  transmitting  sector  and  the 
secondary  sector  is  At,  and  the  secondary  transmitting  sector  is  lo¬ 
cated  in  a  pair  of  (transmitting  and  nontransmitting)  sectors  num¬ 
bered  m. 

Since 

AT(t)=T*(t) •T*{t)  (2.18) 

by  expanding  each  of  the  functions  in  the  right  side  of  expression 
(2.18)  into  a  series,  it  is  possible  to  determine  the  change  in  the 
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transmission  function  of  the  modulating  disk. 


For  simplicity  we  will  assume  x*  =1  and  x* .  =0 , 

max  min 

the  notations  introduced  above,  we  have 


1  if  -5-4.^ 


0  if  0</<iL-j.^, 

+-</<—; 
<i>  {i> 

1  if 


t;w= 


Allowing  for 


(2. 19) 


0  if  ^(m_l)</<g.(n,_i)+  (2.20) 


The  change  in  the  transmission  function  of  a  modulating  disk 
due  to  deviations  of  the  actual  raster  pattern  from  the  ideal  is  rep¬ 
resented  by  the  Fourier  series 


j  m  - 

V  jCx 


In  analyzing  expression  (2.21)  we  should  note  that  the  additional 
signal  generated  at  the  radiant-energy  detector  carries  information 
at  the  oj-carrier  and  ^-envelope  frequencies.  This  signal  will  be 
picked  up  by  selective  elements  of  information  processors  and  will 
affect  the  output  signals  of  the  control  elements  of  the  equipment. 


A  modulating  disk  with  nonuniform  width  of  the  transmitting  and 
nontransmitting  sectors  (Fig.  1.11)  is  used  in  optoelectronic  instru¬ 
ment  making.  To  represent  the  transmission  function  of  such  a  disk 
as  a  Fourier  series,  we  must  determine  the  expansion  coefficients  for 
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u)=(i)(t).  Even  in  the  simplest  case  where  u)=(jOo+Ct 


,  finding  the 


coefficients  of  the  series  reduces  to  calculating  an  integral  of  the 


T 

type  /  sin  t^  dt.  This  integral  in  turn  is  represented  by  a  series 
0 

of  the  type 


I 


sini*c//= 


(2n  +  I)l(4n  +  3) 


2  /  n  r?  ■  m  ris  . 

42  ‘1320  83160  ”^ 


.  .  (2.22) 


The  need  for  a  large  volume  of  computations  forces  us  to  seek 
methods  that  make  possible  a  reduction  of  the  amount  of  time  needed. 
Electronic  simulation  or  modeling  equipment,  in  particular,  may  be 
used  for  this  purpose. 

2.2.  Simulation  of  the  Time-and-Frequency  Characteristics  of 
Modulators  Using  Analog  Computers 

When  analog  computers  are  used  to  simulate  the  time-and-f requency 
characteristics  of  modulators,  sine-wave  voltages  of  different  fre¬ 
quencies  must  be  added.  A  sinusoidal  voltage  can  be  obtained  in  an 
analog  computer  [ABM]  by  solving  the  equation 

y"+ay=0  (2.23) 

under  certain  initial  conditions. 

The  solution  of  equation  (2.23)  is 

y=A  sin  (u)t+\())  (2.24) 


where 


Fig.  2.4.  Block  diagram  of  the 
solution  of  the  equation  y"+ay=0. 


Fig.  2.5.  Symbols  for  decision 
elements  of  an  analog  computer. 
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To  solve  equation  (2.23)  on  an  analog  computer,  we  must  have 
two  integrators,  one  inverter,  and  three  potentiometers.  A  block 
diagram  of  the  solution  is  shown  in  Fig.  2.4.  The  symbols  for  the 
decision  elements  accord  with  Fig.  2.5.  The  transfer  ratios  of  poten¬ 
tiometers  2  and  3  depend  on  the  initial  conditions,  and  the  value 
of  the  transfer  ratio  of  potentiometer  1  is  determined  by  the  value 
of  OJ*  . 


There  are  two  reasons  why  it  is  extremely  difficult  to  use  this 
method  to  add  the  components  of  a  series.  The  first  is  that  there 
must  be  a  large  number  of  decision  units  in  the  analog  computer.  The 
number  N  of  integrators  and  inverters  is 

N=3n+1 

where  n  is  the  number  of  terms  in  the  expansion. 


The  second  factor  that  determines  the  practical  inapplicability 
of  adding  sinusoidal  components  is  the  fact  that  the  addition  of  the 
sinusoidal  components  of  a  series  is  extremely  difficult  work  that 
requires  much  experience  and  knowhow,  even  with  large  analog  computers 
on  which  the  operator  can  work  without  having  to  look  for  reserve 
decision  elements.  The  sinusoidal  components  must  have  multiple  fre¬ 
quencies;  it  is  practically  impossible  to  determine  exactly  the  value 
of  the  transfer  ratios  of  units  and  potentiometers.  And  if  the  op¬ 
eration  of  a  modulator  must  be  analyzed  over  a  single  period,  the 
addition  of  sinusoidal  components  makes  it  difficult  to  correct  the 
transfer  ratios,  which  determine  the  accuracy  with  which  the  required 
sinusoidal  voltage  frequency  is  set. 

In  investigating  the  operation  of  a  modulator  over  a  time  t>>TQ, 
the  method  of  adding  sinusoidal  components  cannot  be  used.  Even  slight 
inaccuracy  (a  fraction  of  a  percent)  in  the  setting  of  the  coefficients 
that  determine  the  frequencies  of  the  sinusoids  being  added  will  lead 
to  phase  mismatch  of  the  terms  in  a  few  periods.  Therefore,  a  number 
of  artificial  procedures  are  used  to  simulate  time-and- frequency 
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characteristics.  The  simplest  signals  to  modulate  are  those  which 
have  passed  through  a  modulating  disk.  This  is  done  by  using  a  set 
of  nonlinear  units  and  inf ralow- frequency  equipment.  For  example, 
a  low-frequency  periodic  oscillator  (such  as  the  NGPK-2  or  the  NGPK-3) 
is  used  to  produce  signals  that  simulate  the  operation  of  the  modulat¬ 
ing  disks  shown  in  Fig.  l.lf  and  l.le.  [When  such  a  generator]  is 
combined  with  a  multiplier,  it  is  possible  to  construct  the  model 
modulator  shown  in  Fig.  2.1b.  A  circuit  that  solves  relation  (2.11) 
is  constructed  for  this  purpose - 

The  use  of  inf ralow- frequency  equipment  imposes  major  constraints 
on  the  frequencies  of  the  simulated  signals.  The  actual  frequency 
values  are  limited  to  a  few  tens  of  hertz.  The  introduction  of  a 
nonunity  time  scale  during  simulation  eliminates  the  possibility  of 
working  with  real  signal  processors  on  signals  picked  up  from  radiant- 
energy  detectors. 

A  number  of  modern  analog  computers  have  frequency  responses 
that  allow  the  investigation  of  the  operation  of  modulators  without 
introducing  a  time  scale.  For  example,  type  EMU-10  electronic  simu¬ 
lation  equipment  can  be  used  to  obtain  sinusoidal  voltages  at  fre¬ 
quencies  on  the  order  of  a  few  thousand  hertz.  Furthermore,  the  ele¬ 
mentary  logic  operations  that  analog  computers  are  capable  of 
performing  make  it  possible  to  construct  model  modulators  without 
additional  sets  of  equipment  such  as  the  NGPK  [low-frequency  periodic 
oscillator] ,  the  use  of  which  always  entails  additional  difficulties 
in  trigger  synchronization. 

Electronic  and  electromechanical  function  generators  make  it 
possible  to  obtain  models  of  very  simple  modulators  relatively  easily. 
For  example,  if  a  function  corresponding  to  the  signals  picked  up 
from  a  modulator  of  the  "uniform  grating"  type  is  superposed  on  the 
field  of  an  electromechanical  function  generator,  signals  corresponding 
to  modulators  of  the  "uniform  grating"  and  "nonuniform  grating"  types 
can  be  obtained  from  the  generator's  output  by  varying  the  value  of 


the  control  voltage  according  to  a  prespecified  law.  If  the  control 
voltage  obeys  a  linear  law,  the  repetition  period  of  rectangular  volt 
age  pulses  that  simulate  the  transmitting  sectors  of  a  modulating 
disk  will  be  constant.  When  use  is  made  of  a  control  voltage  whose 
variation  is  nonlinear,  the  voltage-pulse  repetition  period  at  the 
output  of  the  fvinction  generator  will  not  be  constant  and  will  be 
determined  by  the  law  that  governs  the  control  voltage. 

Exact  assembly  of  the  transmission  function  of  a  modulating  disk 
cannot  be  accomplished,  even  when  the  simplest  raster  pattern  is  used 
because  no  function  generator  makes  it  possible  to  create  a  pulse 
edge  close  to  a  rectangular  one.  A  way  out  of  this  situation  can 
be  found  if  a  sinusoidal  voltage  at  the  fundamental  of  the  maximum 
amplitude  is  produced  in  the  analog  computer  by  using  equation  (2.24) 
and  a  circuit  after  Fig.  2.4  (for  most  analog  computers,  the  peak 
voltage  amplitude  is  100  V) .  By  using  output-voltage  limiters  it 
is  possible  to  obtain  rectangular  pulses  in  practice. 

A  sweep-frequency  generator  circuit  must  be  used  to  model  disks 
of  the  "nonuniform  grating"  type.  This  circuit  can  be  constructed 
on  the  basis  of  equation  (2.24)  and  the  circuit  in  Fig.  2.4.  The 
only  difference  is  that  the  law  governing  the  variation  of  the  coef¬ 
ficient  a  must  be  calculated  according  to  the  law  u)=(jj(t),  and  this 
law  must  be  implemented  in  a  function  generator.  The  quantity  a  (see 
Fig.  2.4)  is  represented  by  potentiometer  1;  therefore,  a  function 
generator  must  be  installed  in  place  of  the  potentiometer.  Since 
a  change  in  a  will  lead  to  a  change  in  the  amplitude  of  the  output 
voltage,  we  will  need  still  another  function  generator  to  ensure  a 
constant  amplitude  of  the  sinusoidal  voltage. 

In  many  cases  the  law  governing  the  radiant  flux  need  not  be 
reproduced  exactly  by  the  modulator.  For  example,  if  selective  con¬ 
verters  operating  on  signals  received  from  a  radiant-energy  detector 
are  tuned  to  the  fundamental  frequency,  then  there  is  no  need  to 
achieve  an  exact  match  between  the  form  of  the  simulating  signal 
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and  the  theoretical  form.  It  suffices  to  use  only  the  first  harmonic 
or  a  signal  close  to  it  in  form. 

2.3.  Simulation  of  Optical  Systems  for  Measuring  the  Angular 

Coordinates  of  Radiating  Objects 

The  simulation  of  optical  systems  for  measuring  the  angular  co¬ 
ordinates  of  radiating  objects  can  be  carried  out  by  two  methods. 

The  first  consists  in  drawin  up  a  general  equation  for  a  system  and 
then  simulating  it.  In  the  second  method,  the  equations  of  individual 
elements  in  the  system  are  put  together,  a  model  is  constructed  for 
each  element,  and  the  resultant  schemes  are  switched  according  to  a 
block  diagram. 

Each  of  these  simulation  methods  has  its  own  advantages  and  draw¬ 
backs.  The  first  method  makes  it  possible  to  get  by  with  a  minimal 
number  of  decision  elements  in  the  analog  computer,  but  at  the  same 
time  requires  correction  of  the  values  of  many  coefficients  of  the 
model  when  even  one  parameter  of  the  real  device  changes.  Although 
less  economical,  the  second  method  does  makes  it  possible  easily  to 
compensate  for  the  characteristics  of  individual  elements  and  of  an 
entire  device,  and  simplifies  the  analysis  and  synthesis  of  the  equip¬ 
ment  as  a  whole.  Simulation  using  the  second  method  actually  leads  to 
an  electronic  breadboard  model  of  the  system  and  may  be  called  elec¬ 
tronic  breadboarding. 

The  need  to  introduce  a  time  scale  due  to  the  limited  frequency 
responses  of  analog  computers  arises  in  the  construction  of  a  bread¬ 
board  or  a  circuit  for  solving  the  overall  equation  of  a  system.  The 
introduction  of  a  time  scale  limits  the  use  of  actual  selective  nodes 
of  a  system, but  all  transients,  like  steady  states,  can  not  only  be 
view  with  a  cathode-ray  display,  but  also  recorded  by  means  of  mirror- 
galvanometer  oscillographs  or  recorders. 

Before  undertaking  simulation,  we  must  determine  the  range  of 
questions  that  must  be  answered  by  the  investigation,  for  the  circuit 
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depends  on  this.  In  many  cases  it  is  advisable  to  linearize  the  equa¬ 
tions  to  be  simulated.  When  linear  equations  are  implemented,  the 
need  vanishes  for  strict  scaling.  This  greatly  expands  the  possi¬ 
bility  of  varying  the  parameters - 

In  simulating  angular-coordinate  measuring  systems,  one  must 
have  the  block  diagram  of  the  system  and  the  transfer  functions  of 
its  individual  elements  or  the  transfer  function  of  the  entire  system. 

The  modeling  of  high-frequency  filters  of  tuned  amplifiers  and 
narrow-pa ssband  filters  is  based  on  the  analogy  of  a  parallel  tuned 
circuit.  The  equation  that  describes  the  processes  which  occur  in 
the  circuit  has  the  form; 

pf 0 + A/f 0  ■{“  /*  ~  f 0 = I  -5 ) 

P 

where  p — the  derivative  symbol; 

and  Eq — the  input  and  output  signals; 
ky — the  gain; 

Af — the  circuit's  passband; 
f — the  circuit’s  natural  frequency. 

The  machine  time  t  =m. t  is  determined  by  the  time  scale  m. ,  which 

Mu  u 

varies  over  a  broad  range.  For  modeling  optical  tracking  systems, 

2  4 

m^  may  have  values  of  10  -10  . 

With  allowance  for  the  time  scale,  the  equation  of  a  parallel 
tuned  circuit  is 

— +  — £o  =  — 

m/  "’t  p  nit 

where  Eq  and  E^^  are  the  output  and  input  voltages. 

A  block  diagram  of  the  solution  of  equation  (2.26)  by  an  analog 
computer  is  presented  in  Fig.  2.6.  The  generally  accepted  procedure 
presented  in  Ref.  23,  in  particular,  is  used  to  calculate  the 


Fig.  2.6.  Block  diagram  of  Fig.  2.7.  Twin-T  filter. 

the  modeling  of  a  tuned 

amplifier. 

constants,  which  can  be  determined  by  using  potentiometer- type  voltage 
dividers,  and  to  calculate  the  transfer  ratios  for  each  input  of 
amplifiers . 

Twin-T  RC  filters  are  widely  used  in  optoelectronic  devices. 

The  transfer  function  of  the  filte.  ?hown  in  Fig.  2.7  has  the  form 

^  ( »)=s  ^  +  T 1  (7*23  4-  T^)  + 

TiTjTi/^  -t  +  Ts)  +  3^2  (rj  +  3*13  +  rj))p*+  ' 

+  (7*1  +  ris+r23  +  r3)  1 

(2.27) 

where  T^=Rj^Cj^,  T2=R2C2»  ^23“^2^3’ 

To  obtain  the  best  selective  properties  in  a  filter,  the  rela¬ 
tions  between  the  parameters  of  branches  of  the  bridge  are  selected 
as  follows: 

R2=R3=R,  C^=C3=C,  R^=|,  C2=2C 

The  nominal  values  of  R  and  C  are  adopted  on  the  basis  of  the  condi¬ 
tion  RC=l/w.  With  allowance  for  the  foregoing,  the  equation  for  a 
twin-T  filter  reduces  to  the  form 

+  (2.28) 

A  block  diagram  of  the  modeling  of  equation  (2.28)  is  shown  in 
Fig.  2.8. 


Fig.  2-8.  Block  diagram  of 
the  modeling  of  a  twin-T 
filter. 


Low-frequency  filters  usually  contain  only  a  capacitor  and  a 
resistor  and  are  first-order  inertial  elements.  The  transfer  func¬ 
tion  of  a  If  filter  is 


W(p)=i^  (2.29) 

where  T=RC. 

Circuits  of  diode-type  level  limiters  or  the  diodes  found  in 
most  machines  can  be  used  as  detectors  in  analog  computers.  The  block 
diagram  of  the  modeling  of  a  detector  and  a  If  filter  is  presented 
in  F ig .  2.9. 


Fig.  2.9.  Block  diagram  of  the 
modeling  of  a  detector  and  a  If 
filter. 


Delay  lines  or  constant-delay  circuits  can  be  implemented  in 
special  devices  called  time-delay  units.  The  relation  between  the 
input  and  output  quantities  of  a  time-delay  unit  is  given,  in  the 
ideal  case,  by  the  expression 


where  Tq  is  the  delay  time. 

The  output  quantity  of  a  delay  unit  usually  exactly  duplicates 
the  input  quantity,  but  with  a  time-shift  of  T^.  Recurrent  networks 
consisting  of  passive  and  active  two-terminal-pair  networks  that  sim¬ 
ulate  long  lines  are  used  as  time-delay  units.  Passive  networks  using 
RC  elements  have  an  extremely  limited  range  of  delay  times  and  do 
not  allow  continuous  variation  of  the  delay  time.  Many  analog  com¬ 
puters  are  equipped  with  special  delay  devices  based  on  capacitors. 

The  present  value  of  the  output  quantity  is  continuously  stored  in 
a  capacitor- type  memory,  and  is  continuously  pulled  from  it  after 
a  specified  delay  time  Tq.  The  capacitors  are  switched  by  means  of 
stepping  switches. 

It  is  advisable  to  use  magnetic  recording  for  delays  of  over 
50  s.  When  the  magnetic  recording  principle  is  used,  the  need  arises 
for  frequency  modulation  of  the  recorded  signal  and  subsequent  de¬ 
modulation.  This  complicates  the  modeling  circuit. 


Fig.  2.10.  Block  diagram  of  the 
modeling  of  an  adjustable  delay. 


A  general-purpose  delay  circuit  is  shown  in  Fig.  2.10  [10].  In 
drawing  up  this  circuit,  use  was  made  of  a  transfer  function  that 
reproduces  the  first  two  terms  of  the  expansion  of  e'^'^O  into  a  frac¬ 
tional  Fade  series; 


g-pr,  ^  yp'  ip)=  1  — 


+  6pro  +  12 


(2. 30) 
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For  the  circuit  shown,  the  transfer  ratios  of  the  eunplifiers 

12  6 

are  related  to  the  delay  time  by  the  equations 

Xq  ig 

12 

^21^11^41^“t^'  ^21^12~^22^13  ‘  time  delay  Tg  can  be  changed  by 

changing  the  transfer  ratios  of  the  amplifiers  at  different  inputs. 
The  circuit  makes  it  possible  to  vary  the  delay  time  from  0.05  to 
100  s. 


A  circuit  that  puts  out  a  sinusoidal  voltage  at  a  frequency 
is  used  to  simulate  a  reference  generator.  To  do  this,  we  construct 
a  circuit  for  solving  the  equation  y‘’+ay=0  (see  Fig.  2.4).  If  two 
reference  voltages  phase-shifted  by  an  angle  11/2  must  be  produced, 
then  an  additional  inverting  amplifier  with  a  transfer  ratio  l/co  is 
used.  A  diagram  of  the  simulation  of  a  reference  generator  that  puts 
out  two  voltages  phase-shifted  by  an  angle  tt/2  is  shown  in  Fig.  2.11. 


Fig.  2.11.  Block  diagram  of  the 
modeling  of  a  reference  generator. 


Phase  detectors  convert  the  difference  of  two  phases  of  harmonic 
oscillations  into  a  voltage.  The  output  voltage  of  a  phase  detector 
usually  is  expressed  by  a  harmonic  function  of  the  phase  difference 
having  a  maximum  value  proportional  to  the  amplitude  of  one  of  the 
input  signals: 


£o=£„-i/oSin  0. 


(2.31) 


One  multiplier  unit  is  needed  to  implement  equation  (2.31)  on  an  ana¬ 
log  computer.  In  most  computer  electronic  multipliers  operate  on 
the  equation 

z=-10^xy 

where  10^  is  the  transfer  ratio  of  the  multiplier  unit. 

A  block  diagram  of  the  modeling  of  a  phase  detector  and  filter 
is  presented  in  Fig.  2.12. 

U^CQS(Ulot«^6) 


Fig.  2.12.  Block  diagram 
of  the  modeling  of  a 
phase  detector  and  filter. 

General  recommendations  for  modeling  the  operation  of  modulators 
were  given  in  the  previous  section.  As  an  illustration,  let  us  con¬ 
sider  a  model  of  an  optical  tracking  system  with  the  modulating  disk 
shown  in  Fig.  2.1b.  The  purpose  of  the  simulation  is  to  determine 
the  effect  of  the  fabrication  imprecision  of  the  modulating  disk  on 
the  output  signals  of  the  tracking  system. 

The  initial  data  for  modeling  are; 

— the  number  of  pairs  of  sectors  of  the  modulating  disk  is  z=6 ; 

— the  rotational  velocity  is  f2=452  rad/s; 

— the  tuned  eimplifier  is  tuned  to  the  first  harmonic  of  the  car¬ 
rier  frequency; 


Fig,  2.13.  Block  diagram  of  a  model 
of  an  optical  tracking  system. 

Key:  (a)  limiter 

— the  passband  of  the  tuned  amplifier  is  A(jj=0.15w; 

— the  radiant-energy  detector  and  the  preamplifier  are  no-lag 
elements; 

— the  error  signal  amplifier  is  equipped  with  a  twin-T  RC  filter. 


Since  u)=2zQ=2x6x452=5420  rad/s,  the  passband  of  the  tuned  ampli¬ 
fier  is  Af=2Aco~1630  rad/s.  The  time  scale  m^=10^  was  introduced  during 
breadboarding  of  the  optical  tracking  system. 


The  model  of  the  modulating  disk  is  created  according  to  formula 
(2.11);  it  is  shown  in  Fig.  2.13.  The  original  relation  for  generating 


Ti(t)  is  realized  at  amplifiers  1,  2,  and  3.  When  t=0 ,  a  quantity 

equal  to  100  V  is  specified  in  the  initial  conditions  at  the  output 
of  amplifier  2.  The  amplitude  of  the  sinusoid  at  the  output  of  ampli¬ 
fier  3  will  be  ^^^m,  .  The  value  of  the  transfer  ratio  of  poten- 

u>  t  5.42 

tiometer  A^  is  ai=0.294,  while  the  transfer  ratios  of  amplifiers  1 
and  2  are  10.  This  ensures  that  the  specified  frequency  oj  is  produced. 

Together  with  potentiometer  A2,  amplifier  4  ensures  the  genera¬ 
tion  of  a  sinusoidal  voltage  with  a  100-V  amplitude  (a2=0.542,  and 
the  transfer  ratio  of  amplifier  4  is  10) .  A  level  limiter  is  used 
to  produce  near-rectangular  signals.  The  limiting  level  is  assumed 
to  be  0.2,  i.e.,  the  reference  voltages  of  amplifiers  6  and  7  respec¬ 
tively  are  -20  and  +20  V.  Amplifier  5  prevents  the  multiplier  unit 
from  having  any  effect  on  the  operation  of  the  level  limiter. 

The  function  ti  (t) ,  which  is  obtained  by  means  of  a  relay  cir¬ 
cuit,  is  generated  at  the  output  of  amplifier  5.  Up  to  a  time 

'^O^t 

t=-^ — ,  the  machine  unit,  i.e.,  100  V,  is  supplied  to  the  input  of 
^  T 

0 

multiplier  unit  Bill;  the  relay  is  triggered  at  a  time  t=— and  sup¬ 
plies  a  zero  potential  to  the  input  of  the  multiplier.  The  actuation 
time  of  the  relay  is  set  by  varying  the  voltage  supplied  to  the  input 
of  the  integrator  that  controls  the  relay  winding.  The  voltage  at 
the  input  of  integrator  8  is  set  by  means  of  potentiometer 

The  product  of  the  functions  Ti(t)  and  12 (t)  is  generated  at 
the  output  of  multiplier  unit  Bill.  The  function  Talt)  is  generated 
by  potentiometer-type  divider  A^  with  a  transfer  ratio  of  0.2.  A 
voltage  corresponding  to  Ti(t)T2(t)  is  supplied  to  one  of  the  inputs 
of  adding  amplifier  9,  while  a  voltage  corresponding  to  13 (t)  is  sup¬ 
plied  to  the  other  input.  The  function  t (t) =t 1 (t) t 2 (t) +t 3 ( t)  is  pro¬ 
duced  at  the  input  of  amplifier  9. 

A  relay  circuit  consisting  of  two  control  relays  and  P^2 
used  to  introduce  the  function  Ax  (t) .  The  actuation  time  of  each 


relay  can  be  set  by  means  of  time  counters  (some  analog  computers 
contain  control  relays  whose  actuation  time  is  set  in  time  counters) 
or  by  means  of  integrators.  The  time  when  Axtt)  appears  is  determined 
by  the  actuation  of  relay  Py2'  time  when  it  disappears  is  de¬ 
termined  by  the  actuation  of  relay  Limiter  II  prevents  the  pos¬ 

sible  appearance  of  an  elevated  voltage  due  to  inaccurate  setting 
of  the  actuation  time  of  relays  ^yl  *  limiting  level  is 

assumed  to  be  40  V  and  is  formed  by  means  of  potentiometer  (a=0.4) 

D 

and  decoupling  amplifier .  10 .  The  signal  at  the  output  of  limiter  II 
corresponds  to  the  function  t (t) =t ( t) +At (t) . 

The  radiant-energy  detector  is  represented  by  amplifier  11,  which 
changes  the  phase  of  the  input  signal  by  180°.  The  signal  supplied 
from  potentiometer  (a=0.4)  provides  a  zero  level. 

The  circuit  for  simulating  a  tuned  amplifier  (see  Fig.  2.6)  uses 
four  operational  amplifiers.  If  we  allow  for  the  digital  data  and 
time  scale  adopted,  the  equation  for  the  amplifier  has  the  form 

/Jfa-f- I,fi3£«4- 29.4  — (2.321 

p 

Potentiometer  Ag  (a=0.294)  regulates  the  tuning  frequency  uu, 
and  potentiometer  Ag  (a=0.163)  ensures  the  required  passband.  The 
transfer  ratio  of  integrator  14  is  100  (the  transfer  ratio  for  the 
input  of  10  and  the  ratio  10  of  the  feedback  circuit  [sic]).  Level 
limiter  III  with  limiting  levels  0  and  +10  V  was  used  as  the  detector. 

The  frequency  responses  of  a  low-frequency  filter  (a  block  dia¬ 
gram  is  presented  in  Fig.  2.9)  are  determined  by  potentiometer  A^^^^ 
(a=0.904),  while  the  gain  is  determined  by  the  values  of  the  transfer 
ratios  of  decoupling  amplifier  16  and  potentiometer  A^g  (a=0.452). 

An  error-signal  cunplifier  with  a  twin-T  RC  filter  is  modeled 
after  the  scheme  in  Fig.  2.8.  The  settings  of  the  transfer  ratios 
of  potentiometers  A^  are  summarized  in  the  following  table: 


■4i2  j 

'^13 

I 

0/ 

0,296 

0.321 

1 

0,272 

0.295  ' 

0.321 

0.544 

The  transfer  ratios  of  the  potentiometers  and  the  operational  ampli¬ 
fiers  were  calculated  with  consideration  for  the  time  scale  m^  and 
the  frequency  Q . 

The  phase  detectors ‘are  two  multiplier  units  to  which  a  voltage 
is  fed  from  the  T  filter  and  the  reference  generator.  Integrators 
with  a  transfer  ratio  of  0.5  were  used  as  the  filters,  while  a  cir¬ 
cuit  for  producing  two  harmonic  voltages  phase-shifted  by  an  angle 
tt/2  (see  Fig.  2.11)  was  used  as  the  reference  generator.  The  transfer 
ratio  of  potentiometer  is  a=0.204,  corresponding  to  Q/m2=0.452 
rad/s.  Amplifier  27  and  potentiometer  A^g  (a=0.442)  ensure  that  the 
quantity  cos  (f2t+i|/)  is  generated.  Potentiometers  P2  and  are  used 
to  introduce  the  initial  conditions;  the  voltage  across  P2  must  cor¬ 
respond  to  SI  cos  Ip,  while  that  across  Pg  must  correspond  to  sin 
Potentiometers  P2  and  Pg  phase  the  reference  generator. 

The  simulation  and  investigation  of  an  optical  tracking  system 
are  based  on  the  output  signals  of  each  channel  (the  system  remains 
open)  successively.  First  a  system  with  an  ideal  modulating  disk 
is  studied.  Phasing  of  the  reference  generator  ensures  that  a  signal 
corresponding  to  the  position  of  the  radiation  source  in  one  of  the 
measuring  planes  is  generated  at  the  output  of  the  tracking  system. 

In  this  case  the  voltage  from  one  of  the  integrators  (25  or  26)  cor¬ 
responds  to  a  unit  error  (mismatch)  signal,  while  the  voltage  from 
the  other  integrator  is  zero.  The  variation  of  the  signals  at  the 
output  of  the  terminal  integrators  (filters)  characterizes  the  effect 
of  the  precision  of  fabrication  of  the  modulating  disk  on  the  opera¬ 
tion  of  the  optical  tracking  system. 

The  availability  of  a  breadboard  of  an  optical  system  makes  it 
possible  tc  r.tilize  relatively  simply  the  methods  of  the  theory  of 
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accuracy  of  computing  devices  to  investigate  the  accuracy  of  the  op¬ 
eration  of  a  system  as  a  whole  and  to  synthesize  such  a  system,  since 
the  physical  processes  under  study  and  the  processes  that  occur  in 
a  breadboard  are  described  by  the  same  equations.  The  theory  of  the 
accuracy  of  computers  manipulates  the  following  basic  concepts; 

— An  Ideal  Device — a  device  that  can  reproduce  a  specified  rela¬ 
tion  between  output  and  input  quantities  with  absolute  accuracy. 

— A  Real  Device — any  actual  device  that  is  designed  to  implement 
a  specified  relation. 

— The  Error  of  a  Computing  Device — the  difference  between  the 
output  quantities  of  a  real  device  and  an  ideal  device  when  the  input 
quantities  have  absolutely  precise  values.  The  total  error  due  to 
the  error  of  the  computing  device  and  to  errors  in  the  input  of  the 
initial  data  is  called  the  total  error  in  the  output  quantity. 

Device  errors  due  to  the  limited  accuracy  of  the  input  of  ini¬ 
tial  data  and  to  errors  associated  with  the  limited  precision  of  fab¬ 
rication  and  operation  of  components  are  called  primary  errors. 

If  we  are  given  an  equation  of  the  following  type 

y=f(x^)  (2.33) 

where  y  is  the  sought  function  and  (i=l,  2,  3,  ...»  n)  is  the  argu¬ 
ment  of  the  sought  function,  and  the  equation  is  to  be  used  for  solu¬ 
tion  by  means  of  a  computing  device,  then  the  structure  of  the  comput¬ 
ing  device  is  selected  so  that  the  processes  that  occur  within  it 
are  described  by  an  equation  analogous  to  the  specified  equation. 

A  breadboard  model  of  an  equation  or  a  system  of  equations  is  pre¬ 
cisely  such  a  device. 


The  equation  of  a  computing  device  is  represented  as 


- 


i 

I 


where  — independent  variables  of  the  computing  device  which  repre¬ 
sent  X .  ; 

1 

— constants  of  the  computing  device. 

The  input  quantities  cannot  be  determined  exactly.  Moreover, 
the  constants  q^  also  cannot  be  sustained  exactly  in  the  implementa¬ 
tion  of  computing  devices. 

For  a  real  computing  device,  equation  (2.34)  may  be  written  as 

V  \  ;  A'/ ' 

where  ip'  is  the  value  of  the  function  ip  that  corresponds  to  the  actual 
and  q^,  and  Ai(/^  and  Aq^  are  the  primary  errors  of  the  computing 
device.  The  difference  between  the  functions  and  ip'  by  definition 
represents  the  error  of  the  computing  device: 


(2.36) 

Equation  (2.36)  is  unsuitable  for  practical  use  because  (p'  is 
a  nonlinear  function  of  a  large  number  of  primary  errors.  This  leads 
to  a  large  volume  of  computations.  The  calculations  use  an  approxi¬ 
mate  formula  for  AV’ obtained  by  linearizing  equation  (2.36)  with  re¬ 
spect  to  the  primary  errors  Ai|;^  and  Aq^.  Let  us  write  the  expression 
for  the  error  Ay7  in  expanded  form: 

?'(<!•/  + At/,  9*+A9i)-?(t/.  ^i).  (2.37) 

Since  At^  and  Aq  are  small  compared  with  and  q  ,  by  expanding 

X  S  IS 

the  function  tp'  q^+Aq^)  into  a  Taylor  series  and  confining 

ourselves  to  the  zero  and  first  terms  of  the  expansion,  we  can  write 
the  equation  for  a  real  computing  device  as  follows: 

M 

f’  (t/+  At,.  <7,+  A?,)=<p(t„  7,)  +  V  (S:)  At,  + 

tn 

+  A7,.  (2.38) 
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Hence,  the  error  of  the  computing  device  is 


where  f^)  ,  f^)  are  the  partial  derivatives  of  the  fiinction  v'  into 

Vd+Jo  U?  Jo 

which  the  exact  values  of  ip .  and  q  are  inserted. 

1  s 

The  primary  errors  and  partial  derivatives  must  be  known  in  order 
to  use  formula  (2.39)  for  calculations.  In  most  cases  the  primary 
errors  are  random  quantities  or  random  functions  of  time.  As  far 
as  the  partial  derivatives  are  concerned,  the  presence  of  a  breadboard 
makes  it  possible  to  seek  them  out  extremely  simply.  A  finite  incre¬ 
ment  of  the  i-th  independent  variable  or  of  the  s-th  parameter  must 
be  specified  to  determine  the  partial  derivatives  from  the  independent 
variables  and  the  constants  q^.  The  increment  in  the  function  (p’ 
that  is  due  to  a  change  in  or  and  that  is  referred  to  the  value 
of  the  increment  in  the  variable  or  parameter  will  be  the  approximate 
value  of  the  partial  derivative: 


Knowing  the  partial  derivative  makes  it  possible  numerically  to 
estimate  the  degree  of  influence  of  a  change  in  design  parameters 
and  independent  variables  on  an  output  quantity.  Furthermore,  if 
the  partial  derivatives  and  the  static  characteristics  of  the  indepen¬ 
dent  variables  and  design  parameters  are  known,  then  it  is  possible 
to  obtain  probabilistic  characteristics  of  automatic  devices  without 
constructing  a  probabilistic  model  of  the  system,  by  investigating 
a  deterministic  model  alone — in  this  case,  the  breadboard. 
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CHAPTER  3 


SPATIAL  FREQUENCY  CHARACTERISTICS  OF  MODULATORS 

3.1.  Analytical  Determination  of  the  Spatial  Frequency  Transfer 
Functions  of  Modulators 

In  using  the  spatial  frequency  method  to  investigate  the  trans¬ 
mission  of  information  through  modulators,  we  will  assume  that  they 
are  linear  components.  The  linearity  condition  presupposes  that  the 
reaction  of  a  given  element  to  the  sum  of  two  two-dimensional  signals 
is  equal  to  the  sum  of  the  two-dimensional  reactions  to  each  signal 
separately.  For  example,  if  the  input  signal  has  the  form 

F(y,  z)=aF^(y,  z)+bF2(y,  z)  (3.1) 

then  the  output  signal  for  all  values  of  a  and  b,  including  zero, 
will  be 


F(n,  C)=aFj(n/  C)+bF2(n,  Z)  (3.2) 

The  frequency-shift  theorem  also  can  be  applied  to  spatial  fil¬ 
ters.  For  a  known  weighting  function  E.  (y,  z)  of  a  "zero"  element 

0  0 

located  at  the  origin  and  a  transparent  element  E^^ (y,  z)  shifted 
by  amovints  Ay  and  Az  respectively  relative  to  the  "zero"  element, 
and  provided  that 


E^  (y,  z)=Eg^(y+Ay,  z+Az)  (3.3) 

the  Fourier  transform  will  be  defined  by  the  expression 

W,  (/li,  /v)=  j’j  f  e.  {y,  Z)  dy  dz.  (3. 4) 

After  replacing  the  variables 

u=y+Ay,  du=dy,  v=z+Az,  dz=dv 


equation  (3.4)  reduces  to  the  form 


/v)x  [2ii(|ia;/+vA2)1.  (3. 5) 


where 


Vi'oih.  <'’)  =  2)*[-2.t  {^y-^vz)]dtidz, 

x{x)—e‘*. 


The  exponential  term  in  formula  (3.5)  defines  the  spatial  shift  of 

the  weighting  function  El  (y,  z)  relative  to  the  function  E,  (y,  z) . 

0  1  0  0 

Rotation  of  the  weighting  function  in  image  space  leads  to  an 
analogous  rotation  of  the  transform  of  this  function  in  frequency 
space . 

Let  the  following  weighting  function  be  given  in  polar 
coordinates: 

where  ps/ya+z^  and  V’=arctg  z/y  respectively  are  the 
of  the  vector  p;  and  a  is  an  angle  that  defines  the 
transparent  element  having  a  weighting  function  E^ 

"zero"  element. 

The  Fourier  transform  of  the  function  E.  will 

0  1 

where  0=/y a+va  and  i|^=arctg  v/y  respectively  are  the 
of  the  spatial  frequency  vector  9. 

These  properties,  which  can  be  used  in  the  analysis  of  electronic 
filters,  facilitate  the  process  of  analytical  determination  of  the 
spatial  frequency  transfer  functions  [n^IW*]  of  spatial  filters. 

The  spatial  frequency  transfer  function  of  a  modulating  disk, 
like  that  of  any  multidimensional  filter,  is  determined  by  the  direct 
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(3.6) 

mcdulus  and  phase 
location  of  a 
relative  to  the 

be  defined  as 
(3.7) 

modulus  and  phase 


Fourier  transform  of  the  disk's  weighting  function,  which  in  this 
case  corresponds  to  its  transmission  function  T(y,  z) .  In  general, 
the  spatial  frequency  transfer  function  of  modulating  disks  can  be 
determined  directly  by  calculating  the  two-dimensional  Fourier  trans¬ 
forms  of  the  transmission  function  T{y,  z)  of  the  entire  disk: 

m 

/v)=jj  T({/,  2)it[  — +  (3,8) 

I  in  a  cartesian  coordinate  system,  and 

i 

I 

'!»)=  ?)*[  — 2n6ecos(<p— (3.9) 


in  a  polar  coordinate  system. 

To  do  this  we  must  know  the  function  r  (y ,  z)  or  t(p,  (/>)  ,  which 
depend  on  the  transmission  coefficient  and  the  size,  shape,  number, 
and  relative  location  of  the  transparent  and  nontransparent  elements 
of  the  modulating  disk. 


In  the  general  case,  the  function  T(y,  z)  varies  according  to  a 
specific  law  along  the  y  and  z  coordinates;  however,  it  usually  is 
a  step  function  of  the  coordinates  and  is  expressed  by  the  relation 


T(y,  z)  = 


1  in  the  region  of  transparent  elements; 
0  in  an  opaque  region; 

1/2  in  a  semitransparent  region. 


Therefore,  a  periodic  transmission  function  is  represented  by  the 
Fourier  series 

AT  M 

+  (3.10) 

*-0  m-0 


1 

0  0 
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. . 


where 


Here  L  and  L  are  the  repetition  periods  (spacings)  of  transparent 

y  2 

elements  of  the  modulating  disk  along  the  y  and  z  coordinates, 
respectively . 

In  many  cases  the  solutions  of  integrals  (3.8)  and  (3.9)  turn 
out  to  be  cumbersome,  even  when  one  harmonic  component  is  determined; 
furthermore,  this  solution  often  is  expressed  in  terms  of  special 
functions  (cylindrical  functions,  gamma  functions,  etc.)  which  are 
inconvenient  to  deal  with  in  spectral  spatial  emalysis. 

A  method  of  determining  the  spatial  frequency  transfer  function 
by  using  the  transform  properties  considered  above  is  the  most  con¬ 
venient  for  modulating  disks  that  have  a  periodic  structure,  i.e., 
when  the  transparent  element  are  periodically  repeated  on  the  disk 
surface.  In  this  case  we  determine  the  Fourier  transform  of  the 
transmission  function  of  one  elementary  transparent  region,  and  then, 
using  the  frequency- shift  theorem  and  the  superposition  principle, 
we  determine  the  spatial  frequency  transfer  function  of  the  entire 
modulating  disk.  For  example,  if  the  Fourier  transform  of  the  trans¬ 
mission  function  To(y,  z)  of  a  transparent  element  located  at  the 
origin  yOz  is  WQ(ip,  iv)  and  the  other  transparent  elements  are  sep¬ 
arated  from  the  "zero"  element  by  nAy  and  mAz,  respectively  (where 

n=l,  2,  3,...,  N;  m=l,  2,  3,  M) ,  then  the  Fourier  transform  of 

the  transmission  function  n  of  the  m-th  element  will  be 

/v)*[2n(iinAf/-|-vmA3'l.  (3.  1 1) 

According  to  the  superposition  principle,  the  spatial  frequency  trans¬ 
fer  function  of  the  entire  modulating  disk  will  be  given  by  the 
expression 


,V— I  .H— I 


w  (i(i.  /v)=vr,(/(i,  h)  2  2 

m^O 


X  [2n  1). 

1 3.  121 

I 


Let  us  consider  several  examples  of  the  determination  of  the 
spatial  frequency  transfer  functions  of  different  types  of  modulating 


Fig.  3.1.  Modulating 
disk  with  alternating 
transparent  and  opaque 
bands  or  strips. 


disks,  assuming  that  they  are  linear  elements  with  known  transmission 
f\inctions  T(y,  2)  or  t(p,  <f>)  . 


A  modulating  disk  that  consists  of  alternating  equidistant  trans¬ 
parent  and  opaque  strips  (Fig.  3.1)  has  a  transmission  function 


|2!<- 


0 

i-any  ,,_any 
value  value 


(3. 13) 


within  a  single  transparent  strip. 


The  Fourier  transform  of  the  transmission  function  of  an  axisym- 
metric  transparent  element  located  at  the  origin  yOz  will  be  defined  as 


A  • 

T  T 


/v)=:  f  J»[— 

"A  ♦ 

-TT 

_ gin  -tM  sin 


(3. 14) 


Expression  (3.14)  is  the  modulus  of  the  complex  transfer  function 
of  the  given  transparent  element,  i.e.. 


I 

T 


(3. 
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Each  successive  transparent  element  of  the  disk  is  shifted  by 
a  distance  2b  along  the  Oz  axis  relative  to  the  preceding  element. 
Based  on  relation  (3.12),  we  find  the  spatial  frequency  transfer  func¬ 
tion  of  the  modulating  disk  as  a  whole; 

M-l 

/v)  =  «/„(!(»,  /v)  V,(2nvr«2i>= 

-rx[2nv(iVI-l)2ftl>.  (3.16) 

where  the  expression  in  brackets  is  a  geometric  progression  the  sum 
of  whose  terms  is 


1— x[2nv(Al— l>2tl 
1  —  x(2«v2i) 


X  [  -  2«v*  (M  -  1)1  -  [2.ivt  (.Vf  -  1)1  x[2nvt(,Vf-l)] 

T  ( —  invt)  —  X  (^2n»b)  X  (2nv4) 

—  J'"  « 12nv^(.U  —  2)].  (3.1"' 

sin  2«v»  i  '  I 


With  consideration  for  (3.14)  and  (3.17),  the  spatial  frequency 
transfer  function  of  the  modulating  disk  is 


w (/[i,  /v)=  9, 


1 


Gotnyi 


itiyMb  n(x/( 

x(2nv6(/M-2)I. 


(3.  18) 


where  oo=MAb  [sic]  is  the  area  of  the  transparent  portion  of  the  mod¬ 
ulating  disk,  and  the  exponential  factor  indicates  that  the  origin 
Oyz  does  not  coincide  with  the  center  of  the  modulating  disk. 

One  quadrant  of  the  modulus  of  the  spatial  frequency  transfer 
function  of  a  modulating  disk  with  alternating  transparent  and  opaque 
band  is  shown  in  Fig.  3.2.  It  follows  from  an  examination  of  this 
figure  that  in  space- frequency  coordinates  the  maxima  of  the  disk's 
spatial  frequency  response  lie  along  the  Ov  coordinate  axis  at 
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Fig.  3.2.  Modulus  of  the  spatial  frequency 
transfer  function  of  a  modulating  disk  with 
alternating  transparent  and  opaque  bands 
(one  quadrant) . 


Fig.  3.3.  One-dimen¬ 
sional  cut  of  the  mod¬ 
ulus  of  the  spatial 
frequency  transfer 
function  of  a  modulat¬ 
ing  disk  with  alter¬ 
nating  transparent  and 
opaque  bands  along  the 
Ov  coordinate  axis. 

v=0  and  and  respectively  are  equal  to  |w()j,  'J)|y_^_Q=oo  and 

|W(u,  v=l/ (2b)  illustrated  by  Fig.  3.3. 


The  maximum  at  the  origin  of  the  frequency  coordinates  charac 
terizes  the  sensitivity  of  the  transfer  function  to  a  uniform 


distribution  of  the  radiant  flux  over  the  entire  plane  of  the  modu¬ 
lating  disk.  The  other  maxima  (at  u=0  and  v=±-^)  determine  the  sen¬ 
sitivity  of  the  transfer  function  to  a  radiant-flux  distribution  such 
that  the  edge  of  the  image  of  the  thermal  field  runs  parallel  to  the 
strips  on  the  disk. 

In  order  to  improve  the  spatial  discrimination  of  extended  back¬ 
grounds  with  sharply  defined  edges,  the  transfer  function  of  the  modu¬ 
lating  disk  must  not  have  maxima  along  the  Oy  or  Ov  frequency  axes, 
but  maxima  that  are  shifted  relative  to  them  by  some  phase  angle. 

This  requirement  is  met  by  a  modulating  disk  with  a  checkerboard  ar¬ 
rangement  of  transparent  and  opaque  squares,  the  dimensions  of  each 
of  which  are  commensurate  with  the  diameter  of  the  optical  system's 
circle  of  confusion. 


Fig.  3.4.  Modulating 
disk  with  a  checker¬ 
board  pattern  of  equal¬ 
sized  transparent  and 
and  opaque  elements. 


For  a  modulating  disk  with  a  checkerboard  pattern  of  equal- sized 
transparent  and  opaque  elements  (Fig.  3.4),  the  transmission  function 
of  a  "zero"  transparent  element  is 

1  when  !  y  1  <  Iz  I  < 

^)=  '  . 

0  when 

The  spatial  frequency  transfer  function  of  the  element  is  a  two- 
dimensional  Fourier  transform: 
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NL 
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T  T 

a  b 

“T*T 

(3.19) 

lyta  nv6 


and  is  the  modulus  of  the  spatial  frequency  transfer  function  of  this 
element. 

Transparent  elements  are  shifted  by  distances  2mb,  where  m*0, 

1,  2,  M,  along  the  Oz  axis  relative  to  the  zero  element.  There¬ 

fore,  based  on  the  frequency  shift  theorem,  the  Fourier  transform 
of  the  transmission  function  of  a  row  of  transparent  elements  with 
center  coordinates  y=0  and  z=2mb  will  be  defined  as 


W,(t,x,  M^Mba 

2RvAfi  njut 


} 


cos 


«(2«v6(.Vf—  Hj. 


(3.  20) 


by  analogy  with  relation  (3.18). 


Transparent  elements  of  the  second  row  are  shifted  by  an  amount 
b  along  the  Oz  axis  relative  to  the  first  row,  and  by  an  amount  a  rel 
ative  to  the  Oy  axis.  The  Fourier  transform  of  the  transmission  func 
tion  of  two  rows  of  a  modulating  disk  will  be  written  as 


^v) = vr,(/i»,  «v)  { 1  -1-  X  (2n  (na  -f  v6  ij ) .  ( 3. 2 1 ) 


Let  us  transform  the  expression  in  brackets: 

I X  [2«  (|»a + '’*)]=*  I”  ^ 

x[  —  npia  +  ^*)1  -(•  X [w(Ha  +  y»)l  2_ 

X  2 

= 2  cos  +  ''*)!•  (3.22) 


The  spatial  frequency  transfer  function  of  two  rows  of  a  modulating 
disk  is 


I 


/v)=2Af6a  «”?^^***  iLWx 
’  2nyM6  ttfut 

^  «t"t^+"v»(2^Vl  + 1)1.  (3.23) 

The  first  two  rows  of  transparent  and  opaque  elements  of  the  modu¬ 
lating  disk  repeat  along  the  Oy  axis  with  a  period  2a.  Hence,  the 
spatial  frequency  transfer  function  of  the  modulating  disk  as  a  whole 
will  be 


W(t\x,  /v)=Uf',(/,i,  lv){l-f.*(2nf.2a)+«(2n|*4a)+ 

-f  •  .  .  +*l2n|»(iV— l)2a]}.  (3.24) 

After  inserting  the  value  of  the  function  W2(iy/  iv)  and  perform¬ 
ing  some  simple  transformations,  we  get  the  final  form  of  the  spatial 
frequency  transfer  function  of  a  modulating  disk  with  a  checkerboard 
pattern  of  transparent  and  opaque  sectors: 


X 


W{t^,  /v)  =  9. 


sin  ‘2i>(kNa 
TxfkNa 


sin  2nv  Mb 
2nvM( 


Cos  n  (ya  +  v>) 
cos  nfta  ■  cos  nvA 


*  ( .1  {|i(  2 . V  1 )  a  -f  V  ( 2-M  1 )  ft) ) ,  ( 3. 25 ) 


where  ao=2MNab  is  the  area  of  the  transparent  portion  of  the  modulat¬ 
ing  disk. 


One  quadrant  of  the  modulus  of  the  spatial  frequency  transfer 
function  of  the  modulating  disk  in  question  is  shown  in  Fig.  3.5. 
The  maxima  of  the  spatial  frequency  transfer  function  correspond  to 
the  spatial  frequencies  y=0 ,  v=0 ,  and  y=±-^,  and  have  values 

of  Oa  and  0.405ao.  The  normalized  amplitude  spatial- frequency  re¬ 
sponse  is  shown  in  Fig.  3.6  in  a  cross  section  with  respect  to  the 
Ov  axis,  which  is  at  an  angle  a=arctg  a/b  to  the  Oy  axis. 


As  for  the  previous  modulating  disk,  the  transfer  function  of 
a  disk  with  a  checkerboard  pattern  of  transparent  and  opaque  squares 
is  very  sensitive  to  a  uniform  distribution  of  the  radiant  flux  over 
the  entire  area  of  the  modulating  disk.  The  presence  of  maxima  at 
y=±-j^  and  indicates  the  good  sensitivity  of  the  transfer 
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Fig.  3.5.  Modulus  of  the  spatial  fre¬ 
quency  transfer  function  of  a  modulat¬ 
ing  disk  with  a  checkerboard  pattern 
of  equal- sized  transparent  and  opaque 
elements  (one  quadrant) . 


Fig.  3.6.  One-dimen¬ 
sional  cut  of  the 
modulus  of  the  spa¬ 
tial  frequency  trans¬ 
fer  function  of  a 
modulating  disk  with 
a  checkerboard  pattern 
of  tremsparent  and 
opaque  elements  along 
the  O-d-  coordinate  axis , 
directed  at  an  angle 
of  45°  to  the  Oy  and 
Ov  axes. 

function  to  sharply  defined  edges  of  the  background  formation,  which 

are  positioned  at  an  angle  a»arctg  a/b  to  the  axis  of  the  image  plane. 


Fig.  3.7.  A  specially 
shaped  modulating 
disk. 


If  the  modulating  disk  considered  above  scans  along  the  Oy  or 
Oz  axis,  then  the  incidence  of  a  radiant  flux  from  an  extended  back¬ 
ground  with  clearly  defined  edges  parallel  to  these  axes  does  not 
lead  to  modulation  of  this  flux.  This  cannot  be  said  of  a  modulating 
disk  with  alternating  transparent  and  opaque  bands. 

The  modulating  disk  shown  in  Fig.  3.7  can  be  represented  as  a 
set  of  N  modulating  disks  with  transparent  and  opaque  bands.  These 
bands  are  expanded  by  an  angle  a a  with  respect  to  each  other  and  are 
at  the  following  distance  from  the  center  of  the  disk: 

where  y^^  and  Zq^  are  the  coordinates  of  the  coordinate  centers  of 
the  n-th  set  of  bands. 

In  Oy  z„  coordinates,  the  transfer  function  of  the  n-th  set  of 
n  n 

bands  will  be  given  by  the  expression 

In  the  Oyz  coordinate  system,  the  spatial  frequency  transfer  function 
of  the  same  n-th  set  of  bands  is 

(a  27) 

where 

(»' = j»  sin  na^  v' = v  cos  «o,. 

According  to  the  superposition  principle,  the  spatial  frequency  of 
the  entire  modulating  disk  is  defined  by  the  sum 

W(l^,  /v)*2  ^v)  =  V  wr^disin  >ie,+vcos«a,)X 

a-t 

X>(— 2x/?(i»sin /10,+vcosAo,)].  (a  28) 


For  N=6  the  spatial  frequency  transfer  function  of  the  modulating 
disk  is  expressed  by  a  modulus  and  has  the  form 


W'  (/[i.  /v)=  2aj  (  cos  -t-  2 


sin  nvMi 
f  3nvM6  ^ 


"X'M  .  -X-r. 

2 _  cos  y  SnvJi  cos  n[i/? 

—  (U  cos— nvS 


(3. 29) 


where  ao=MbA  is  the  total  transparent  area  of  one  n-th  set  of  bands. 


For  odd  N,  the  spatial  frequency  transfer  function  of  the  modulating 
disk  is  a  complex  function.  For  exaunple,  for  N=7  we  have 


where 


L  «A(‘ 


2  /»»"  2ainA<tv  sin  Cin^  cos  2ai.'i;?v 
\  2a]nA(iv  e\nA^  cosajiiiv 
Sin  7a%nMby  sin  e^nA^  cos  2a^/iv 
2a^AU^  CjnAit  cos  a^itiv 


-  x(-2rinAfiA)-t- 

(-2fjnAet*)  + 


sin  2asniVlfrv  sin  csWiV  cos  QasnJiv  /  r>  d 
QagnMiv  csn^  cosasnfrv  ^^s^  Wj* 

(3..30) 


a,=.sin-^; 

oo 

fj=cos-^; 

a,=sin  a,; 

r,=coso,; 

a,=sin-^o^ 

f,=cos-|-a,. 

In  the  general  case,  for  any  odd  value  of  N  the  spatial  frequency 
transfer  function  of  the  modulating  disk  in  question  is  given  by  the 
relation 


r  (/^ /v)=»,  , (- 2n^)+ 


+ 


sin  2a^vAf»  sin  nfte^A  cos  2nva„fi 


*  (  —  j. 

(3.31) 


where 
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fUl 


=  sin-^;  f,=cos^;  «=1.2 . N. 


N 


U 


The  modulating  disk  considered  differs  favorably  from  a  "check¬ 
erboard”  -type  modulating  disk  in  that  it  has  better  spatial  filter¬ 
ing.  On  this  disk  individual  modulating  bands  are  inclined  toward 
the  radii  that  pass  through  the  centers  at  different  angles;  there¬ 
fore,  each  band  cuts  off  the  target  image  at  an  angle  somewhat  dif¬ 
ferent  from  its  neighbor.  When  a  straight  edge  of  the  background 
is  incident  on  the  plane  of  the  disk,  just  one  band  in  each  group 
can  turn  out  to  be  parallel  to  this  edge,  and  the  other  bands  will 
intersect  it  an  inclined  angles.  This  enhances  the  process  of  spa¬ 
tial  discrimination,  since  the  degree  of  background  suppression  is 
directly  proportional  to  the  angular  deflection  of  the  bands  from 
the  corresponding  radial  lines. 

For  circular  modulating  disks  the  analytical  representation  and 
transformations  of  the  original  transparency  functions  are  carried 
out  in  a  polar  coordinate  system.  Since  the  image  plane  has  coordi¬ 
nates  p  and  yj,  which  are  related  to  the  cartesian  coordinates  by  the 
equations 

(3.32) 

and  the  Fourier  transform  of  the  transmission  function  is  expressed 
in  terms  of  0  and  tf,  so  that 


t*aretg  — , 


(a  33) 


where  6  and  ij;  are  the  modulus  and  phase  of  the  spatial  frequency  vec¬ 
tor  9,  respectively. 

As  an  example,  let  us  determine  the  spatial  frequency  transfer 
function  of  a  modulating  disk  with  alternating  transparent  and  opaque 
sectors  (see  Fig.  l.le).  The  transmission  function  of  the  "zero" 
transparent  sector  is  written  as 
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to(o.9)=; 


1  \9\<f:  Q<^. 

0  l9\>f:  0>^. 


(3.34) 


where  R — the  radius  of  the  modulating  disk; 

— the  central  angle  of  a  transparent  sector. 

The  Fourier  transform  of  the  transmission  function  of  a  trans¬ 
parent  sector  having  a  central  angle  ^Q='n’/N  has  the  form 


\r,(/0.<|.)= 


I  s 


x[  — 2n8ecos(9— <!»)]d;p.  (3.35) 


We  know  [2]  that  the  function  k [-z  sin  a]  is  the  generating  function 
for  the  Bessel  functions: 


*(-*sin  a)  -./^*)+22/|i{*)cos2*o- 

m 

-2i^/^^(^)sin  (2k- 1)  a. 


(3.36) 


where  I,  (z)  is  a  k-th-order  Bessel  function  of  the  first  kind, 
k 


For  our  example. 


t|  — 2n6lp cos (f —!;»)]==*  ^  ill  j 

Jit 

-  /, (2n9«)+2  ^  /*,(2.ieB)  cos  2* 

m 

-21  J /*>j(2a«B)sin(2*-l)(i+f-<^j.  (3.37) 


Inserting  expression  (3.37)  into  formula  (3.35)  and  integrating 
over  (p,  we  obtain 


g 


+4 
(-1)* 


^  cos  2*t*5in  A»,  j  r 


*»(2n8e)erfe- 


-  *  -1JiS-c08(2*~  l)«;.sin  (2*- 1)^  f  I^(2v»Q)9d^. 


i 


(3.38) 


The  integrals  in  formula  (3.38)  reduce  to  the  tabular  values 
[2]  through  substitution  of  the  variable  p=Rx. 


/«(2«ee)Qde=s/?*  ^'/■^{2n6Rx)xdx  — 

=  [2*(2ne/?)  /^2jie/?)-(-Sft»_,(2ne/?)~ 

-  2n6;?/a»_,  (2jie;?)  5,,^  (2ne/?)  +  2  t*)  1.  (3. 39) 

^  \*/  J 

A  1 

7,=^  /,t_i(2n6e)Qdo=^  J  /^i(2nORx)xdx= 

[(2*-  l)(2ae/?)/^(2a9/?)+5o.»-,(2ne/?)-. 


-  2ae/?/^(2a8/?)  {2ne/?) + 2 


, 

(-t) 

,  rj 

(3.40) 


where  S  (x)  are  Lommel's  functions,  and  F (y)  are  the  gamma  functions. 
P  /<3 


The  rotation  of  a  function  in  image  space  is  known  to  lead  to 
an  analogous  rotation  of  the  transform  function  in  frequency  space. 
Therefore,  in  determining  the  Fourier  transform  of  the  tranmission 
function  of  the  n-th  transparent  sector  from  the  origin,  the  argu¬ 
ments  of  the  cosines  of  expression  (3.38)  must  take  into  considera¬ 
tion  the  phase  shift  ,  where  n=0,  1,  2,  ...,  N-1  and  N  is  the 
number  of  transparent  sectors.  The  quantity  is  the  angular  repe¬ 
tition  period  of  transparent  sectors  and  is  measured  between  radii 
drawn  through  the  centers  of  adjacent  sectors. 


The  Fourier  transform  of  any  n-th  transparent  sector  has  the 

form: 


-97- 


«^,(/8.  y  f-^^sin  >t,cos2A(1.- 

ICt 

sin  (2*-  l)'^cos(2*-l)(’l-2/*<Po)A  ]. 

(a  41) 

where  the  integrals  and  J2  are  determines  from  formulas  (3.39) 
and  (3.40),  respectively. 

After  applying  the  superposition  principle,  we  obtain  the  spatial 
frequency  transfer  function  of  a  modulating  disk  with  alternating 
transparent  and  opaque  sectors; 


v-i 

r  (id,  =  (10, «  =  + 

fl— 0 

'S  *90 tos  2k  l<!»- (.V  - 1 )  90I  h - 

K  S{n  4^fpo 

-<  (2*_l)^cos{2^-l)l-;.-(.V-l)?„]X 


^  »m(2*--l)2Airto  y 

sin(2ft-l)2fo  *' 


(3.42) 


The  first  term  of  expression  (3.42)  characterizes  the  sensitiv¬ 
ity  of  the  transfer  function  to  a  uniform  distribution  of  the  radiant 
flux  over  the  entire  area  of  the  modulating  disk  and  does  not  have 
a  significant  effect  on  the  transmission  of  small  signals.  The  most 
characteristic  features  for  analysis  of  a  transfer  function  are  its 
maxima,  which  are  at  spatial  frequencies  0  that  are  remote  from  the 
center  of  the  frequency  coordinates  and  that  are  determined  by  the 
second  and  third  terms  of  formula  (3.42). 


As  we  see  from  the  calculations  presented  above,  the  process 
of  determining  spatial  frequency  transfer  functions  entails  consid¬ 
erable  difficulties,  even  for  modulating  disks  that  have  transparent 
elements  of  comparatively  simple  configuration.  Furthermore,  the 
expressions  obtained  for  the  spatial  frequency  transfer  function  are 
not  very  convenient  when  used  for  further  spectral  analysis  of  the 


4 


i 
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transmission  of  information  through  optoelectronic  devices.  It  will 
be  shown  below  that  the  location  of  the  maxima  of  the  spatial  fre¬ 
quency  transfer  function  in  frequency  space  must  be  known  in  analyz¬ 
ing  the  transmission  of  information  on  radiating  objects  through  mod¬ 
ulating  disks.  It  is  difficult  to  obtain  such  data  from  the  theore¬ 
tical  expressions  of  the  spatial  frequency  transfer  functions  of  most 
modulating  disks.  Therefore,  experimental  methods  are  used  to  inves¬ 
tigate  the  spatial  frequency  transfer  functions  of  modulating  disks 
having  elements  of  complex  configuration. 

3.2.  An  Experimental  Method  of  Determining  the  Spatial  Frequency 

Transfer  Functions  of  Modulators 

According  to  the  classical  theory  of  radiation,  a  light  field 
can  be  described  by  using  the  vectors  of  electric  field  strength 
S(r,  t)  and  magnetic  field  strength  B(r,  t) .  In  experimental  studies 
of  a  light  field  the  instruments  react  to  the  electric  field.  A  Fou¬ 
rier  series  or  integral  is  used  to  identify  the  four-dimensional  equa 
tion  for  electric  field  strength,  and  this  strength  is  treated  as 
the  sum  of  two  complex  components; 

«  (F,  (F.<)+i?'-(r,0.  (3.43) 

Both  components,  which  are  complex-conjugate  to  each  other 

»-(7, (r.  4 

contain  equivalent  physical  information  and  are  called  the  complex 
strength  or  the  complex  amplitude  of  the  light  wave.  These  defini¬ 
tions  will  be  used  below  in  examining  diffraction  effects  when  modu¬ 
lating  disks  are  irradiated  with  a  parallel  beam  from  a  monochromatic 
source  of  coherent  radiation. 

Let  us  imagine  that  the  modulating  disk  is  a  very  simple  dia¬ 
phragm  with  an  aperture  having  an  area  to  and  that  it  is  located  in 
a  plane  (Fig.  3.8).  When  a  plane  monochromatic  wave  from  a  point 
source  is  incident  on  the  modulating  disk,  a  diffraction  pattern  can 
be  observed  on  a  screen  located  in  plane. 02* 
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Fig.  3.8.  Toward  calculating  the  dis¬ 
tribution  of  field  strength  in  a  dif¬ 
fraction  pattern. 

Kirchoff's  formula  [11]  is  used  to  Calculate  the  complex  strength 
g(x,  y,  z,  t)  in  plane  Q2  for  known  strength  of  the  light  wave  inci¬ 
dent  on  the  modulating  disk; 


*(-*/•)  dgg.ii.c.f)  1  . 

r  dn  I  ’ 


(3.44) 


where  k=2iT/\ — the  wave  number; 

n — the  normal  to  the  plane  of  the  modulating  disk; 
r — the  distance  between  points  A  and  P. 

This  formula  makes  it  possible  to  calculate  the  electromagnetic  field 
at  the  observation  point  P  when  the  field  of  the  light  wave  is  speci¬ 
fied  on  an  arbitrary  surface. 

Integral  (3.44)  is  taken  within  the  area  of  the  transparent  aper¬ 
ture  on  the  assumption  that  the  strength  behind  the  opaque  portion  of 
the  modulating  disk  is  zero.  It  can  be  seen  from  Pig.  3.8  that  the 
quantity  r  is  defined  by  the  expression 


r*=(jc  -  5)»-h(j>- (z  -  0*. 


where 
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In  view  of  the  fact  that  observation  is  done  at  a  considerable 
distance  from  the  modulating  disk,  we  may  assume 


—  r,)  (r  -  r,) « ‘_>r  {r-r„). 

Then 

r=r.-^±-^^  -f- -+  *''-+  .  (3.45) 

ro  'ir^ 

K  plane  light  wave  incident  on  the  plane  of  the  modulating  disk 
has  the  complex  amplitude 

«(5.  n,  C.  C)*[-(u./-*5)l.  i3.46) 

and  its  derivative  with  respect  to  the  direction  of  the  normal  n  in 
the  plane  of  the  aperture,  i.e.,  for  5=0,  is  defined  as 


rd»(5.ii.C.O 

L  4/1 

6-0  [  45 

=  C)x(  — (1)0. 

(3.47) 


The  other  partial  derivative  in  Kirchoff's  formula  (3.44)  is 
expressed  by  the  relation 


d  r»(  — fer)  -j  _ 

I  r  Jt-o 

-)-gC  +  n 

rp  2ro  l\ 

t-o 

(3.48) 

For  small  diffraction  angles  9  and  large  r  that  are  much  greater 
than  the  dimensions  of  the  modulating  disk  (which  always  is  the  case 
in  practice),  the  quantity  r  in  the  denominator  of  expression  (3.48) 
may  be  replaced  by  rQ.  Then  derivative  (3.48)  will  be  written  as 


d  (»(-*/•)  ■ 

-k(r  - 

»*l  +  1 

(A# 

t-o  L  M" 

■  ..  ■I 

2/-0  )J' 

(3.49) 


and 
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!!i±iL  <il, 

2ro  To  ' 

where  d  is  the  dicuneter,  or  the  linear  dimension  of  the  aperture  of 
the  modulating  disk. 


If  the  condition 


is  observed,  then  the  last  term  in  expression  (3.49)  may  be  ignored; 


d  f  X  ( —  *r) 
dc  L 


^)\ 


(3.50) 


With  allowance  for  expressions  (3.47)  and  (3.50),  the  complex 
amplitude  of  the  light  field  intensity  in  the  plane  of  the  diffraction 
pattern  will  be  given  by  the  relation 


2nr 


n. 


(Jo. 


’(3.51) 


In  experimental  research  it  may  be  assumed  that  the  light  field 
incident  on  the  modulating  disk  is  steady.  Steadiness  does  not  at 
all  mean  that  the  field  is  time -independent;  on  the  contrary,  the 
field  usually  oscillates  extremely  rapidly.  This  means  that  our  data 
on  the  field  do  not  change  with  time.  In  this  case  by  steadiness  we 
mean  equal  probability  of  the  statistical  description  of  the  light 
beam  with  respect  to  a  shift  in  time. 


When  a  monochromatic  light  source  is  used,  so  that  k=27r/Xo=const , 
the  expression  in  front  of  the  integral  is  independent  of  the  inde¬ 
pendent  variables  y  and  z,  and  since  we  are  only  interested  in  the 
spatial  distribution  of  the  complex  cunplitude  of  the  light  field's 
intensity,  we  will  rewrite  formula  (3.51)  as 

if  ( y,  z ) = M  ^  (T).  g  *  f  - f n  -  ; )  1  t/n d'. 

J  I  \  >0  ^0  ^0  ^0  '  J 

«• 

(3.n.’ 
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where 


hfo 

Let  us  convert  to  the  new  independent  variables  by  substitution: 
^  =v .  We  will  divide  expression  .(3.52)  by  A,  and  we  find 

that  the  normalized  complex  amplitude  of  a  light  wave  in  the  plane 
of  the  diffraction  pattern  is  the  analog  of  the  spatial  frequency 
transfer  fxinction  of  a  very  simple  modulating  disk: 


/v)=W'(/ii,  /v)  = 

= (3.53) 

In  formula  (3.53)  the  limits  of  integration  extend  to  infinity  on 
the  assumption  that  the  function  ^(n»  is  equal  to  zero 

outside  the  area  of  an  elementary  transparent  element. 

The  modulating  disks  used  in  optoelectronic  devices  have  a  more 
complicated  arrangement  of  a  large  number  of  transparent  and  opaque 
sectors.  However,  when  a  coherent  radiation  source  is  used  as  the 
irradiator,  each  transparent  element  may  be  represented  as  a  separate 
coherent  source. 

When  elementary  waves  are  added,  a  time-stable  pattern  of  illum¬ 
inance  distribution  that  has  the  form  of  alternating  maxima  and  minima 
is  produced  in  the  plane  of  the  diffraction  pattern.  Thus,  when  a 
wavefront  emanates  from  a  coherent  radiation  source,  the  diffraction 
effects  represent  the  interference  of  the  elementary  waves  originating 
from  each  elementary  sector  of  the  surface  of  the  original  wave  [sic] . 

To  prove  the  validity  of  using  multiray  interference  with  wave- 
front  division  to  obtain  the  spatial  frequency  transfer  functions 
of  modulating  disks,  we  will  show  how  the  light  field  is  distributed 
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Fig.  3.9.  Modulating  disk 
with  diffracting  elements 
that  are  identically 
oriented  in  the  plane  of 
the  disk. 

in  the  interference  pattern  if  the  modulating  disk  has  on  its  surface 
N  transparent  diffracting  elements  that  are  identically  oriented  in 
the  plane  of  the  disk  (Fig.  3.9). 

The  coordinates  of  some  diffracting  point  of  the  n-th  transparent 
element  are  defined  as 

y=y,-]r^,  (3.54) 

The  complex  amplitude  of  the  light  wave  in  the  diffraction  pat¬ 
tern  from  the  n-th  transparent  element  is,  according  to  (3.53), 

m 

/v)=JJ«{-2n(p(y,-}-T))-l-v(z,4-i:)lWT)<<2.  (3.55) 

Factors  containing  the  quantities  y^  and  may  be  removed  from 
the  integral  sign,  and,  as  a  result  of  the  interference  effect,  the 
complex  amplitude  of  the  light  wave  is  equal  to  the  sum  of  the  complex 
amplitudes  obtaining  upon  diffraction  from  each  of  the  N  elements 
of  the  modulating  disk: 

9  (ii»,  /v) = 2  *  I  ~  -f-  vz,)l  X 

4*1 

m 

xJJ«l-2n(!»n+>(C)ldtjdC.  (3.56) 


Thus,  the  distribution  of  the  complex  amplitude  of  a  light  wave  is 
defined  as  the  product  of  two  factors: 

tv),  (3. 57) 

one  of  which,  ^  (iy,  iv) ,  depends  solely  on  the  shape  and  size  of 
the  n-th  individual  transparent  element,  while  the  other,  F(iy,  iv) 
characterizes  the* relative  position  of  the  transparent  elements  of 
the  modulating  disk  with  respect  to  each  other,  i.e.,  it  reflects 
the  frequency-shift  and  superposition  principles.  Hence,  the  use 
of  diffraction  effects  and  multibeam  interference  with  wavefront  di¬ 
vision  makes  it  possible  to  obtain  experimentally  the  spatial  fre¬ 
quency  transfer  functions  of  any  modulating  disks  by  measuring  the 
illuminance  in  the  interference  pattern  while  the  modulating  disks 
are  irradiated  with  a  monochromatic  parallel  beam  from  a  coherent 
radiation  source. 

We  know  that  the  complex  sunplitude  of  the  light  field  intensity 
cannot  be  measured  with  recording  instruments.  Photoelectric  devices 
that  measure  the  energy  parameters  of  the  field  usually  are  used  as 
meters.  In  other  words,  they  measure  the  square  of  the  modulus  of 
the  light  wave's  amplitude: 

£(l».v)=iy(fi*.HP.  (3.58) 

The  sensor  is  located  in  the  plane  of  the  interference  pattern  and 
obeys  a  specific  law  in  scanning  over  this  plane.  To  improve  resolu¬ 
tion,  during  measurements  the  area  of  the  sensor  is  made  as  small 
as  possible,  and  the  distance  r^  from  the  modulating  disk  to  the  in¬ 
terference  pattern  is  made  as  large  as  possible.  The  scale  of  the 
spatial  frequency  is  determined  by  the  value  of  rg  and  the  radiation 
wavelength  according  to  the  equation  m=l/(XorQ),  which  is  constant 
for  a  monochromatic  radiator. 

Since  the  sensor  measures  the  intensity  of  the  light  field  in 
y  and  z  coordinates  (Fig.  3.8),  we  must  multiply  the  coordinates  of 
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plane  Q2  by  a  scaling  factor  in  order  to  obtain  a  function  that  de¬ 
pends  on  the  independent  variables  u  and  v: 

V=m2 

To  eliminate  errors  in  the  measurement  of  the  peak  values  of 
the  function,  we  must  have  a  radiation  source  with  a  time-constant 
light  beam.  It  was  demonstrated  in  the  previous  section  that  if  the 
maximum  value  of  the  transmission  fvinction  of  the  modulating  disk 
is  unity,  then  the  spatial  frequency  transfer  function  has  a  maximum 
equal  to  the  total  area  of  the  disk's  transparent  elements.  In  proc¬ 
essing  the  experimental  results  it  is  best  to  deal  with  the  normalized 
spatial  frequency  transfer  function  of  the  modulating  disk,  which 
can  be  defined  as 


(3.59) 


where  E(u,  v) — the  measured  value  of  the  intensity  function  of  the 
light  field  in  the  diffraction  pattern  for  a  given 
experiment; 

E  — the  maximum  value  of  the  intensity  of  the  light  field 

maX 

in  the  given  experiment. 

A  continuous-wave  laser  with  a  time-stable  light  be6un,  a  high 
degree  of  coherence,  and  a  high  degree  of  monochromaticity  is  used 
as  the  light  source.  A  block  diagram  of  the  experimental  setup  for 
determining  the  spatial  frequency  transfer  function  of  modulating 
disks  is  shown  in  Fig.  3.10.  In  processing  the  measurements,  the 
normalized  value  of  the  transfer  function's  modulus  is  defined  as 
the  mathematical  expectation  over  discrete  values  of  the  illuminance 
at  points  of  a  plane.  An  estimate  of  the  mathematical  expectation 
is  the  arithmetic  mean  of  the  realizations  of  N  successive  measure¬ 
ments  : 
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Fig.  3.10.  Block  diagram  of  the 
experimental  setup  for  determin¬ 
ing  the  spatial  frequency  trans¬ 
fer  function  of  modulating  disks: 
(1)  radiation  source;  (2)  optical 
collimator;  (3)  modulating  disk; 
(4)  sensor;  (5)  amplifier; 

(6)  instrument;  (7)  oscillo¬ 
graph;  (8)  scanning  mechanism; 

(9)  power  supply  and  control 
unit. 


where  U^(u»  v) — the  value  of  the  voltage  of  the  instrument  during 

the  n-th  measurement  at  a  given  point  in  the  plane; 
^r^max — maximum  value  of  the  voltage  at  the  central  maxi¬ 
mum  during  the  n-th  measurement. 

We  will  use  several  examples  to  demonstrate  the  experimental 
results  and  the  degree  of  agreement  between  the  experimental  and  the¬ 
oretical  data.  If  the  modulating  disk  is  a  diaphragm  with  a  circular 
aperture,  then  the  diffraction  pattern  has  a  number  of  concentric 
circles  after  the  disk  is  irradiated  with  a  beam  of  parallel  rays 
from  a  monochromatic  coherent  radiation  source.  The  theoretical  value 
of  the  spatial  frequency  transfer  function  of  such  a  modulating  disk 
is  given  by 

(»)-*.  (3.61) 


where  oo — the  aperture  area; 

rg — the  aperture  radius; 

0 — the  modulus  of  the  spatial  frequency  vector. 


In  Fig.  3.11  the  solid  lines  indicate  the  theoretical  moduli 
and  the  dashed  lines  the  experimental  normalized  moduli  of  the  spatial 
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Fig.  3.11.  One -dimensional  cuts  of 
the  spatial  frequency  transfer 
functions  of  diaphragms  with  cir¬ 
cular  apertures;  (a)  d=»0.3  mm; 

(b)  d®0.6  mm. 

frequency  transfer  functions  of  modulating  disks  with  apertures  hav¬ 
ing  diameters  of  0.3  and  0.6  mm.  The  experimental  data  are  in  good 
agreement  with  theoretical  calculations  and  confirm  the  feasibility 
of  experimental  determination  of  the  spatial  frequency  transfer  func¬ 
tion  by  the  method  described.  Fig.  3.12  shows  the  distribution  of 
illuminance  in  the  interference  pattern  at  the  output  of  a  modulating 
disk  with  12  transparent  and  opaque  sectors. 

The  interference  pattern  at  the  output  of  a  modulating  disk  with 
a  checkerboard  pattern  of  transparent  and  opaque  square  elements  (Fig. 
3.13)  shows  the  distribution  of  maxima  in  the  space-frequency  domain. 
Shown  in  Fig.  3.14  are  the  one-dimensional  normalized  spatial  fre¬ 
quency  responses,  as  calculated  theoretically  (curve  a)  and  as  ob¬ 
tained  experimentally  (curve  b)  at  an  angle  of  45®  to  the  coordinate 
axes. 


For  circular  modulating  disks,  it  is  convenient  to  represent 
the  spatial  frequency  transfer  function  as  the  product: 


-108- 


Fig.  3.12.  Distribution  of  illum¬ 
inance  in  the  interference  pattern 
at  the  output  of  a  modulating  disk 
with  alternating  transparent  and 
opaque  sectors. 


Fig.  3.13.  Distribution  of  illuminance 
in  the  interference  pattern  at  the  out¬ 
put  of  a  modulating  disk  with  a  check¬ 
erboard  pattern  of  transparent  and 
opaque  elements. 

\r(rt,t)=vr(/9)-r  (i).  (3.62) 

where  W(i9)  and  respectively  are  the  amplitude  and  phase  spa¬ 

tial  frequency  responses. 


Thus,  for  the  modulating  disk  shown  in  Fig.  l.lj,  which  is  a 
circle  with  an  asymmetric  arrangement  (relative  to  the  center)  of 


Fig.  3.14.  One-dimensional  normalized 
amplitude  spatial  frequency  response  of 
a  modulating  disk  with  a  checkerboard 
pattern  of  transparent  and  opaque 
elements:  (a)  theoretically  calculated; 

(b)  experimental. 

of  transparent  and  opaque  elements,  the  normalized  modulus  of  the 
spatial  frequency  transfer  function  reduces  to  the  form: 

(9.  Wlw- Icos  2n9Y! exp  ( - -nfl -L J  (A6)|  U?', 

(3.  63) 

In  formula  (3.63)  the  first  term  determines  the  shape  and  width  of 
the  maximum  near  the  origin  of  frequency  coordinates,  while  the  second 
term  describes  the  amplitude,  shape,  and  width  of  the  maxima  of  the 
spatial  frequency  transfer  function  of  the  disk  at  higher  spatial 
frequencies.  The  functions  in  formula  (3.63)  can  be  approximated 
by  some  simple  relations: 

|r.(Aa)|  =  ^expl-Y(9-9,)*l,  (3.64) 

i3.65) 

where  the  parameters  y  and  A  depend  on  the  dimensions,  shape,  and 
number  of  transparent  elements;  N  characterizes  the  number  of  sectors 
on  the  checkerboard  half  of  the  modulating  disk;  and  0o  is  the  modu¬ 
lus  of  the  spatial  frequency  at  which  the  maximum  is  located. 


3.3.  Conversion  of  Random  Optical  Signals  by  Multidimensional  Filters 


Radiating  objects  having  various  linear  dimensions  ranging  from 
the  minimum  discernible  to  dimensions  commensurate  with  the  field 
of  view  are  seen  through  optoelectronic  devices.  Amidst  this  diver¬ 
sity  of  objects,  the  device  must  obtain  maximal  information  on  the 
required  object,  which  differs  from  the  other  objects  in  shape,  size, 
or  energy. 

The  task  of  filtering  is  known  [12]  to  consist  in  determining 
the  statistical  characteristics  of  a  filter's  output  signal  from  given 
statistical  characteristics  of  one  or  several  random  processes  that 
are  the  input  variables.  The  methods  of  statistical  analysis  used 
in  the  investigation  of  electronic  devices  can  be  used  to  evaluate 
multidimensional  spatial  filters. 

Let  a  random  signal  X(y,  z)  with  an  autocorrelation  function 
y*'  2')  and  a  spectral  density  S^(vi,  v)  be  supplied  to  the 

input  of  a  multidimensional  filter  with  a  weighting  function  E^(y,  z) . 
Then  the  autocorrelation  function  and  spectral  density  of  the  output 
signal  respectively  are 

Rs,  (t|.  rf  ;  c,  C')= JJ  dydz  JJ  E,{%  y;  C,  z) £, (r,',  y’; Q',  j')  x 

X  (y,  z,z’)d y'  dz',  (3. 66) 

(I*,  v) = \W  (/}»,  /v)|»  ((»,  V).  (3, 67) 

To  determine  the  dispersion  of  the  output  signal  of  a  multidi¬ 
mensional  filter,  it  suffices  to  set  n'=n  and  i;'=^  in  formula  (3.66). 
Then  we  obtain 


D,.  (ti.  C)-/?,.  (t),  g= 

m  m 

y;  c,  z)Et{r\^y'\  c,*')x 
XRj\{y,y’\  z,z^)]dy'd;d.  (3.68) 

Let  us  determine,  for  example,  the  dispersions  of  the  output 
signals  after  spatial  white  noise  passes  through  filters  having 
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bell- shaped  and  rectangular  weighting  functions.  As  we  know,  the 
spectral  density  of  spatial  white  noise  is  a  constant,  i.e., 

S  (y,  v)=D  =const.  After  the  signal  passes  through  a  multidimen- 

X  X 

sional  filter  having  a  bell-shaped  weighting  function  whose  spatial 
frequency  transfer  fxinction  has  the  form 


r  v)=exp  (3.69) 

the  spectral  density  of  the  output  signal  is 

5,.  ((»,  v)=D^exp  (3. 70) 

To  find  the  dispersion  of  the  output  signal,  we  must  integrate 
the  spectral  density; 


Dx,  =  f  J Sjc^  ((»,  y)d^d\— 


*  r  2 


But  if  the  white  noise  is  converted  by  a  filter  having  a  rec¬ 
tangular  weighting  function,  then  the  spectral  density  and  dispersion 
of  the  output  signal  respectively  are 


5„(i».  v)=-^^sin*|i  (3.72) 

m 

Q  }^£x^  Tf  — 5in*j»i^sin*  V—  d^d\=Djiy^z^, 

’  (2n)J  J  3  (i2v3  2  2  * 

(3.  73) 

i.e.,  the  larger  the  area  of  the  rectangle,  the  greater  the  disper¬ 
sion  of  the  output  signal. 

The  quality  of  a  spatial  filter  always  can  be  evaluated  if  we 
know  the  signal-to-noise  ratio  at  the  filter's  output,  i.e.,  if  we 
have  data  on  the  valid-signal  and  noise  spectra.  Let  a  valid  signal 
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E  (n»  5)  and  a  background  or  other  noise  signal  E  (ri/  C)  be  defined 
c  n 

in  the  plane  of  the  modulating  disk  in  the  form  of  spatial  functions 
of  irradiance.  In  the  general  case,  these  signals  are  random  func¬ 
tions  of  coordinates.  We  will  assume  that  these  functions  describe 
steady-state,  ergodic,  isotropic  processes  that  are  not  intercorre- 
lated.  Then  the  overall  irradiance  distribution  in  the  plane  of  the 
modulating  disk,  which  determines  the  total  input  signal  of  the 
filter,  is 

EQ(n,  5)=E^(n,  5)+Ejj(n,  O  (3.74) 

A  spatial  filter  must  isolate  the  valid  signal  E^(n,  c)  from  the  over¬ 
all  signal  E^Cn,  ^) ,  and  this  can  be  done  only  when  and  E^  have 
different  properties. 

In  analyzing  the  filtering  properties  of  modulating  disks,  we 

must  know  the  spectral  power  densities  of  the  valid  signal  E  (u,  v) 

c 

and  the  noise  signal  E^{u,  '>)  . 

The  signal-to-noise  (S/N)  ratio  often  is  used  as  a  criterion 
of  valid-signal  detection.  Hereafter  by  signal-to-noise  we  mean  the 
ratio  of  the  signal  power  from  a  legitimate  source  to  the  root-mean- 
square  value  of  the  noise  signal  power  at  the  output  of  the  modulat¬ 
ing  disk. 

The  signal  power  from  a  valid  source  at  the  output  of  a  spatial 
filter  can  be  represented  as 

"  2 

J J  5,  (i*.  y)W  (|»,  V) » 12n{it,  n + vC)l  dv  . 

(3.75) 

If  the  Khinchine-Wiener  spatial  spectrum  of  the  background-noise  power 
in  the  plane  of  the  modulating  disk  is  defined  by  the  function 
Sjj(u,  v)  ,  then  after  the  spatial  filter  this  noise  spectrvim  will  be 

v)  W  (I*.  v)|*.  (3. 76) 
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After  the  spatial  filter  the  rms  value  of  the  background-noise  power 
will  be 


Pn(%  C)=JJ 5,(|».  v)|r  (3. 77) 

In  cartesian  coordinates,  the  square  of  the  ratio  of  the  radiant 
power  of  a  valid  source  to  the  rms  noise  power  is  given  by  the  follow¬ 
ing  expression; 


^n(l.C)- 

m 

Se  (|i.  v)  IT  (i(t,  iv)  « [2n  (fti)  +  vC))  dv 


jT  S„  (|i,  v)l  IT  (i(i. 


(3.78) 


In  a  polar  coordinate  system,  the  same  formula  will  be 


R  u  2 

f  8  </e  f  Se  (fl)ir  (rt ,  4.)  » [2«8r  cos  (r  ~  4,)) 

0  b _ 

_  _  , 

fs,(8)8(<8|nr(/a.4<)|*i/4. 


(3.79) 


The  absence  of  a  dependence  of  the  functions  and  on  the 
independent  variable  ijj  is  explained  by  making  the  assumption  of  iso¬ 
tropic  radiation  from  noise  and  from  a  valid  source,  and  of  an  iso¬ 
tropic  distribution  of  the  irradiance  from  them  in  the  plane  of  the 
modulating  disk. 


The  S/N  ratio  is  a  characteristic  that  reflects  the  pott.  itial 
for  separating  a  signal  from  backgroiand.  A  comparison  of  the  filter¬ 
ing  properties  of  modulating  disks  based  on  the  S/N  criterion  makes 
it  possible  to  obtain  filtering  performance  that  is  acceptablv  in 
practice. 

The  identification  of  valid  signals  by  optoelectronic  devices 
is  known  to  take  place  against  a  radiating  background  whose  radiance 
distribution  in  space  is  a  random  function  of  coordinates,  radiation 


wavelength,  and  time.  For  many  devices  that  have  a  small  field  of 
view  and  that  operate  for  a  short  period  of  time,  we  may  assume  that 
a  fully  specific  background  formation  that  undergoes  practically  no 
change  in  shape  in  time  falls  within  the  field  of  vision. 


The  images  of  valid  and  spurious  radiators  in  the  plane  of  the 
modulating  disk  have  different  dimensions  and  shapes  with  different 
energy  levels.  The  field  of  the  background  image  brightness  is  seen 
as  a  two-dimensional  random  process  B(y,  z)  which  is  a  family  of  ran¬ 
dom  values  of  brightness  that  are  placed  in  correspondence  with  the 
nonrandom  argioments  y  and  z.  When  the  background  is  considered,  it 
is  assumed  to  be  steady.  This  is  justified  by  the  specific  features 
of  optical  devices;  by  the  small  field  of  vision  and  low  probability 
that  the  device  will  operate  in  regions  of  clear  statistical  nonuni¬ 
formity — the  horizontal  line,  the  zone  adjacent  to  the  sun,  etc.  The 
assumption  that  the  random  process  B(y,  z)  is  isotropic  presupposes 
that  its  statistical  properties  are  constant  in  all  directions.  This 
is  justified  by  the  fact  that  when  work  is  done  far  from  the  horizon 
and  sun,  there  are  no  factors  to  cause  a  change  in  the  probabilistic 
characteristics  of  the  bright  background  field  in  any  direction.  The 
following  expression  for  the  energy  spectrum  of  the  natural  background 
is  given  in  Ref.  5; 


•S'*  (I*.  v)= 


(3. 80) 


where  a,  b — quantities  that  characterize  the  anisotropy  of  the 
background; 

g — a  quantity  that  determines  the  rate  of  decrease  of  the 
spectrum  with  increasing  spatial  frequency; 

Sno=Sn(0'  0) • 


Actual  radiation  sources  that  an  optical  device  must  identify 
from  among  the  diverse  radiators  generally  are  point  sources  that 
occupy  an  insignificant  fraction  of  the  field  of  vision.  Because 
of  the  electromagnetic  nature  of  light  and  the  finite  dimensions  of 
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the  field  of  vision,  the  image  from  a  point  source  has  the  form  of 
a  diffraction  mottle  in  the  general  case.  For  an  idealized  optical 
system,  the  spatial  frequency  spectrum  of  the  image  of  a  point  source 
is  given  by  the  formula  [19]; 


5.(1*.  (3. 81) 


K.3T; 

where  r — the  radius  of  the  diffraction  mottle; 
a=-Tr* ; 

Eq — the  ir radiance  amplitude. 

For  a  fixed  wavelength  Xo 


where  Eq=Eq/Xo. 


(3. 82) 


As  the  diauneter  of  the  image  of  a  radiation  source  increases,  the 
energy  distribution  in  the  spectrum  draws  closer  and  closer  to  the 
origin  of  frequency  coordinates. 


£(•) 


Fig.  3.15.  One-dimen¬ 
sional  spatial  fre¬ 
quency  spectra,  (a)  val¬ 
id  signal;  (b)  background. 


Comparing  the  spatial  spectra  of  the  background  and  a  valid  sig¬ 
nal  (Fig.  3.15),  we  see  that  a  significant  fraction  of  the  energy 
is  distributed  at  low  spatial  frequencies  in  the  spectra.  Therefore, 


we  may  conclude  that  it  is  difficult  to  isolate  a  valid  signal  at 
low  spatial  frequencies.  It  is  best  to  determine  the  filtering  prop¬ 
erties  of  modulating  disks  at  high  frequencies. 


All  modulating  disks  are  characterized  by  the  fact  that  their 
spatial  frequency  transfer  functions  have  subordinate  maxima  at  higher- 
order  frequencies.  We  will  call  these  maxima  the  primary  maxima  and 
we  will  excimine  the  filtering  properties  of  some  modulating  disks 
in  the  region  of  the  primary  maxima.  To  do  this,  let  us  simplify 
expression  (3.78)  by  applying  the  Cauchy-Schwarz  inequality  to  the 
numerator ; 


Z*< 


JJ  Sc  (it.  v)4»rfv 


\r (i|i,  iv)dfidv 


J  |  5,  ((».v)|ir  (it,v)l2diidv 


(3.83) 


and  we  determine  all  functions  for  values  of  the  spatial  frequencies 
Uo  and  Vo  that  correspond  to  the  location  of  a  primary  maximum.  Then 
the  approximate  S/N  expression  will  be  converted  to  the  simple 
formula; 


X< 


Sc  (w>.  V 
l5«(W.vo)l‘  *  ■ 


(3.84) 


Formula  (3.84)  seemingly  does' not  take  into  account  the  charac¬ 
teristics  of  the  modulating  disk,  since  its  transfer  function  is  miss¬ 
ing.  However,  this  does  not  indicate  that  the  structure  of  the  mod¬ 
ulating  disk  does  not  affect  the  S/N  ratio.  Indeed,  the  primary 
maxima  of  the  transfer  functions  of  modulating  disks  are  located  at 
specific  spatial  frequencies  whose  values  reflect  the  arrangement, 
size,  and  shape  of  the  transparent  and  opaque  elements  of  a  given 
modulating  disk.  These  factors  also  determine  the  number  of  primary 
maxima. 


As  an  exaunple,  let  us  compare  the  filtering  properties  of  two 
modulating  disks:  one  with  alternating  transparent  and  opaque  sectors 
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(disk  I)  and  one  with  a  checkerboard  pattern  of  transparent  and 
opaque  sectors  (disk  II) . 


To  simplify  the  solution  of  the  problem,  we  will  assume  that 
the  background  is  isotropic,  i.e.,  that  its  spectrum  is 

and  that  the  signal  spectrum  is  described  by  expression  (3.82). 

We  will  assume  that  the  width  a  of  a  transparent  band  of  modulat¬ 
ing  disk  I  is  equal  to  a  side  of  a  square  of  disk  II,  and  that  the 
valid-signal  image  has  a  radius  r=a/2.  Then  the  S/N  ratio  for  modu¬ 
lating  disk  I  is 


K+oY/*. 


and  for  modulating  disk  II 


(3.86) 


^  (2jvy  i4. + vj. 

(3.87) 


where 


If  we  assume  that  (ui  ^ )  =3^^  (U2  j ,  V25),  i.e.,  that  the  peak 
values  of  the  background  spectrum  at  these  frequencies  are  equal, 
then  the  estimate  of  the  effectiveness  of  the  filtering  properties 
of  the  modulating  disks  can  be  determined  as  the  ratio 


H  /i(2yT«IH,r)  gi  y'»  1 

Xi  /i(2w)  \  l4+  O*  / 


(3.88) 


When  6*3/2,  a<<uo»  and  2r*a,  this  ratio  is  equal  to  1.163,  i.e.,  mod¬ 
ulating  disk  II  ensures  that  a  valid  signal  will  exceed  the  background 
by  approximately  16%  compared  with  modulating  disk  I. 


For  an  anisotropic  distribution  of  irradiance  in  the  image  of 
a  background  process  in.  the  plane  of  the  modulating  disk,  the  spatial 
frequency  spectrum  of  the  background  is  written  as 


v)=- 


(Ojil  +  4v*  +  o*)^ 


(3. 89) 


Let  us  examine  the  effect  of  background  anisotropy  on  the  esti¬ 
mate  of  the  effectiveness  of  modulating  disks  I  and  II.  We  will  as- 
s\ime  that  in  formula  (3,89)  a=l  and  b-*-0,  i.e.,  a  background  process 
represents  the  edge  of  a  solar-illianinated  cloud  or  the  skyline  that 
is  perpendicular  to  bands  of  modulating  disk  I.  In  this  case  the 
formula  for  estimating  the  effectiveness  is 


e^ll^iijgiiWL^const. 

Xt  /2/,(2w) 


(3.90^ 


i.e.,  for  a  specific  value  of  the  radii  of  the  image  of  a  valid  ob¬ 
ject  is  a  constant  [sic] . 


If  the  edge  of  a  cloud  or  the  skyline  is  parallel  to  transparent 
bands  of  modulating  disk  I,  i.e.,  b=l  and  a-^O,  then  the  estimate  of 
effectiveness  will  be  given  by  the  formula 


_ X2_  /i(2KW) 


K2/i(2n|ior) 


(3.91) 


If  a<<vi,  then  since  the  value  of  xi  tends  toward  zero.  The  re¬ 

sults  obtained  are  fully  consistent  with  the  physics  of  the  process 
of  modulation  of  the  light  from  radiating  objects.  Modulation  of 
the  background  is  not  observed  in  the  first  case,  where  a=l  and  b->0, 
and  the  values  of  xi  ^^d  X2  remain  constant,  depending  solely  on  the 
cunplitudes  of  the  luminous  fluxes  of  a  valid  radiator  and  the  back¬ 
ground,  The  second  case,  with  a-*-0  and  b=l ,  leads  to  a  situation  in 
which  modulating  disk  I  obtains  an  identical  percentage  modulation 
for  the  light  from  a  valid  radiator  and  the  background,  since  the 
dimensions  of  the  valid-object  and  background  images  are  commensurate. 


To  investigate  the  transmission  of  information  from  radiating 
objects  through  circular  modulating  disks,  the  Khinchine-Wiener  spec¬ 
trum  of  background  radiation  is  expressed  by  the  relation 


[•*  (a  co»*  i{i  -(-  (  sin2  i|i  4-  oS)J^ 


(3. 92) 


Let  us  consider  the  effect  of  anisotropy  of  a  background  process  on 
the  filtering  properties  of  the  modulating  disk  shown  in  Fig.  1.11, 
using  the  primary  maxima  of  its  spatial  frequency  transfer  function 
and  the  phase-frequency  response: 


N 


(3.93) 


If  the  background  process  is  isotropic,  then  the  S/N  ratio  based  on 
formula  (3.84)  will  assxame  the  form 


(3.94) 


For  extreme  anisotropy,  i.e.,  when  the  edge  of  an  illuminated 
cloud  or  the  skyline  is  located  in  the  field  of  vision  of  the  optical 
device  (a=l,  b->0)  ,  the  S/N  ratio  is 


X  — 


-Il« 


7  li~co«jy4i»rf4 
ie*(cot*4  +  o»)r 


(3.95) 


where  the  superscript  *  means  that  the  investigation  is  done  in  the 
presence  of  an  anisotropic  background. 


After  the  integrals  of  formula  (3.95)  are  calculated,  the  S/N 
ratio  for  the  extreme  case  of  anisotropy  turns  out  to  be  lower  than 
the  analogous  ratio  for  an  anisotropic  background  process  by  a  factor 
of  approximately  3.  This  indicates  that  background  anisotropy  great¬ 
ly  impairs  a  modulating  disk's  ability  to  discriminate  a  valid  signal 


3.4.  Some  Problems  of  the  Optimization  of  Spatial  Filters 

If  a  random  valid  signal  m(y,  z)  and  noise  n(y,  z)  are  supplied 
to  the  input  of  a  multidimensional  filter,  then  it  is  desirable  that 
the  filter  reproduce  the  valid  signal  with  a  minimal  rms  error  at 
its  output. 

We  will  use  the  following  conventions: 

(y,  z)=m(y,  z)+n(y,  z) — the  total  input  signal; 

g(y»  2) — the  required  output  signal,  which  is 
equal  to  the  valid  signal  m(y,  z)  in 
the  case  in  question; 

Stf,  S  ,  S  ,  S  --the  spectral  densities  of  signals  <», 
r  m  n  g 

m,  n,  and  g,  respectively; 

S  ,  S  --the  Fourier  images  of  the  cross-corre 
mn  nm 

lation  functions  of  the  valid  signal 
and  noise. 

The  spectral  densities  of  spatial  valid  signals  and  noise  that 
are  encountered  in  practice  usually  are  rational  functions  of  the 
spatial  frequencies  u  and  v.  For  example,  the  spectral  density  of 
a  valid  signal  can  be  represented  as 


v)= 


q» 

l»(<m  <v)l»  ’ 


(3. 96) 


and  the  noise  spectrum  can  be  represented  by  the  expression 


5.(P.  v)= 


le(<p.'v)r 


(3.97) 


In  the  absence  of  a  cross  correlation  between  the  valid  signal 
and  the  noise,  the  cross-correlation  functions  are  equal  to  zero; 
therefore,  their  Fourier  transforms  also  are  equal  to  zero: 

S  «S  sO 
nm  mn 
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The  solution  of  the  problem  of  an  optimal  filter  that  was  found 
by  Wiener  and  Kolmogorov  also  can  be  extended  to  spatial  filters  if 
we  introduce  the  assumption  of  steadiness  and  isotropy  of  the  radia¬ 
tion  from  background  noise  and  observed  objects  that  the  filter  iso¬ 
lates  from  the  overall  light  flux  that  falls  within  the  field  of 
vision  of  the  device's  optical  system. 

By  analogy  with  an  electronic  filter,  the  expression  for  the 
optimal  transfer  function  of  a  spatial  filter  has  the  form 


« (p,  «)  L  —  J 


(3.98) 


where  Wj^(p,  q) — the  desired  transfer  function; 
p=iu,  q=iv — operators. 

v)+S,(|»,  v).  (3.99) 

We  will  assume  that  the  filter  must  reproduce  the  input  signal 
at  the  output,  i.e.,  that  the  desired  transfer  function  of  the  filter 
is  Wj^  (iu ,  iv)  *1 . 

To  determine  the  optimal  transfer  function  of  a  spatial  filter, 
we  will  convert  the  bracketed  expression  in  formula  (3.98): 


Vi(P.  -J-) 

0(  — A  —9) 


Up,  9)A{  —  p.  —q) 


where 


=1  u+i  I— 

j  k  and  1_ — constants  that  can  be  determined  by  the 

n  n 

method  of  coefficient  comparison; 
q)  and  q) — functions  of  the  operators  p  and  q. 


In  the  theory  of  electrical  filters  the  frequency  transfer  func¬ 
tion  is  determined  uniquely  by  the  locations  of  two  singularities: 
poles  and  zeros.  For  a  physically  feasible  system,  i.e.,  for  one 
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whose  weighting  function  is  g(t)=0  when  t<0,  the  poles  of  the  trans¬ 
fer  function  are  located  in  the  left  half-plane.  Therefore,  to  de¬ 
termine  the  optimal  transfer  function  from  (3.98)  we  use  only  that 
portion  of  expression  (3.100)  which  is  marked  with  the  symbol  [ 

For  spatial  filters  the  weighting  function  Et  ly,  z)  exists  when 
y,  zSO,  therefore,  the  criterion  for  physical  feasibility  of  the  fre¬ 
quency  response  of  a  spatial  filter  is  that  the  arrangement  of  the 
poles  must  be  symmetric  with  respect  to  the  coordinate  axes.  Hence, 
expression  (3.100)  is  used  in  its  entirety  to  determine  the  optimal 
transfer  function  from  formula  (3.98).  Then  the  expression  for  the 
transfer  function  of  an  optimal  multidimensional  filter  will  have 
the  form 


'*'o*r(*l*.  iv)s 


v) 


S«0*.  »)+.Sii(m'v) 


(3. 101) 


We  will  illustrate  the  use  of  formula  (3.101)  by  solving  several 
examples. 

Let  a  valid  signal  be  supplied  to  the  input  of  an  optimal  multi¬ 
dimensional  filter.  The  spectral  density  of  the  signal  is 


v)=0,exp|— ^(l**-|-v*)j 


Let  the  noise  be  spatial  white  noise  with  a  spectral  density 

Sj^(y,  v)=D^ 

The  spatial  frequency  transfer  function  of  an  optimal  filter  is 


W'-rd*.  v)=- 


1 


o,  '■? 

1  +  -TP  e*P—  (p*  +  v2) 


(3.  102) 


When 


•^>>1  v)=l; 
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if 


■^<1  v)=-^exp|--^(,.a  +  v2)j  . 

For  real  backgrounds  and  valid  signals,  the  most  frequent  case  is 
D^<<D^,  in  which  the  optimal  spatial  frequency  transfer  function  of 
a  spatial  filter  has  the  following  form  if  the  background  noise  is 
spatial  white  noise: 

v)= ^ 

Since  this  function  is  symmetric  with  respect  to  the  origin  of  fre¬ 
quency  coordinates,  we  will  write  it  in  one-dimensional  space: 

'>^o«(»)=^exp(-a»e»).  (3.104) 

where 


The  one-dimensional  weighting  function  of  an  optimal  spatial 
filter  will  be  determined  by  the  inverse  Fourier  transform  of  its 
spatial  frequency  transfer  function; 


^  exp(— o*»*)exp(2nM»)dd.  (3.105) 


The  integral  in  formula  (3.105)  is  the  tabular  one  [2].  It  is 
equal  to: 


Yn  o. 

a  Dg 


or,  if  the  weighting  function  and  the  transmission  function  are 
equivalent. 


II 
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The  function  x  _(ri)  takes  on  its  maximum  value  when  ri=0.  It  was 
onT 

noted  above  that  the  maximum  value  of  the  transmission  function  of 
a  modulating  disk  does  not  exceed  unity.  In  the  example  considered 


Oj  ft 


whence 


Dm  ^ 

Da  iQiVn  ‘ 

Given  pi=0.4679  and  r=0.1  mm,  we  obtain 

^-0,06  «!, 

This  confirms  the  legitimacy  of  the  assumption  made  above. 


If  the  spectral  characteristics  of  the  valid  signal  and  back¬ 
ground  are  expressed  by  the  relations 


and 


i.(,,  (3.107) 

2ar  ]/|*i  +  V* 


5,(1*,  v)= 


(^  +  va  +  a*)^ 


(3. 108) 


then  the  spatial  frequency  transfer  function  of  an  optimal  spatial 
filter  is 

/V)  - - \ - - 

_ 2.7r  »  ■«»-!- _ 

Dm  ! ,  (2nr  \  i»^  +  vi)  (^i2  +  v2  4.  a2)^ 

(3.  109) 

and  it  is  a  symmetric  function  with  respect  to  the  origin  of  coordi¬ 
nates  Ou  and  Ov. 
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Let  us  write  this  function  in  a  one-dimensional  version  and  on 

the  assumption  that  D  <<D  : 

m  n 

Ij/i  2jtrv 

The  one-dimensional  transmission  function  of  an  optimal  modulating 
disk  will  be  given  by  the  expression 

(»*-l-g*)*exp(2n<n»)^ft-  i3.  Ill) 

The  transmission  function  of  a  modulating  disk  is  real;  there¬ 
fore: 


fo.T(’l)=-^^J  /i(2«r»)i^±5^cos2.m»d».  (3.112) 

We  will  designate  the  integral  in  formula  (3.112)  as  Jq  and  express 
it  as  the  sum  of  two  integrals  and  J2  on  the  assumption  that  g=l 

J  /,(2nr0)9cos2nn*'d#-j' 

«i 

+0*  J  /i(2nr9)  i cos 2.^^l9d9.  (3.  U3i 

Both  integrals  are  tabular  and  can  be  calculated  as 


-) 

V  2 

J _ 

O.Oor 

y, = 2a^x)sarcsln  -1-  =*  —  (r* — Tf  )'/*, 


3.  114i 

(3. 115) 


where  r(|)s0.89  is  the  gaunma  function  of  the  argviment 
Hence 
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Then  the  transmission  function  of  an  optimal  modulating  disk  for  in¬ 
put  signals  having  spectral  densities  (3.107)  and  (3.108)  will  be 
expressed  by  the  formula: 


(Then  n=0,  the  function  is  tom (0) 1*  i.e., 

whence 


- - - .  (3.118) 

Dm  a2(/-»  +  0.05) 

It  follows  from  expression  (3.118)  that  in  order  to  preserve 

the  conditions  t_,„=1  and  D  <<D^,  the  radius  r  of  the  image  of  a  valid 
max  m  n 

radiation  source  must  be  increased  or  the  quantity  a  must  be  larger 
than  unity.  For  exaunple,  for  r=0.1  ram  and  Dj^=0.05D^,  a=3.16. 

To  determine  the  dimensions  of  the  transparent  apertures  of  a 
modulating  disk,  we  set  expression  (3.117)  equal  to  zero: 


0.05r« 


Hence  we  obtain 


n^tr 


i.e.,  once  again  we  find  that  the  dimensions  of  the  apertures  of  the 
modulating  disk  must  be  equal  to  those  of  the  image  of  a  valid 
radiator . 
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